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ANALYSIS OF STALL FLUTTER 


OF A HELICOPTER ROTOR BLADE 

By Peter Crimi 
Avco Systems Division 


SUMMARY 


A study of rotor blade aeroelastic stability was 
carried out, using an analytic model of a two-dimensional 
airfoil undergoing dynamic stall and an elastomechanical 
representation including flapping, flapwise bending and 
torsional degrees of freedom. Results for a hovering rotor 
demonstrated that the models used are capable of reproducing 
both classical and stall flutter. The minimum rotor speed 
for the occurrence of stall flutter in hover was found to 
be determined from coupling between torsion and flapping. 
Instabilities analogous to both classical and stall flutter 
were found to occur in forward flight. However, the large 
stall-related torsional oscillations which commonly limit 
aircraft forward speed appear to be the response to rapid 
changes in aerodynamic moment which accompany stall and 
unstall, rather than the result of an aeroelastic instabil- 
ity. The severity of stall-related instabilities and re- 
sponse was found to depend to some extent on linear stabil- 
ity. Increasing linear stability lessens the susceptibility 
to stall flutter and reduces the magnitude of the torsional 
response to stall and unstall. 


INTRODUCTION 


Aeroelastic stability of a helicopter rotor blade is 
a multifaceted problem because of the extreme variations of 
the aerodynamic environment within the flight envelope of 
the aircraft. In hovering flight, a blade can undergo 
classical binary flutter (Ref. 1) or stall flutter (Ref. 2). 

In forward flight, the linear instability experienced by 
systems with periodically varying parameters can occur 
(Ref. 3). While these types of instability are not normally 
encountered with blades of current design, due to the rela- 
tively low disc loading and weak coupling of translational 
and rotational degrees of freedom, they are certainly not 
precluded from new designs, particularly those intended to 
extend present performance capabilities . Of immediate 
concern, however, in both design and operation, is the 
occurrence of large -amplitude torsional oscillations and 
excessive control-linkage loads associated with blade stall 
on the retreating side of the rotor disc at high forward 
speed or gross weight, effectively limiting aircraft per- 
formance. This problem has prompted a number of recent 
studies of dynamic stall and the effects of stall on blade 
dynamics (Refs. 4-8). 

While stall has been identified as a causal element of 
the problem, the nonlinearity of the stall process, coupled 
with the unsteady aerodynamic environment, has precluded an 
analysis to the depth required to gain a thorough under- 
standing of the mechanisms involved. In particular, it has 
not been clear whether the blade undergoes a true aero- 
elastic instability, a simple forced response, or some 
hybrid phenomenon which takes on the character of one or 
the other extreme, depending on flight conditions and blade 
vibrational characteristics. 

Stall flutter for axial flight is amenable to analysis 
by empirical methods similar to those developed for analyz- 
ing stall flutter in cascades (Ref. 9). The flutter 
mechanism for that case has been identified as deriving 
from the extraction of energy from the free stream by the 
periodic variation of the aerodynamic moment . Analogous 
methods applied to the forward-flight problem (Refs . 10 
and 11) have been inconclusive, however, the primary diffi- 
culty possibly being in applying empirical methods without 
a clear definition of the underlying mechanism of the problem. 

A method was recently developed for analyzing dynamic 
stall of an airfoil undergoing arbitrary pitching and 
plunging motions which provides an ideal tool for analyzing 
the stall problem in forward flight. The method, which is 
described in detail in Ref. 7 , employs models for each of 



the basic flow elements contributing to the unsteady stall 
of a two-dimensional airfoil. Calculations of the loading 
during transient and sinusoidal pitching motions are in 
good qualitative agreement with measured loads. Dynamic 
overshoot, or lift in excess of the maximum static value, 
as well as unstable moment variation, are in clear evidence 
in the computed results . 

This study was directed to analyzing the aeroelastic 
stability of a helicopter rotor, particularly as it relates 
to stall, using the method of Ref. 7 to compute aerodynamic 
loading. The representation of the elastomechanical system 
includes flapping and flapwise bending degrees of freedom 
as well as torsion. A listing of the computer program used 
to perform the calculations is given in Appendix A. 


4 


SYMBOLS 


blade semichord, m 

mean lift coefficient, ratio of time average 
of 1 to pfi 2 R 2 b 

lift coefficient, C-j. = C 1 /( P U 2 b) 
moment coefficient referred to quarterchord, 

C m c/4 = ra cA /(2 p U b > 

blade chord, m 

mode shape of first uncoupled torsional mode, 
unit tip deflection 

mode shape of first uncoupled flapwise 
bending mode, unit tip deflection 

tip deflection due to flapping, semichords 

tip deflection due to bending, semichords 

translational coordinates of 2-D system 
(i = 1, 2), semichords 

moment of inertia of 2-D system about pitch 
axis, kg - m 

blade moment of inertia about elastic axis 
per unit span, kg - m 

translational spring stiffnesses of 2-D 
system (i = 1, 2), N/m^ 

torsional spring stiffness of 2-D system, N/rad 

lift per unit span at aerodynamic reference 
radius, N/m 

offsets of springs from pitch axis of 2-D 
system (i = 1, 2), m 

total blade mass, kg 

blade mass per unit span, kg/m 

aerodynamic moment per unit span at aerodynamic 
reference radius, N 


masses of 2-D system, kg/m 
rotor radius, m 

inner radius of blade lifting surface, m 

aerodynamic reference radius, m 

instantaneous free-stream speed at aerodynamic 
reference section, m/sec 

reference speed, U Q = r R , m/sec 

distance aft of elastic axis of blade 
section mass center, m 

distance aft of pitch axis of mass center 
of m^, m 

generalized coordinate of 2-D system, equivalent 
to , semichords 

generalized coordinate of 2-D system, equivalent 
to h 0 , semichords 

angle of attack, deg 

flapping hinge offset, m 

collective pitch angle, deg or rad 

blade tip torsional deflection, rad 

angle of zero restraint of 2-D system torsion 
spring, rad 

advance ratio, ratio of forward speed to XL R 

free-stream density, kg/m^ 

dimensionless time, t = U Q t/b 

blade azimuth angle measured from downwind 
direction, deg or rad 

rotor rotational speed, rad/sec 

dimensionless rotor speed, ft* = ftR /( w @ b) 

flutter frequency, rad/sec 



frequency of first uncoupled, 
torsion mode, rad/sec 


nonrotating 


frequency of first uncoupled, nonrotating 
flapwise bending mode, rad/sec 



PROBLEM FORMULATION 


Aerodynamic Loading 

In the flutter analysis, only leading-edge stall was 
considered, so the following relates specifically only to 
that type, even though the basic method can treat trailing- 
edge stall as well. When the airfoil is not stalled, the 
flow elements represented are (see Figure la): (l) the 

laminar boundary layer from the stagnation point to separa- 
tion near the leading-edge, (2) the small leading-edge 
separation bubble; and, (3) a potential flow, including a 
vortex wake generated by the variation with time of the 
circulation about the airfoil. When the airfoil is stalled, 
as indicated in Figure lb, the flow elements are: (l) the 

laminar boundary layer, (2) a dead-air region extending 
from the separation point to the pressure recovery point; 
and, (3) a potential flow external to the airfoil and 
dead-air region, again including a vortex wake. The 
analytic representations of these elements are described 
briefly below. Details are given in Ref. 7 * 

Potential Flow. — Given the airfoil section character- 
istics and motions, together with the distribution of 
pressure In the dead-air region if the airfoil is stalled, 
the flow and pressure over the airfoil must be determined 
to compute the integrated load and analyze the boundary 
layer. The problem was formulated by imposing linearized 
boundary conditions of flow tangency and pressure, using 
a perturbation velocity potential derived from source and 
vortex distributions. The resulting coupled set of singular 
integral equations is solved by casting the singularity 
distributions in series form and solving for the unknown 
coefficients by imposing boundary conditions at prescribed 
points . 

Boundary Layer. — Because the relative importance of 
the individual elements of the boundary layer flow as they 
affect dynamic stall could not be established in advance, 
the representation in Ref. 7 was made as general as possible. 
The method of finite differences for unsteady flows with 
variable step size in both streamwise and normal directions, 
was employed, with the error in each finite-difference 
approximation the order of the square of the step size. 

It was determined from preliminary calculations performed 
for this study that, at least for leading-edge stall, 
results are virtually unaffected by assuming quasi-steady 
flow in the boundary layer. That assumption was therefore 
employed for all flutter computations, to take advantage 
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(1) LAMINAR BOUNDARY LAYER 

(2) DEAD-AIR REGION 

(3) AIRFOIL AND VORTEX WAKE 



(b) LEADING -EDGE STALL 
Figure 1 FLOW ELEMENTS 
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of the resulting substantial savings in computer storage 
requirements and computing time. 

Dead-Air Region . — The function of the model of the 
dead-air region is to define the streamwise distribution 
of pressure in that region, given the locations of the 
separation and recovery points and the pressure at the 
recovery point. The dead-air region is assumed to consist 
of a laminar constant -pressure free shear layer from sepa- 
ration to transition, a turbulent constant-pressure mixing 
region, and a turbulent pressure-recovery region. The 
laminar shear layer is analyzed by the method of Ref. 12, 
assuming quasi-steady flow. The turbulent mixing and 
pressure-recovery regions are analyzed using the steady-flow 
momentum integral and first moment equations. Profile 
parameters in these regions are assumed to be universal 
functions of a dimensionless streamwise coordinate, with 
those functions derived from an exact viscous-inviscid 
interaction calculation. Matching of approximate solutions 
for the mixing and pressure-recovery regions at their inter- 
face completes the analysis. 

Leading-Edge Bubble . — The leading-edge bubble on an 
unstalled airfoil is analyzed using the same basic relations 
employed for the dead-air region. Given the boundary-layer 
parameters at separation, the length of the bubble and the 
amount of pressure rise possible, for that length, in the 
pressure recovery region, are computed. That pressure 
rise is compared with the rise in pressure in the potential 
flow over the length of the bubble. If the latter is 
greater than the former, the bubble is assumed to have 
burst, and the stall process is initiated. 

Loading Calculation Procedure . — Calculations proceed 
by forward integration in time, using the blade motions 
derived by integrating the equations of motion of the 
elastomechanical system. If, at a given instant, the 
airfoil is not stalled, the potential flow is computed, 
and the boundary layer and leading-edge bubble are analyzed 
to check for bubble bursting. If the airfoil is stalled, 
the pressure distribution in the dead-air region is com- 
puted, the potential flow evaluated, and the boundary layer 
is analyzed to locate the separation point. The last two 
steps are repeated iteratively until assumed and computed 
separation points agree. Rate of growth of the dead-air 
region is determined from an estimate of the rate of fluid 
entrainment derived from the potential-flow solution. 

Unstall is determined by first postulating its occurrence 
and analyzing the leading-edge bubble which would then 
form to ascertain whether that event did in fact occur. 
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During unstall, the dead-air region is washed off the 
airfoil at the free-stream speed. 


Elastoraechanical Representation 

The equations of motion for a rotor blade with flapping, 
flapwise bending and torsional degrees of freedom can be 
written in the form (Ref. 3 ) 
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where h p and hM are tip displacements due to flapping 

and bending, respectively, in semichords, ej. is torsional 
displacement at the blade tip and the frequencies* are the 
following functions of rotational speed: 
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The inertial and centrifugal-force coefficients are 
given by 
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♦Barred quantities are dimensionless frequencies, U n /b 
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The complexity of the aerodynamic representation pre- 
cludes evaluation of the generalized forces Fp , Fd and 
P@ by the usual strip approximation. It was felt essential, 
however, to retain both translational degrees of freedom in 
the investigation of the forward-flight problem, so a 
simple two-dimensional model of the dynamics could not be 
used. Therefore, a two-dimensional airfoil suspended in 
such a way as to have three degrees of freedom was analyzed. 
Inertial and stiffness parameters were assigned to make the 
coupled natural frequencies of the two-dimensional system 
match those of the rotor blade. 
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The system analyzed is shown schematically in Figure 2. 
The matching of the two-dimensional system with the blade 
dynamics proceeds as follows. Three generalized coordinates 
are first defined to correspond to those of the blade. 
Clearly, angular displacement ©i should correspond to blade 
torsional displacement at the blade tip. The counterparts 
of flapping and bending, Z R and Zd, respectively, are 
defined by p 


Zp = A n h n + Bh c , Z,t = A c h n - Bh, 
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Figure 2 TWO-DIMENSIONAL E LASTOMECH AN ICAL SYSTEM 
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With the above definitions, Zp + Z$ = -h, , to give the 
correct translational correspondence. It can further be 
shown that the uncoupled natural frequencies of the two- 
dimensional system match those of the blade, provided 
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By comparing the generalized masses of the two systems, it 
follows that 

m i b2 A 0 " - h b2 /( M se r2 > 

r2 ) 5 . 

The last relation, together with Eqs . (1) and (2), fixes 

mg/n^ : 


m 2 / m i 


(1 + + X m 

* n n 1 

( \ m + s/r 


^ 4 ) 


- 1 




Equating the corresponding coefficients of the characteris- 
tic equations of the two systems provides three additional 
relations, which can be solved for the coupling parameters 
x, l Sl , lg 2 . That calculation is outlined in Appendix B. 
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To complete the matching, quasi-steady approximations 
to the damping terms of the flapping equations are equated 
with the result that 


m l R /(~ A i) 



R a [1 - (r 0 / R ) 4 ] 


u/o 0 



1 - ( r 0 / R ) 2 
1 - (r 0 / R)4 


pi sin ^ 


where ^ r^ = U 0 • The aerodynamic reference radius r^ was 
selected to be .75H. 

The angle of zero restraint in torsion was varied 
periodically to approximate the effects of cyclic pitch 
variation in forward flight, according to the formula 

e = e 0 [l - 2 (R/r R ) n sin*] 

This variation gives nominally constant lift . 

The equations of motion were solved by integrating 
analytically, using linear extrapolations to approximate 
the variation of lift and aerodynamic moment over the 
interval of integration. This scheme was found to give 
satisfactory results, provided the time interval of integra- 
tion is no longer than about one fifth of the period of the 
coupled mode having the highest natural frequency. 


RESULTS OP COMPUTATIONS 


Configurations Analyzed 

Vibrational and aerodynamic characteristics of the 
blade analyzed were selected to correspond to those of the 
model rotor blade described in Ref. 2. That blade is un- 
twisted, of constant chord, with offset flapping hinge. 
Pertinent dimensionless parameters of the model blade are 
listed in Table 1. 


TABLE 1 

BLADE PARAMETERS FOR NOMINAL CONFIGURATION 
Parameter Value 


b/R 

8/ R 

r o/ R 

"©</ "«f 0 

PR b S /M, 

x n/ b 

m R/M b 

*e ^ M b R 


.0435 

.0543 

.174 

3.69 

.00431 

.216 

1.055 

3.51 x 10' 4 


Two elastomechanical configurations in addition to the 
nominal one were analyzed. One of these had Wq = 2.5, 

with all other parameters as listed in Table 1. ° 0 

The third configuration had x m /b = .108, with the remaining 
parameters as listed in Table 1. 

The bending mode shape, which was computed by a 
finite-element method, was found not to vary appreciably 
over the range of rotational speeds of interest. The mode 
shape for « 0 /ft = 1.26, which is plotted in Figure 3, 
was used for all computations. The torsional mode shape 
for the nonrotating blade, also shown in Figure 3, was used 
to evaluate torsional inertia parameters. 
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The test blade had a NACA 23012 section. The 
variation of static lift and moment coefficients with angle 
of attack for this section were computed from a series of 
transient pitch calculations, and are shown in Figure 4, 
together with the measured section characteristics, from 
Ref. 13. The aerodynamic model is seen to give nearly the 
correct maximum lift, but at a slightly lower angle of 
attack, and, as indicated from the variation of C m c /i±, the 
computed center of pressure is somewhat further aft than 
that of the actual airfoil section below the stall angle. 


Stability in Hover 


Initial calculations were performed for hovering flight 
with the nominal configuration, to allow a direct comparison 
with the test results of Ref. 2. First, rotor speed was 
varied parametrically, with the collective pitch at a value 
well below the stall incidence. A classical bending-torsion 
instability was encountered at 12 * = 12 R /(«@ 0 b) = 5.3 

with 03 f / 03 Qq - .803. The variation of bending, flapping, 
and torsional displacements with azimuth angle at flutter 


onset are shown in Figure 5 . By way of comparison, tests 
(Ref. 2) yielded classical flutter at about 12* = 7.1 
with « f 7 wq 0 = . 72 . 
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It should be noted that since the system stability was 
analyzed by direct simulation, a precise point of linear 
instability was not computed. The values of 12* at onset 
of a linear instability, both for hover and forward flight, 
were obtained by successively increasing or decreasing rotor 
speed, in small steps, until the transient response changed 
from convergent to divergent, or visa versa. The maximum 
error in the value of flutter speed, for the results pre- 
sented here, is estimated to be about three percent. 


Susceptibility of the system to stall flutter was in- 
vestigated next. It was found that a torsional limit cycle, 
at approximately the highest coupled natural frequency of 
the system, could be triggered for 12 * as low as 3 . 4 . 
Computed blade motions for stall flutter at 12* of 3.5 
are shown in Figure 6 . 

For 12* below 3.4, a limit cycle could not be set up, 
regardless of the initial conditions or the collective pitch 
angle. Severe oscillations involving repeated stall and 
unstall could be made to occur by imposing a large initial 
bending deflection. However, the flapping response modu- 
lated the torsional response, and caused continuous stall 
and/or unstall of the blade over a significant portion of 
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a .degrees 


Figure 4 AIRFOIL SECTION CHARACTERISTICS FOR NACA 23012 
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Figure 5 DISPLACEMENT TIME HISTORIES AT CLASSICAL FLUTTER ONSET 
n*-5.3, 0 O =* 11 deg, m ” 0 







a revolution, due to the large plunging rate generated by 
the flapping motion. An example of this occurrence is 
shown in Figure 7. Thus, while stall flutter involves only 
the rotational degree of freedom, the results obtained 
indicate that the minimum speed for its occurrence is deter- 
mined by coupling with a translational degree of freedom. 

Results for the hovering case are summarized in 
Figure 8, which compares computed and measured flutter 
speed and frequency, plotted against collective pitch angle. 

No upper limit in collective pitch angle for the occurrence 
of stall flutter was calculated, since that limit would 
depend strongly on initial conditions, and so would be 
arbitrary. Quantitative differences between the computed 
and measured stability boundaries of Figure 8 can be attrib- 
uted in large part to the use of a two-dimensional aerodynamic 
model, which cannot precisely reproduce the aerodynamic 
coupling between the rotational and translational degrees 
of freedom. 

From the basic similarity of the computed and measured 
stability boundaries and the character of the computed in- 
stabilities (Figures 5 and 6) it can be concluded that the 
aerodynamic and dynamic models formulated are capable of 
reproducing both classical and stall flutter as experienced 
by a rotor blade, and so can be employed to investigate the 
forward-flight problem. 


Stability in Forward Flight 

The nominal configuration was analyzed next for an 
advance ratio of .1. Computations were carried out in the 
same sequence as for hovering. First, the rotational speed 
at which classical flutter occurs was determined. Then, 
stall-related instabilities were investigated. 

A linear bending-torsion instability of the Floquet 
type (Ref. 14) was encountered at £2* =5.2. Blade 
motions as a function of azimuth angle at flutter onset 
are shown in Figure 9. The torsional and bending displace- 
ments are seen to display the aperiodic character typical 
of this type of instability. The flapping motion is the 
steady-state response to the cyclic pitch variation. 

An instability analogous to stall flutter in hover was 
found to occur for £2* as low as about 4.4, with collec- 
tive pitch angle greater than 12 deg. Blade motions for 
£2* =4.8 are shown in Figure 10. The torsional displace- 

ment time history, while not strictly periodic, is nonetheless 
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COLLECTIVE PITCH ANGLE, 9 0 . deg 


Figure 8 FLUTTER SPEED AND FREQUENCY VARIATION WITH COLLECTIVE 
PITCH ANGLE FOR A HOVERING ROTOR 
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brought about by successive stall and unstall. The azimuth 
positions at which those events occur are marked by (S) and 
(U), respectively, on the \p -scale. 

The blade motions for the type of instability shown in 
Figure 10 are not of the same character as those of particu- 
lar concern in the limiting of helicopter performance, in 
that the excessive torsional displacements shown in Figure 
10 persist over a complete revolution of the blade. The 
control load time history, taken from flight test (Ref. 6), 
shown in Figure 11 illustrates the type of stall-related 
blade motions usually encountered at a thrust level or 
forward speed near the upper limit of an aircraft. Large 
oscillations in the control loads, presumably deriving from 
blade torsional oscillations, are seen from Figure 11 to 
persist only between about \p = 270 deg and = 400 deg, 

rather than throughout a complete revolution of the blade. 

A torsional displacement time history closely resembling 
the variation of control loads in Figure 11 was obtained 
for less than 4.4, for collective pitch angles between 

12 and 13 deg. Results for two typical cases are shown in 
Figures 12 and 13. The occurrences of stall and unstall 
are indicated on the abscissas. The large oscillations in 
torsion are clearly related to stall, but their persistence 
is not the result of successive stalling and unstalling, 
as would be the case for true stall flutter. The blade 
appears to be responding to the sudden changes in aerodynamic 
moment at stall onset and unstall, as can be seen by compar- 
ing the variation of moment coefficient shown in Figures 12 
and 13 with that of torsional displacement, and noting the 
azimuth positions at which stall and unstall occur. There 
is some cyclic stall-unstall within the stall zone evident 
in the results, particularly at the higher rotor speed 
( ft* = 4.15, Figure 13). However, the major contributors 
to the oscillations appear to be the initial and final 
pulses associated with stall and unstall upon entering and 
leaving that zone. There are, in general, two cycles of 
torsional oscillation of excessive amplitude after the blade 
unstalls the last time on a given revolution. The response 
can be regarded as transient, on a localized time scale, or 
forced, when viewed on a scale of several rotor revolutions. 
The severity of the response is apparently due in part to 
the suddenness of load changes at stall and unstall, and 
partly to the relative lack of aerodynamic damping in pitch, 
particularly when the blade is not stalled. 

If the collective pitch angle is increased, the blade 
does undergo stall flutter, as seen from the time history 
plotted in Figure 14. "These results are for the same rotor 
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LINK LOAD, » 



AZIMUTH, ^ , degrees 

Figure 11 VARIATION OF PITCH LINK LOAD IN FLIGHT 
TEST OF CH47 AT 123 KNOTS 
(from Ref. 6) 


29 




Figure 12 DISPLACEMENT AND 
TORSI 
12* = 3.8! 




AERODYNAMIC MOMENT 
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speed as those of Figure 12, but with 0 O increased from 
12 deg to 14.3 deg. Successive stall and unstall persists 
over the whole revolution of the blade for this case . 

It could be argued that the blade torsional oscillations 
of Figures 12 and 13 are still a manifestation of stall 
flutter, even though successive stall and unstall Is not 
taking place, since the aerodynamic moment can undergo 
unstable variations when the blade remains stalled through- 
out a cycle (Ref. 4). It may, in fact, be the case that 
the large deflections do result partly from that effect, so 
choosing to term them as simply a response may be somewhat 
misleading. On the other hand, the solutions are distinctly 
different from what is definitely stall flutter obtained 
both in hover (Figure 6) and In forward flight (Figures 10 
and 14) so that label would seem to be even less appropriate. 
Further, the persistence of the oscillations after exit from 
the stall zone is clearly symptomatic of a response, so, fGr 
lack of a more precise term, solutions of the type shown in 
Figures 12 and 13 are identified in what follows as exces- 
sive response . 


Linear Stability Boundaries 

The value of SI* at the onset of linear instability 
was determined for the three configurations considered, 
for advance ratios of 0, .1, .2, and .3. The effects of 

advance ratio and torsion-bending frequency ratio on linear 
stability are shown in Figure 15, where ft* is plotted 
against n for two different frequency ratios. Increasing 
advance ratio is seen to cause some decrease in flutter 
rotational speed, with most of the decrease occurring 
between advance ratios of .1 and .2. The substantial de- 
crease in frequency ratio, from 3.69 to 2.5, caused only 
about a 4 percent reduction in flutter speed over the range 
of advance ratios considered. The insensitivity to frequency 
ratio can be attributed to the large chordwise mass imbalance, 
which produces the same effect in classical binary flutter 
of a wing (Ref. 15 ) . 

The effect of chordwise mass imbalance on linear sta- 
bility is shown in Figure 16, where ft* at flutter onset 
is plotted against M for values of x m of .216 and .108 
semichords. As one would expect, the reduction in x m , and 
hence in the coupling between bending and torsion, causes 
a substantial increase in the flutter rotational speed. 



♦ 



Figure 15 EFFECT OF ADVANCE RATIO AND 
TORSION-BONDING FREQUENCY RATIO 
ON LINEAR STABILITY -Xm/b - 0.216 
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Stall Flutter and Response Boundaries 

The effect of forward speed on stall-related Instabili- 
ties for the three configurations was investigated by 
systematically varying the collective pitch angle and 
advance ratio, with S2* equal to 3.89. In order to relate 
the results to rotor performance, a mean lift coefficient 
(Tr is defined, according to 


1 

Cy Z 5 * 

h P SI* b 

where 1 is the time-averaged lift per unit span at the 
aerodynamic reference radius. This coefficient is, to a 
good approximation, directly proportional to the thrust 
coefficient (see Ref. 16). The two-dimensional aerodynamic 
model does not provide a good measure of ZTr when the rotor 
is partially stalled, so Cl was computed assuming it varies 
linearly with the collective pitch angle, using the formula 


C L = a( n )(e c + .0217) 


The slope a and zero-lift collective pitch angle of -.0217 
rad were obtained from calculations of Cr for the nominal 
configuration with stall precluded. The variation of a 
with p is shown in Figure 17. 

The results obtained for the nomi-nal configuration are 
summarized in Figure 18 as a plot of Cr vs p . As thrust 
is increased at a given M , the rotor is seen to first 
encounter a region of excessive response, of the type dis- 
cussed previously, and then, for p of .2 or less, a 
region where stall flutter occurs. Increasing advance 
ratio has the effect of suppressing the tendency_for stall 
flutter. At p = .2, stall flutter occurs at Cr = .85, 
but a further increase in Cr results in excessive response 
again. At P = .3 a limit-cycle type of oscillation 
could not be triggered at all. As a result, stall flutter 
is confined to a region somewhat as indicated by the shaded 
area in Figure 18. 

The suppression of stall flutter at high advance ratio 
is apparently caused by an effect similar to the one en- 
countered at low rotor speed in hover, whereby the flapping 
motion prevented a limit cycle from occurring. This can be 
seen from the blade motions obtained for p = .3 and 
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Cl = . 78, plotted In Figure 19. On the first revolution, 
as the blade enters the stall zone on the retreating side, 
it appears that a limit cycle is being set up, with repeated 
stall and unstall occurring. However, at about \Jj = 420 
deg, the flapping motion has built up in response to the 
large cyclic pitch changes, producing a negative plunging 
rate sufficient to keep the blade unstalled over the remain- 
der of its passage on the advancing side. Then, when the 
blade again enters the stall zone, the large positive flap- 
induced plunging rate precludes unstall until exit from 
the stall zone at about \p = 670 deg. As a result, the 
blade subsequently undergoes excessive torsional response, 
rather than stall flutter. 

The effect of torsion-bending frequency ratio on stall- 
related instabilities can be seen from Figure 20, where Cl 
is plotted against H for w e 0 / w ^ 0 = 2.5. No instance 
of excessive torsional response occurred with this config- 
uration for an advance ratio of .2 or less. Instead, 
limit -cycle type oscillations were set up, with almost no 
evidence of suppression by the flapping motion, even at 
relatively high values of Ul with M = .2. At a* = . 3, 
however, only excessive response was obtained, similar to 
the results for w 6(/ u )?So = 3.69. 

The marked deterioration in stability at the lower 
frequency ratio is apparently associated with the lessened 
linear stability of the system. The configuration with 
Xm/b = .108, which is more stable. In the linear sense, than 
the nominal one, exhibited a trend opposite to the one re- 
sulting from a decrease in frequency ratio. The results 
for the smaller mass center offset, shown in Figure 21, are 
similar to those of the nominal configuration. Figure 18, 
but the region in which stall flutter occurs is somewhat 
reduced, there being no occurrence of stall flutter at an 
advance ratio of .2. Also, the amplitude of the torsional 
oscillations in the region of excessive response is con- 
siderably reduced, as evidenced by comparing the blade 
motions plotted in Figure 22, which are for M = .1, 

Cl = .95 and Xm/b = .108, with those of the nominal config- 
uration plotted in Figure 12. 
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Figure 19 DISPLACEMENT TIME HISTORIES AT HIGH ADVANCE RATIO - 
J2*«3.89, C L = 0.78, m = 0.3 
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Figure 20 STALL STABILITY BOUNDARIES FOR ft* = 3.89, co ft /om = 2.5 

AND Xm/b = 0.216 0 ° 
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ADVANCE RATIO, ft 


Figure 21 STALL STABILITY BOUNDARIES FOR ft* = 3.89, uq / W 0 O = 3.69 

AND Xm/b- 0.108 
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Figure 22 DISPLACEMENT TIME HISTORIES FOR EXCESSIVE TORSIONAL RESPONSE. 
n* = 3.89, C L = 0.95, n = 0.1, AND Xm/b = 0.108 
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CONCLUSIONS 


An analysis has been performed of the aeroelastic 
stability of a helicopter rotor blade in hovering and 
forward flight. An analytical model of an airfoil under- 
going unsteady stall and an elastomechanical representation 
including flapping, flapwise bending and torsional degrees 
of freedom were employed in the study. The following con- 
clusions can be drawn from the results obtained. 

1. Analysis of aeroelastic stability for a 
hovering rotor demonstrated that the aero- 
dynamic and dynamic representations developed 
are capable of reproducing classical and 
stall flutter. 

2. While stall flutter is an instability 
involving a single rotational degree of 
freedom, the minimum rotational speed 

for its occurrence, in hover, is determined 
from coupling with a translational degree 
of freedom. 

3. In forward flight, the rotor can undergo 

a linear instability analogous to classical 
flutter and a stall-induced flutter which, 
while not manifested by a strictly periodic 
limit cycle, has the same basic mechanism 
for its occurrence as stall flutter of a 
hovering rotor. 

4. The large stall-related torsional oscillations 
which limit forward speed and thrust are 
primarily the response to the rapid changes 

in aerodynamic moment which accompany stall 
and unstall, rather than the result of an 
aeroelastic instability. 

5. Linear stability is relatively insensitive 
to advance ratio for advance ratios as 
large as .3. 

6. While excessive response due to stall occurs 
at high advance ratio, stall flutter is pre- 
cluded by the large flap-induced plunging 
rates . 
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7 . 


The severity of stall-related instabilities 
and response depends to some extent on 
linear stability. Increasing linear sta- 
bility lessens the susceptibility to stall 
flutter ; and reduces the magnitude of the 
torsional response to stall and unstall. 
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APPENDIX A 


PROGRAM LISTING 


APPENDIX A 
PROGRAM LISTING 

A listing of the FORTRAN coding of the computer program 
follows . The program was written in FORTRAN IV for use on 
an IBM 360/75 computer. 



|oo i I do no ! ! ioqo j o O o oo 


PROGRAM TO ANALYZE UNSTEADY AIRFOIL STALL 


COMMON /flL 1 / NT I ME * NDIMC , ISTD 

COMMON /CLCMBL / CLVB , CMVB , CMPAVB 


COMMON / I NPTVB/ 

FTVBI64I , 

FPVBI64) 

, FPPRVBI 641 , D10RVBI64I, 

4 

XMVB C 64) , 

OELVB, 

XMUVB, 

FOVB, 

XMU4VB, 


B 

ATOVB v 

ATCVB, 

ATSVB, 

ROVB 

, RVBI64I, 

C 

MVBC6*) , 

NVB 





COMMON /INPUTS/ 

NSBL, 

NZ, 

NO FF, 

NGAM, 

NSIG, 

A 

NC 01 f 

NCORO, 

LOWER, 

MSTOP, 

MAX T , 

MDTR, 

B 

NOTBL « 

INOV, 

ELSIG, 

OXI, 

REB, 

ROBB, 

C 

FRZ, 

ARR, 

AMPLU* 

FREQU , 

ALPH1, 

ALPH2 , 

0 

HEAVE, 

AROT, 

FREQF, 

PHIH, 

NY, 

RYi, 

E 

DRY, 

y i ioo ) , 

TEST* 

UPRIM, 

XU< 30), 

YU( 30 ) 

E 

XL C 301 , 

YU33) , 

EIU, 

ER2, 

ER3, 

B03R, 

G 

RRDBR 






H t 

C MPA , C MPAS , 

BARG, EMI, 

HVOR, NVOR, SSPA, 

SVOR, TtlRF, 

X IVOR 

I t 

PLOTOP, PSILOW, 

PSIUP 






J , NOUT 

COMMON/ III/ 2 <3 I 


MAIN 2 
MAIN 3 
MAIN 4 


SETUPS17 
SETUPS1 8 
SE TUPS19 
SETUPS20 
SETUPS21 
SETUPS22 
SETUPS2 3 
SETUPS24 
SE TUPS2 5 
SETUPS26 
SETUPS2 7 
SETUPS2 8 
SETUPS29 
SETUP S3 0 


SETUP S3 1 


OIMENSION US AV ( 3 00 , 1001 ,SCALS(300) MAIN 5 

DIMENSION US AVI l ,1 I .SCALSI300) ~ '' MAIN 5 

DIMENSION CAMBRI24I ,THICKI24> MAIN 6 

DIMENSION XGAHI3O»,XSlGIldOJ,XSIGAadO» t xSieft(ldd),XCI3MTTxno0D7MArN — 7 

1SBLI300) (XBSIGIIOO) MAIN 8 

Dt ME NS I ON AC API 30*3) , BC API 100 ,3T, ASZ f30 > /AS <30, 30I,BS(30. 30I,ASHZMAIN 9 

IllOO) ,ASHf 30,30 l,BSHI30,30 I , ARI 30) , ARHI 1001, UEI 300,3 » MAIN 10 

DIME NSI ON ALAMI30) » VZIPI 30 I »FPRESI 100 ), GAMAWTIOOOF. XTvTC 1000 » MAIN 11 

OIMENSION BLAMI30) , FLAMI 10 J , X FL AMI 101 MAIN 12 

DIMENSION $CALE(300,2I,U(l,l,l t.UCI 100, 3), VI 100,2) MAIN 13 

1 , PI200,7> 


DOUBLE PRECISION CMATI60.60 ),RMAT< 1301 ~ " MAIN 15 


DATA IN, MOUT, NF/ 5,6, 24/ 

DATA PI ,TIME,UINF,RENEL.UST0P/3.l4l59,d.7T7;47r5E4,"?vB7 MAIN 18 

DATA FLAM /I .75 ,1 .75, 1 .724, 1.527, 1.354, 1., .663, .452, .25MAIM 19 

I4777T7 ~ ~~ " HA IN 2D 

DATA XFLAM /-100. , -1 1, 26.-7.0 1,-3. 48 ,-l. 766, 0., I .888 , 4. MA IN 21 

103,6.77,7.19/ : MAIN — 22 

OATA OEGRES /1.74 53292 51994 3300-2/ SUPPL 38 


EOUIVALTNCfc ICHATIl I ,US AV (11 1 ,( ASHI1I , SC ALSU) 1 


HAIM 16 


r> o n n 


I STD* l 

RAD * 180. /PI 
l L* 8888 


NDIMC^ 60 


CALL '5FTTF5 

IFUSTO . EQ. II GO TO *0 

00 100 J I t 300 

SCALStJI *0j» 

DC 1 00 I * 1 ♦ i 00 

100 LSAV (Jfl) *0 
40 CONTINUE 

C 

CALL RFADIN C IL,G 60 1 

NOTE - CFFSETS“A“RE PUT IN AS LISTED IN THEORY OF WING SECTION St 
AS A FRACTION OF TOTAL CHORD, XI BEING MEASURED FROM THE 

leading edge* Bfc Sure Nf Is an even number.. 

1 1 ME *0 • , 

NTI ME*0 

NWAKE * 999 

I SEP*0 

isFFr^D 

I WASH *2 

uiw «n 

L*0 

1 NO V*I NO V*1 
WRI TE I MOUT #6 1 

jT t TC # » AUPHl 

IFIINDV ♦ MOTR .LE. 21 PITCH « PITCH - ALPH2 

f F n Wv~TF57~TI 

X AMPLU * 1.33333* XMUAV8 • tl.-R0VB**31 / < 1. - R0VB**4) 

I F ( 'TNDV. FQ. 21 FREOU* WBR/ftRMTT 
IF C INDV .GE. 21 GO TO 343 

WKt TE I MO'UT ,251’ WdftVSvOft f HV0R,S ARG,X ivaR,lMT f T 0 RF, S SPA 

R Y*R Yl 

HVOR^VOR**2 
BARG *BARG/6.2832 

^43 CALL SECT i XU , YU/JO^TVLTN'OIT^tWF^RWBTTflOSBftM 0 BB, TH ICK,CAMBR) 


MAIN 


MAIN 
I .EMAIN 
MAIN 
MAIN 
MAIN 


I MAIN 
MAIN 


MAIN 

MAIN 

MAIN 

MAIN 

MAIN 


OG 7875 N*l ,NF MAIN 

FAMBR ( Ni *C A MHRTN'I *C W8S ' MAIN 

7E75 THICK! NI»TH!CKCN)*T MDBB MAIN 

“ WRITE (MW,41 ; “ MAIN 

WRITE! MOUT ,7) AMPLU ,FR6QU,ALPHi , ALPH2, HEAVE, AROT , FREQF, ROBB #REB MAIN 

girnTTTOUTTBT : r ‘ MAIN 

WRI TF I MOUT, 91 ( N ,CAMBR( N) ,THI CKC N) t N* 1, NF I MAIN 

__ m Xa:NSRl ^ M2! - I ' MAIN 

CALL SCAL! SBL,NS8L,FRZ, ARR,RDBBI MAIN 

CALL C ORDX I NSBL",NZ , RT3BB # TBL,X,XCT MAIN 

00 2420 M*'l ,MX £ MAIN 

TFTXC TMT-l . 1242072 4 19,2419 MAIN 

2419 MEND=M-i MAIN 

CT'TO'2421 MAIN 

2420 CONTINUE MAIN 

2421 WT*HFND MAIN 


4 


65 

66 

67 

68 
69 


59 

72 


75 
64 

76 

77 

78 

79 

80 
81 
82 
33 
84 
35 
86 

87 

88 

89 

90 

91 

92 

93 

94 
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MXM1=MX-1 
PE (MX+1 ,l»=l. 

EPSLE=2.*!X<NZl-X(NZ-l) j 

FPSTF = X! MX) -X 1MX-1 ) 

ALTC=8.36F4/SQRT!REB) 

I F ( I STO.EQ, l) GO Tn 50 
00 2422 M* 1 t MX 
SCALE I M, 1 ) =0. 

SCALE ( M, 2) =0. 

00 2422N=1 ,_NY 

u(m,n,u=6. 

2422 UIM,N f 2)=0. __ 

5 C CONTINUE 

NSIGA = NSIG 

NSI 3R=NSt 3 

NS1G1=NS!G*1 

M0TR = M0TR«-1 
N0TBL=NQTB L*1 
XMA X* 1 .-ELS I G 

CCNA* .375«PI/OXI 

ANGS=PI /FLOAT(NSIG) 

_ CALL SETSXINSIGI >l»lt2« »XSIGt ANGS ) 
XSEP«l.l 

00 2430 N=1 tNSIGI 
XSIGB(N)=XS!G (N) 

2430 X SIGA( N)*XSIG(N» 

^DO 2431 N*i f NSIG 

00 2431 NU-1,3 

2431 BCAP!N f NU)*0. 

PINT*2. /FLOAT! NCORO) 

NCP1*NC0RD*1 

THXI *1.5/0 XI 

NGPl *NG AM* 1 

NWM1 *NWAKE— 1 

C CUNT=0. 

00 B456 N*1 *NWAKE 
~ GAMAN! N) *0« 

XIWIN)=1. ♦COUNT 
8456 COUNT*COUNT«-OXI 

ANGLE=PI /FL0AT1NGAM) 

C0PNT*0. 

DO 1002 M*1 ,NGP1 
PHI M=C OUNT* ANGLE 
XGAM!M)-COSIPHIM» 

D0UNT*2. 

00 1001 N“2 tNGAM 
AS(M,N) =COS!DOUNT*PHIM> 

10C1 D0UNT«D0UNT«-1. 

1002 COUNT*COUNT*l. 

CALL WASH! XGAM,NGAM,TIME,ALPHl, ALPH2, 
rTi?B»»NF,VZtP,l ,1» 

00 8458 M*1 «NGP1 

rwr r TWTT T ^n 

TEMP-2. *VZ!P!M) 

RMATIHI-TEHP 
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S 



! o o n n o n 


C MA T I M , 2 J-XGAMfM) 

00 8457 N*3 ,NGPl 
8457 CMATIM, N)=AS(M,N-i) 

B 458 CONTINUE 

CALL ALSOU NGP1 ,CMAT,RMAT) 
DO 8459 N= 1 , NG PJ^ 

ACAP(N f 1 )=RMAT IN) 
ACAP(N t 3)=RMAT JN1 
8459 AC A P ( N , 2 1 * AC AP ( N , 1 1 
00 2 704 M = 1 , M X 


2774 

2715 


SIGN=1 • 

IFIM-NZ1 2774,2775,2775 
$IGN*-ST5N~ 

CALL QECALf l SEP , NGAM ,NS l G, NF, XS IG,ACAP, BC AP, THI CK , R C88 , GAMA 1) 


1 NF , XC ( MfTuFTMal ,S I GNJ 
2784 UF I M f 2 ) -UE(M,l) 

00 r004"M=2,'NiGAM 

1004 BLAM(M)g(l« 125* XG AM ( Ml » , 1 875*1 1 . +XGAM I M ) I *1 l.-3.*XGAM<M ) )*ALOG( 

i ♦ xga m ( m) ) / n . - xg a h r wi nr/ oxt 

BLAMING PI) *— 1 #125/0X1 

C ALL C LC M ( NCOI , I SEP , NGA M, X SIG , NS IG, XS IGA, NS I GA, XS IGB,NS IGB, ACAP 
1AP , THI CK ,R0BB,GAMAW tUINF ,UOOT ,DX I , AROT, CMPA) 

TF 


1CALL 


tt mr 

SUPPL 


iFQ'.Tr 


INDEXING IN TIME IS CARRIED OUT AT THIS POINT. 


9SS9 CCNTINUE 

era: 


IF ( 


'Acocpor 

IACU * LT. 


TACOT 

35000 


) GO TO 99 


NOTE - FOR READ-IN CF FCIL 
SLPH2 


MOTIONS, MAKE ALPH1 * ALPHA, 
AUP R A'- D OT i " AND HE Ay E 


"H-D5TT 


l SB 


IFIMCTR .EU. 21 

XREAO (l N,2 , END* 8 9 89) 

NITS=1 “ 

TIME = TIME+DXI 


ALPHl,ALPH2, HEAVE 


NTT ME - * N TTMF+1 
NWAKE*NTl ME + 2 

IF I NWAK'E-998 )' 

201 NWAKE*998 


202 ,2017201 


202 
8 800 


IRMAXT-NTI MET 
SAVE U-Ul NF 
C*“ 


29297B80072200- 


P(L 

PST3F0 


or 
,11 * 


BCBR 

WfiDT” 


/ RRDBR 

PCL,TT 7 


TIME 

I6CT.T 


RAO 


Ul NF*1* + AMPLU*SI N(FREQJ*T I ME) 

IFiTNOV 


■7F0T~2T 

XC ALL SUPP1IUINF1 

PTTTTH * ALPHI 

IF ( ! NOV ♦ MOTR .LE. 2) PITCH * PITCH 

■UDOT=FRE OU*AHPLO*CO’S i FREQU*TTWT 

STEPX*. 5*DXI*( UINF+S AVEUJ 

T50 XOOT “7*2 ,TI WATTE 


- ALPH2*C0S( FREQF* T IME I 


MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIM 
MAIN 
MAIN 
MAIN 
MA IN 
MAIN 
, UI MAIN 
MAIN 
MAIN 
MAIN 
( 1*MA IN 
MAIN 
MAIN 
,BCMAIN 
MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 


MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 


MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 


150 

151 

152 

153 

154 

155 

156 

156 

157 

158 

159 
16 0 
16f 
162 

163 

164 

165 

166 

167 

168 

504 

505 


169 

170 

171 

172 


175 

176 

177 

178 

174 

182 

183 

184 

185 

186 

187 

188 
189 


190 


475 

191 

192 

193 
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JC=NWAKF-J+2 MAIN 

GAMAWI JC)=GAM4WI JC-l) MAIN 

ICC? XIW< JC)=XIW(JC-l)+STFPX main 

IFUSEP) 2 009 ,2 OOP, 2007 MAIN 

2CC7 DG 2008 N=1,NSIG ~ MAIN 

BCAP(N,3)=BCA P(N ,2) _ MAIN 

20C8 BC4P (N,2 ) =BC AP (N»l ) MAIN 

00 4433 N= I ,NSIG1 MAIN 

XSIGRINI -XSIGA (N ) MAIN 

A A 33 XSIGA1N )=XS IGIN) _ MAIN 

GO to" 2010 " " MAIN 

2C09 DEADL-O. _ __ MAIN 

E LOnT=UI NF - MAIN 

2010 00 1014 M= 1 | MX _ MAIN 

UE <M,3)*UE IM,2> MAIN 

101 A UE(M,2)-UF (M,l) MAIN 

DEA01-DEADL " MAIN 

E LD 1*EL00T MAIN 

ALAM II ) - I 1 • 125+. 75* ALOG I STEPX4.5 I I /DX I MAIN 

00 1 005 M- 2 , NG PI MAIN 

1005 ALAM(M) -BLAM(M) *••75* II •♦I 1.— X GAM I M ) ) / STEPX )* ALOGI I l.+STEPX- XGAMIMIMAIN 

1) /Il.-XGAMIM)) )/OXI MAIN 

DC 2006 M=l,NGPl MAIN 

ACAP<M,3)*ACAPIM,2) MAIN 

20 C 6 "AC API M,2)=ACAPIM*1) MAIN 

AFACT=8. *( ACAPI1,2)*.5*ACAP12,2))-2.*I AC API 1, 3 > ♦ . 5* ACAP ( 2, 3 II MAIN 

ALPHS-VZIPIl) MAIN 

CALL WASH! XGAM , NGAM ,T I ME , ALPH1, ALPH2, HEAVE,AROT , FREQF.PHIH, UINF ,CAMAI N 

1MBR,NF,VZIP* MOTR » I NOV ) ~ MAIN 

DO 1006 M=l,NGPl MAIN 

ASZ(MI-l.*2.*ALAM(MI MAIN 

AS(M,l)-XGAMtM)*ALAMIM| _ MAIN 

SUM-0. ~ ’ MAIN 

DO 4343 J-2,NWM1 MAI N 

4 343 SUM-SUM* (G A M AH I J I «■ ( GA MA M ( J*l I -G AM AW I J ) ) *! X GAM I M )-X I MI J 1 >/! X I WI j* U MA IN 
l-XIWI Jl) l*ALOGI IXIWI J+U-XGAMIMI )/ IXIWI J J-XGAMt M ))) MAIN 

ELX-l.-XGAMIM) MAIN 

IFIM-1) 1006,2130,1006 _ MAIN 

2130 ELX-l. MAIN 

1CC6 AR(M)-2.*VZIP(MI*ALAM<M)*AFACT/3.*ISUM-GAMAW<2)*< l.-X&AMIM) )*ALOG(MAIN 

1 1 l. + STEPX- XGAMI M) )/ ELX) /STEPX 1/ PI ' MAIN 


: 1HE FOLLOWING CALCULATIONS, THROUGH STATEMENT 4444, ARE PERFORMED 

ONLY IF THE AIRFOIL IS STALLED. THE AIRFOIL IS DESIGNATED TO BE 
; STALLED IF INTEGER ISEP IS NONZERO. 

: IF(ISEP) 3247,4444,3247 

3247 GO TP 13344,3345) ,1 WASH 
3 344 X$EP-XSEP*OXl 

IF { XSEP-XMAX) 3248,3347,3347 
3 347 I WASH *2 
I SEP-0 


DO 3015 K* 1 ,3 
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n H 


3015 RCAP ( N, K) -0, 

GO TO 4444 

! F a NOT) 334 3 ,3348 , 3248 
IF(NITS-I) 3248,3349,3248 


3 345 
3348 
3 349 


Iff I NDV.EQ.2) GO TO 6349 
tF<V7FP(l)-ALPHS) 6349,6348,6348 

6348 NITS*? 

GO TP 3248 

6349 CALL UNPOP ( NGAM , AR » ALAM, AFACT , RM AT , CM AT , X GAM, AS , ACAP ,MX ,N2, 
1 ,BCAP,THICK,RPBB,UINF f XC,UE) 

GO TO 2785 

3248 XATT=XSFP+DEAD1**5*( ELD1 +ELDOT ) *DX I 
DFADL=XAf T— XSEP 
0 IFF =1 . -XATT 

XTFST * XSEP ♦ 3* * EPSLE 
CALL SE T SX( NSI G l ,XSEP,XATT ,XS IG, ANGSI 
DO 4434 N-l ,NS I G 

4434 XBSIG(NI*«5*( XSI G (N ) +XS IG ( N+l ) ) 

DO 3 086 M*i,NGPl 
DO 3086 N«l , NSI G 
3 G86 BS(M,NI*d; 

DO 308? M- 1 , NG P i 


3 C 09 
3C9i 


IFTxGA'NT'HT-xSFPr" 

IFIXATT- XG AM (M ) ) 

no 


3187,3087,3091 


3093,3092,3092 


3092 I*i,N$lGl 
IF (XGAM(M)-XSIG(l) ) 

3093 MARK«I 

GO TO 3094 

“3X97 C CNTT FFUE 

3CS4 W f DE S-XSIG I MARK) -XS l G (MARK-1 I 

B S TM , M"A ft K-I r*l X STG I M ARX1 H<GYMTmT T7 W 1 0 E S 
B SIM,MARK) *( XGAMf Ml -XSI G ( MARK— 1 ) )/W IDES 
B STM ,1T « SQRTT i XG m CM) -X S E P17 f x ATT -X GA M (M))) 

3088 IFIDlFF-i.E-6) 3087,3098,3098 

3 C 58 ' B S { M ;n « R SIM ,T JTDTF F**T- 1 . 5 J ’♦'S'flRT'l DE AOL )"* t 27*15 IFF* { WTTf'I . 
1) ) /( XATT- XGAM f M) ))-!•»♦( 4. *XG AMI Ml -U-3,«X ATTI) 


MAIN 
MAIN 
MAIN 
MAIN 
main 
MAIN 
MAIN 
MAIN 
I F , X S I GMA IN 
MAIN 
MAIM 
MAIN 
MAIN 
MAIN 

MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MA IN 
MAIN 
-XGAM(MMAIN 
MAIN 


25 i 

252 

253 
? 54 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 


GO TO 3087 MAIN 286 

3187 BS(M f i)*DIFF**f-l«5) *SQRT I DEADL ) *1 3«4 X ATT-4. *XGAMf Ml ) MAIN 287 

'3067" CONTINUE MAIN 288 

C MAIN 289 

SET-UP' OF THE SECOND SET “OF EQUATIONS START'S HERE • MAIN 290 


D0"435(TK*I7NSTG 

IF I XBSIG (K)-l* 1 4348,4349,4349 

^ C CSK = XB"SIG'm 

SINK^SQRTCl CQSK*COSK) 

THE T K* A RC TIC OS'KJ 

TANT=SIN(.5*THETKI /COS I .5*THETK I 
A SHTCK} « TA NT *C O N A *"!T •♦CTJSXI’^’IT'V = 
ICGNA* 1 1 • +COSK) *SI NK**2 1 /UINF 
A 5HTK, 1 ) *;5t C ASNZTKT-T ANTT^STNK" 
C CUNT-l • 


MAIN 291 
MAIN 292 
MAIN 293 
MAIN 294 
MAIN 295 
2<?6 
MAIN 297 

-3 . *C0 SK J/UINF+THX T* { P I - THE T K ♦ SI NK ♦ M A I N 298 

MAIN 299 

__ MAIN 30Q 

MAIN 301 


D’0 ’43'55 “N=2 , NGAM' MAIN 302 

COUNT-COUNT+i, MAIN 303 

4355 ASH IK, N)=S1 NIC 0UN T * T H ET K1 4. 7 5 »( S I Nil COUNT 41.mHETKI7fCPUNT+ 1,1- SI MAIN 304 
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1N( (CrUNT-1.) ♦THETK) / (COUNT-1. ) )/ ( DXI*UINF) 
GO TC 4350 
4349 A SHZ ( K) =0. 

DO 4359 N*1 f NG AM 


4 359 
4350 


A SH ( K f N) s 
CONTINUE 


0 . 


MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MA IN 


IF (0 IFF-l • E-6) 5005 ,5006 ,5006 
50C5 PREC -0. 

GO TC 5007 

50C6 CALL A TTPRI PRF C # XSI GfNSJGt ASZ , AS , AR L CMAT , RMAT, N GAM, NF, A CAP , THICK ,RMA IN 
IDBB ,GAMA W t Ul NF ,UOOT ,DXI ,BCAP) ' MAIN 

50C 7 CALL MI XER CYPRES ,PREC,U I NF , JDOT, THICK , NF, XBS IG, NSIG, INDT,DEL1?TH€TMAIN 
i 1 1 RFB , USEP » X4 tCPi ) MAIN 

C PCT*C.P1 MAIN 

DO 4800 K*i,NStG MAIN 

CORD = XB$ IG f K) MAIN 

MAIN 
MAIN 
1) tXSMAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 


8 SH ( K* i ) *-1 • +THXI *8 1 NT I XSEPt X ATT , CORO )/U INF 
DO 4808 N-2 tNSI G 

4 8C8 B SH C K , NT *F B < XS l G < N- L ) * X S IGCNI , XS I G( N*l ! , CORD ) +T HX I*GB ( X $ IG( N- 
1IG(N) ,XSIG (N+n , CORD) /U INP 


CALL 

CALL 


ESIGI (2,NSIGA, XSIGA , 9CAP, CORD, VAL 1 I 
ESIGI (3,NSIG8,XSIGBtBCAP,CORD,VAL2) 


ARH(K) =FPRES(K)M2. *VALl-. 5*V AL2 )/ ( DX I*U INF ) 

IFCCORD-UI 5008,4800,4800 
5CC8 CALL EGAMl ( 2 , NG AM , AC AP , BC A P ( 1 ,21 , X S I G A( 1) ,XS IGA ( NS IGA+1 ) , GA MAWI 2 ) ,MA IN 

_ 1C0RD ♦ VAL l ) _ MAIN 

CALL EGAMl ( 3 , N G A M , A C A P , B C A P ( 1 ,3) *XSIGB( l),XS IGB(NSIGB+l ) ,GA MA W( 3 ) , MA IN 
1C0RD , VAL2 ) MAIN 

ARH( KHARH(K)^(2.*VALt-.5^VAL2l/ (OXI^UINF I + .062 5*AFACT* P I*( 1 . KOROMA INI 
l) *(1.-3.*C0RD*THXI*(1.-C0RD*C0R0) )/< 0XI*UINF1 MAIN 

4 SCO CONTINUE MAIN 

4444 CONTINUE MAIN 

C " V MAIN 

C CALCULATIONS FROM THIS POINT ON COMBINE THE MAIN 

MAIN 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 


C 

C 


CASES OF STALLED AND UNSTALLED AIRFOILS 

DO 6 5 00 M*TTNG P L 
RMATIM) * AR ( M) 

CMAfIM,!) 


00 64 85 


ASZfMf 
N=l tNGAM 


l>\t¥-tnXYTW,H+l i * AS IT , N) 

IF ( I SE PI 6486,6500,6486 
6416 DO 6499 N*l,NSIG 
NGG*N+NGPl 

6499 C MAT ( M,NGS 1 * B§ ( M ,Nl 

6500 CONTINUE 

“ IF ( I $E PJ 6"S'02 t6501 t 6 S 62 
6 5 Q1 NTOT -NGPl 

‘ GTrTO" 6T51 — — ” 

65C2 00 6750 K*1,NSIG 


305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

31 7 

318 

319 

32 0 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

33 5 
336 
33 7 
33 8 
33 9 

340 

341 

342 

343 

344 


MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 


345 

346 

347 
34 8 
34 9 

350 

351 

352 

353 

354 

355 

356 


RMAT IKK) *ARHIK) 

C MA T ( KK JTT *A5H Z ( K ) 
00 6748 N*1 , NGAM 


MAIN 

MAIN 

MAIN 

MAIN 

mnr 

MAIN 


w 

MAIN 

MAIN 


T5T 

358 

359 
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on 6750 N= 1 , NSI G 


MAIN 


NGG=N+NGP1 IN 

6150 CMAT(KK f NGG)=ftSH(K»N) MA I(M 

NTOT=NSIGf NGP1 ma !,N 

6151 CALL A L SOL ( NTOT , CMAT ,RM AT ) ' MAIN 

DO 6800 N= l , NGPl MAIN 

68CC ACAPtN, l) = RMAT t N I ' MAIN 

IF(ISEPI_ 6805,6320,6805 MAIN 

66C5 DO 6810 N= l ,NSIG MAIN 


NGO =N+NS Pi MAIN 

6610 BCAP (N, l 1 *RMAT ( NGG) " MAIN 

6e20 CONTI NUE _ MAIN 

GAMA AMI (ACAPfOXItPl ) MAIN 

IF< PS 1360 .GE. PSI LOW .AND. PSI360 .IE. PS 1UP ) GO TO 1736 
DO 1785 M=1 ,MX ~ MAIN 

SIGN*! • MAIN 

IF ( M-NZ ) 1780,1785,1785 ' " MAIN 

1780 SIGN*- SI GN MAIN 

1 785 CALL QEC AL (I SEP , NG A M , N$ I G, NF , XS I G, AC AP, BCAP, THICK, RDBR, GAMA *l( 1) ,UIMAIN 
INF, XC(M) ,UEJ^M,U , SIGN) MAIN 

2165 00 8886 1=1 ,2 " MAIN 

US2 = UEIl,l) MAIN 

DO 8886 M* 1 , M X M 1 MAIN 

US1 =UE ( M ,1 ) MAIN 

uftm ,n^nrsr^2 ♦uetm*i f 1 n / 5 . main 

6886 US2*USi MAIN 

GO ""TO (8351 ,8353) , I WASH " MAIN 

8351 DO 8352 M=1,MX MAIN 

83T? SCAlS(*T=0. ~~ MAIN 

GO TP 1786 MAIN 

8353 CALL YSTT l R Y1 , Y ( 2 ) , NY,Y) ' MAIN 

RY*RY1 MAIN 

■""'60 8354 M= 1 , MX '' MAIN 

8354 SCALSIM>=0. MAIN 

IT<lNDv7EQ.2) GO TO" 8370 ' ' MAIN 

IFUSEP.EQ.O.ANO.VZIPin.LT.ALPHS) GO TO 1786 MAIN 

8370 CALL ST AG < MX, NY , MSTUP,*ST ,DXI , RY“, 0RY7X, Y,OT,UC, V,OSAV, SCAL S , ISEP ) MAIN 
LAMQ*i MAIN 

XSEPS*XSEP ' MAIN 

DXX=DXI MAIN 

IF U SEP. EQ. 1 . AND. IS EPT • EQ.7J. AND. N ITS. EQ. 1 ] D>CX=L.FJO~ MAIN 

8 367 CALL BLC (X,Y,MST, MEND, NY, RY, DRY , OXX, RE8,UPRIM, Ft AM, XFL AM , TE Sf ,U, SC MA IN 

lALE,UEVOC,V,XSEP,USEP,OISF",THETA,LOWlR,LAMQ,TrsiP,KC, USA V , SC AL S, N I T MA I N 

IS, NT! ME, NOTBL , XTEST , N 2, NOUT ) 

IT ( XSE P- XMAX ) t 736 , 7735 ,7735 MAIN 

7135 IF 1 1 SEP) 1786,1786,7736 MAIN 

7736 DELl*DrSP ~ KAm 

THE T1 = TMET A MAIN 

I NOT *1-1 AMO ~ MAIN 

IF ( I NDT.EQ. 1 • AND.NOTBL. EQ.2) GO TO 1786 MAIN 

WRI TETMOUT , 2 3T XSTG TT ) , CPCT, XSFT" ’ “ * MAIN 

IFCINDT) 8462,8462,8463 MAIN 

8462TTIISEP) 8562", 8562T8563 ™IN 

8563 IFINITS-I) 8562,8562,8662 MAIN 

5662 IF (I SEPT) 7742,7 742,8562 MAIN 


36 0 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

37 5 

376 

377 

378 

379 

380 

381 

382 

383 

384 

386 

387 

388 

389 

390 

391 

392 

393 

395 

396 
39 7 

398 

399 

400 

401 

402 

403 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 
41 5 
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85 62 CAU. RU9B(DEU ,THFT1 ,RFB,XSEP ,USFP 

IF ,altc*rfnel,ustopi 

USE P = USF P* # 602046*USFP**3 
PIMFf « (JOSE P-U5) * JUS E PHJ 5 > 

WR I TFTMOUT , 22 ) PD I F F , DC P 
I F ( DCP-PDt FF ) 82 63 * 8366 , 8365 
8263 ISEPT*0 

GO TO 8463 

8 366 IF (I SEP I 8 368,8368,8369 
8369 IFUSFPT) 8467,8467 ,836 8 
8 66 7 I WASH = 1 


,XC5,DCP,DEL5,X,XC,MX ,NZ, X5,U5,UMAIN 416 

MAIM 417 
MAIN 418 
MAIN 419 
MAIN 42 0 
MAIM 421 
MAIM 422 
MAIN 423 
MAIM 424 
MAIM 425 
MAIN 426 


NI TS=2 
GO TO 3344 

8 368 GO TO < 3 a786 ) t_NOTRL 

8168 CALL REA tT (UC, V, X i Y ,MX*NY, RY, DRY ,UE, X5, DEL 5, MS T , R EB I 
L AMQ=0 
GC TO 8367 

8463 IF(ISFP) 7741,7741,7742 
7741 I SF P=l 


MAIM 42 7 
MAIM 428 
MAIN 429 
MAIM 430 
MAIM 431 
MAIN 432 
MAIN 433 
MAIM 434 


NIT S-NI T Sfr 1 


IF(INDT) 7743,7743,7643 
7643 I SEPT*1 


0 XSE P=1 • — XSEP 
XSEP=.6*XSEP+,4 


CALL CPC (I SEP,NGAM,NF,XS IG,NS IG, XS IGA, NS IG A, XS I GB , NS IGB 
1 THI CKtRDBR ,GAMAW ,UI NF, UDOT , 1« , XSEP, OX I, CPI) 

GO TO 3248 

7742 CALL FLDER ( BCAP , XS I G ,NS I G ,U I NF, ELDCT , S IGSUM, YMX 1 

I FTT SF P . E Q ♦ 1. AND . I S F PT . E 0 • 0 . A NO . N I TS . EQ. U GO TO 9210 
IF < XSEP*. 5) 7841 ,7842,7842 

7841 EPS*EPSLE 
GO TO 7843 


7842 E PS*EPSTE 

7843 PXSEP=AB SI XSEP-XSEPS) 


IFIDXSEP-EPS) 7834,7834,9210 
7834 IF ( XSF P- XMAX1 1786,1786 ,78 35 
7635 I SE P*0 
I SE PT a O 


00 7836 K=1 ,3 
00 7836 N»1 ,NSIG 
7 636 BCAPI N,lO*0 • 

GO TO 1786 


9210 NI T S*NI T S* 1 

IF(NITS*EQ.2 .AND. I NOT • EQ* 0 ) XSE PS* XSEP 
IF INI T^-4) 9211,9211,1786 


9211 IF I XSEP- XSE PS ) 9305,9305,9306 

43C5 x$FP*.6*x5FP$*.4*X$eP 

GO TO 9307 

ttCZ “X3ET>^76*XSEP*.A*XSEPS 

9307 IFIXSEP-XMAX) 9212,9212,7835 

9212 CALL CPC ( I SEP,NG AM, NF,XS IGtNS IG, XSISA, N$ lGA,X S f W, NS IGB 
1 THICK, ROBB, GAMAW,UINF,JOOT,1. ,X$EP,DX I, CP1I 

TFI NOTBL «EQ« 1 .AND. XSEP .GT. 0.) XSEP— .48 

GO TO 3248 

7 743" IF(NITS^T) 7737,7737,3248 


MAIN 

MAIM 

MAIM 

MAIN 

MAIN 

, ACA®, BCAP, MAIM 
MAIM 
MAIN 
MAIM 
MAIN 
MAIN 
MAIN 
MAIN 
MAIN 

MAIM 

MAIN 

MAIN 

MAIN 

MAIN 

MAIM 

MAIM 

MAIM 

MAIN 

MAIM 

MAIN 

MAIM 

MAIN 


43 5 
436 
43 7 
43 8 

43 9 

440 

441 

44 2 

443 

444 

445 

446 
44 7 
44 8 

44 9 

45 0 

451 

452 
45 3 

454 

455 

456 

457 
45 8 

459 

460 

461 

462 

463 

464 

465 

466 

467 


MAIN 

MAIN 

MAIM 

MAIN 

, AC AP iUC AP ,"M A fN 
MAIM 


MAIN 

MAIM 


46 8 
469 
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7737 N I T S = N I TS+1 
ELDOr-ELDl 
GO TO 3248 
1786 WRITE ( MOUT, 20) 


NTT ME 


WR t TE ( MO UT , 26 ) X1VOR 
PI TC = PITCH * 180. / PI 
2 C9 WR l TE (MOUT, 10T TIME, UINF,XSEPfX ATT, PIT C 
ALOFG= ALPH1/DEGRES 

WRITE <6, 9001") Z f ALDEG, ALPHl , ALPH2, 
IF! P SI 360 .GE. PSILOW .AND. PS 1360 
IF( NOUT * EO. 0) 

1WRI TE IMOUTfli) 

IF< NOUT • EQ. 6 ) 

1 WRI TE I MOUT « 12) 

2N=1 , NGPi ) 

IF (I SEP) 7432,7433,7432 


MA IN 
MA l\| 

main 

MAIN 

MAIN 


47 0 
*71 

472 

473 
477 


HEAVE 

.LE. PS IUP) 


GO TO 101 


MAIN 473 
SUPP 1.34 9 
SUPP L 350 


MAIN 479 


EQ. 0) 

( N, XG AM< N ) «VZIP|N), AR(N), ACAPIN, l ), X lW( N ) ,GA 1AW< N ) ,MA IN 

MAIN 
MAIN 


480 
48 L 
482 


7432 m NOUT .EQ. 0) 

1WRI TF ( MOUT ,13) 

I FT NOUT . EQ. OT “ 

1 WRITE (MOUT, 17) ( N,X BSIG I N ) , FPRES ( N ) , A RH( N ) , BCAP ( N, l),N=l,NSIG) 

e i n n-r 


MAIN *,3 3 


7433 


WRITE (MOUT ,14) 
WRITE (MOUT, 13) 


WR I TE ( MOUT , 15 ) 
XPC=-1. 


ELOOT 

XSIG(1),CP0T,X4,CP0T,XATT, PREC 


MA If) 
MAIN 
MAIN 
MAI)( 
MAIN 
MAIN 


DO 7102 N*1 »NCPl " MAIN 

CALL OECAL(I SEP, NGAM,NS I G»NF» XS I G, ACAP, BCAP, THI CK , R CBB, GAMA W( 1) ,UIMAIN 
INF ,XPC ,QFL,-i.) ' WAIN 

CALL QECAL(ISEP,NGAM,NSIG,NF,XS I G, AC AP, BCAP, THICK , RCBB, GAMAWI 1) ,UIMAIN 
INF , XPC , QEU * l . ) “ " ' " MAIN 



7102 

TOT 


cpl=Cpu 

DLIFT=CPL-CPU 
WRITE (MOUT, 16) 
XPC=XPC*PI NT 


X PC , OEL , C PL , fliEu , CW.DLTFT 


MAIN 

MAIN 

MAIN 

MAIN 


CONTINUE 

C MPA S=C MPA MAIN 

CALL CLCM( NCOl , I SEP , NGAM.KSI G,NS IG7XS IGA.NSTGA, XTIGB.NTTGB, ACAP , BCMA IN 
1 AP, THI CK, ROBB, GANAW,UlNF,UOOT,DX I,AROT,CMPA) MAIN 

pm SI = ntrr " “ 


thtit 

P<L,3> 


wrc 

Z (3) 


PTC74T 

P( L , 5) 


TTTT 

Z (2 ) 


P(L,6) = tLV* 

P I L , 7) * C MPA 

IF I L .LT. 200 ) GTT 

CALL PLOTSB ( PLOT OP 


TO WIT 


L ) 


5 8 


CONTINUE 


iFUsrn .EQ. rr go to 

DO 7950 M*1 ,MX 
SCALETM,"?! * SC ALE (M, 11 


484 

48 5 
466 

487 

488 
439 

49 0 

491 

492 

493 

494 

495 

496 

497 
49 8 

499 

500 

501 

502 

503 

504 

505 


MAIN 

PAIN 


506 

507 
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n o n n 


SC ALE ! M 1 1 ) == SC ALS ( M) 

DC 7950 N=1 f NY 
U<M f N,2) *U(H,N, n 

7550 U<M f N,l) -USAV(M »N) 

GO TO 9999 

8589 CONTINUE 

99 CONTINUE 

CALL PLOT SB! PLOTOP , P , L ) 
CALL AC UCPU ( I ACU ) 

IF < I ACU «LT. 350 00 I GO TO 60 
GO TO *0 

60 CONTINUE^ 

I F ( PLOTOP • EQ* 0.) CALL EXIT 
CAU PLTN_0_ 

CALI EXIT 

RETURN 


MAIN 503 
MAIN 509 
MAIN 510 
MAIN 51 L 
MAIN 512 
MAIN 513 


F ORMAT ( 1315) 

F OR MAT ( 3F1 0« 4 ) 


FORMAT ( 2F 1 0« 4 J 
FORMAT (1H1 // ) 


MAIN 

main 

MA IN 
MAIN 
MAIN 
MAIN 


FORMAT(6F10*4) 

FORMATdH l , 50X >34HAN ALYS I S OF UNSTEADY AIRFOIL STALL///) 

FORMAT ( 8X*6HUBAR = El 3 .5 / 7X ,7HUFRFQ =F13*5//3X» 1 1HALPHA ONE =E13«.5/MAIN 
13X,11HAIPHA TWO =El 3.5/8X »6HH8A R = E 13 . 5/ 1 IX , 3HA =E 1 3 .5/ 8X» 6HFRE Q *MA IN 
1E13.5//8X,6HR0/B =E 1 3. 5//9X ,5 HREB *E13.5/// ) MAIN 

8 F0 RMAT(29X,lHNt25X t 4HC! M y 26X t 4HT < N) / ) MAIN 

9 EORMAT(!30»2E30*5) MAIN 

10 FORMAT! 5 X f 3HT =E 13. 5/ 5X ,3HU =E13 .5/4X* 4HXS = E13 ♦ 5/4X , 4HX0 *E13*5/4MAIN 

1 X # 4H PA -E13.5////I MAIN 

1 1 FORMAT ( ///4X »1HN tl 1 X tlHX , 14X»5HVZ (X) tUXtSHRNlX )♦ 1 2 X » 4HA(N) t 21X , 3H MA IN _ 

1X1 W,14X,5HGAMMA/1 MAIN 

12 F0RMATU5 f 4EI7.5,8X,2El7.5) MAIN 

13 FORMAT (IHl ,8X f lHN,20X»iHX f 21X *5HFP|X) ,22X, 5HRH! N) f 21X» 4HBC N I /I MAIN 

14 FORMAT! //54Xf9H L-OCT * E 13. 5///5 IX, 2 7HPRESSUR ES IN SEPARATED FLOWMAIN 

X//55 XtlH Xf 19 Xf 2HCP/ ) MAIN 

l 5 FORMATdHl f 1 l X , 1 HX, 16X* 3HQELt 15X, 3HCPL* 15X» 3HQEU* 15X» 3HCPU* i3X»9HCMAIN 

' ” ~ “ ' main 

MAIN 
MAIN 
MAIN 
MAIN 
MAIN 


l PL - CPU/) 

16 F0RMAT(6E18.5) 


17 FORMAT! I10»4E25«5) 

1 8 FORMAT! 3 ! 40X ,2E20. 5 / )) 


19 FORMAT! I5 f 5F10*4l 

20 FORMATdHl »50X t 12HTlME STEP N013//J 


~22 FORMAT! / //40Xf26HI NCREASE IN CP REQUIRED ISE 13 . 5//4GX* 2 6HINCREASE MAIN 
1IN CP POSSIBLE ISE13.5) _ MAIN 

T3 FORMAT! / //45X,23HP0f ENT I AL^rCW XS *El2* 4/60X, 8HCP( X SI *ET2*4/MAIN 

i /45X 1 23HBOUNDARY LAYER XS *E12.4) MAIN 

24 FORMATII5f4Fl0.4/5Fl6.4l “ MAIN 

25 FORMAT! 12 Xf4HNV «I2»3X t 3HS «E12«4t3X t 3HH *E12*4» 3X» 3HG *E12 »4t 3X>4MAIN 


23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
_48 

49 

50 

51 

52 

53 

54 


lHXl »Cn;4//12X f ^RB! »EI?^,3Xf4HwT"-ei2r47Tx7"4HPA -E12.4//71 MaTn 53“ 
26 FORMAT (4 X f 4HX1 *E13.5I MAIN 56 

96Ci FdftMATPO 1 •EQUIVALENT ROTOR BLADE RESPONSE* SUPPL380 
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9CCIA 

// 

T 5 i 

•FLAP DISP »• t G14.5 



SUPPL381 

9CC1P 

♦ 

T47 f 

• BENDING MS P =• t G 14 ,5 



SUPPL3H2 

9KIC 

♦ 

T1 9» 

• TORSIONAL MSP G14.5 



SUPPL3R3 

91C1D 

/ 

T38, 

•SECTION PITCH ANGLE »», 

F9.3, • 

DEGREES HRS 

SUPPL3R4 

90C1E 




F 9 • 4 » f 

RADIANS • 

SUPP L3R*> 

9CCIF 

/ 

T21, 

•SECTION PITCH RATE «• , 

G14.5 


SUPPL3R6 

90 Ci G 

f 

T71, 

•SECTION PLUNGING RATE =' 

t GL4. 

5 // ) 

SUPPL387 


END MAIN 515 
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SUBROUTINE SUPPL 
1 MPLIC l T RE AL*8 tA-H,0-Z> 

REAL*8 FRlSt FR2S * FR3S , 4NSX, QMS 


REAL*4 CLVB, CMVB » CMPAVB 

1 , DUMMY, PLOTOP 

REAL FTVB, FPVB, FPPRVB, DIDRVB, X 

A FO VB , XMUAVB, AT 0 V B , ATCVB, ATS 
C WOXI , PSI , UINF 

REAL ELS1G, D XI , RE Bj _ RDBB , FRZ , 

A ALPHl , ALPH2, HEAVE, AROT, FREQF, 

B X, TEST, UPPIM, XU, YU, XL, YL, 

C RRDBR 

REAL SUM ( 8 ) , YCLD(8), YNEW ( 8 ) , DELC3,3 
A Z , ZPRI3) , SMALLG(3I , Y(3,3), YPRC 
COMMON /BL1 / NT I ME, NDI MC 

COMMON /CLCMSL/ CLVB, CMVB, CMPAVB 

COMMON/ III/ zm 


DIDRVB, XMVB, DEL VB, 
ATCVB, ATSVB, ROVft, 


AMPLU, 


PHIH, 

ERl, 



SUPPL 

1 


SUPPL 

2 


SUPPL 

3 


SUPPL 

4 

XMUVB, 

SUPPL 

5 

RVB, MVS. 

SUPPL 

6 


SUPPL 

7 

FREQU, 

SUPPL 

8 

3YI, ORY, 

SUPPL 

9 

ER3, BOBR, 

SUPPL 

10 


SUPPL 

II 

CL ( 3 ) » G( 3) * 

SUPPL 

12 

3) 

SUPPL 

13 


CMPAVB 


FTVBI64 ) , 

DELVB, 

ATCVB, 

NVB 

NSBL, 

NCORD, 

I NDV , 

ARR, 

AROT, 
XUOO) , 

YL (30 S t 


FPVB164) , 
XMUVB, 
ATSVB, 


FPPRVB1 64), DIDRVB! 64) , 
FOVB, XMUAVB, 

ROVB, R VB (64), 


COMMON / I NPTVB/ FT 
A XMVB ( 64 ) , DE 

B ATOVB , AT 

C MVB ( 64 1 , NV 

COMMON /INPUTS/ NS 
A NCOI, NC 

B NOTBL, IN 

C FRZ, AR 

0 HEAVE, AR 

E DRY, X( 

"P XL (30) ; YL 

G RRDBR 

H , DUMMY CIO) , PLOTOP 
DIMENSION DE LTA (3 ,3 ) 


1 , TOT ( 2 ) 



Z2IX)MFZ/FRlS*FRiS^CF4!X)-F2+( I .-C6*X*X ) *( R2-BZ/FR IS )) *X*X 
SI I X) * ( 2 . *hS*CF 4 CXT7G8**r*( FR I^FR2ST^ *X ) *GA 
S2(X)*(FR1 S-FR2S ) *G A*X^X 


SUPPL 15 
SUPPL 16 
SUPPL 1 7 
SUPPL 18 


NZ, 

NOFF, 

NGAM, 

NSIG, 

SUPPL 

19 

LOWER, 

MSTOP, 

MAX T , 

M3TR, 

SUPPL 

20 

ELSIG, 

DXI, 

REB, 

RDBB, 

SUPPL 

21 

AMPLU, 

FREQU, 

ALPHl, 

ALPH2, 

SUPPL 

22 

FREQF, 

PHIH, 

NY, 

RYI, 

SUPPL 

23 

TEST, 

UPRIM, 

XUI 30), 

YU( 30 1 , 

SUPPL 

24 

TRT, 

ER2, 

ER3, 

BD3R, 

SUPPL 

25 

... 




SUPPL 

26 





SUPPL 

27 

> 3 1 , GAMMA! 3, 3 ) , OMS < 3 ) , OM EGA! 3) 

,CMK<3> 

SUPPL 

28 

20), ANTC 20), AAXf 10) 

, AN SX( 20 ) , 

SORT* 3) 

SUPPL 

29 


SUPPL 

SUPPL 

SUPPL 

SUPPL 

SUPPL 


iXI *S2 ( XI I 


. , WIT H I 


“ANU HAD TOST 


C AS REFERENCES. NONROTATING NATURAL FREQUENCIES ARE 


SUPPL 
SUPPL 
SUPPL 
SURF L 
SUPPL 


C AND S2LB ARE FRACTIONS OF SEMICHORD. X B AR, S IL , AND S2L ARE 


SUPPL 

SUPPL 

SUPPL 


00 63 K » i, 8 


SUPPL 47 


63 YNEWfK) * 0. 


SUPPL 4 9 
5UPPL 50 


60 


66 


67 

69 


DO 66 K = 1 » 8 
YOLD(K) = YNEWtK) 
CALL YV3lYNFW,TT 
IF (l .LE. 1) GO TO 
6T 


69 


DO 67 K 
SUM < K) = 

C CNTI NUF 
F Ml 1 = SUM ( 1 ) 
EM22 = SUM (2) 
EM33 = SUM 1 3) 


~U~Q 
( YNEWIK) 


♦ YDLD< K) ) * (RVRm - RVB( I-ni / 2. ♦ SUM IK) 


EMI 3 = SUM (4) 
E M2 3 = SUM ( 5 ) 
Mil = 

H22 


SUM(6> 
= SUM! 7) 


H33 = - EM33 
H13 * -E Ml 3 
H23 = SUM(8) 
BDBRR=BDBR/RRDBR 
8 D S =BD B R R* *2 
Til =H1 1 *BD S 
T22=H22*BD S 
T33=H33*BD S 

TI3=HT3*BfT5 

T23 =H23 *BD S 

F R l S=BD S *E Ri **2 -T1 1 / EMI 1 
FR2S=ER2**2*BDS-T22/EM22 


FR3S*FR3 **2*BOS-T33/ EM33 
FR1=DSQRT(FRIS ) 

Fft2-D$QRT(FR2'S> 

FR3=DS0RTCFR3S) 

RSTM^E Mil AEM72 

ZETA*(l#+RATM)*( RAT M*FRIS**2*FR2S**2 )/(RATM*FRlS + FR2S)**2 

RM=7FTA-T. ” 

SUMS*FRl S+FR2S 

HlGMS=(SUM5^S0ftTrST7W**2-V^2ETA^lS*ER2sy)/f 2V*ZETaT 

SMALS*FR1$*FR2 S/HIGHS 

~TJEN=F"R2 S-F RI'S 

Ai=-<HISHS-FRI$)/DEN 

A2*-l'i-4i 

A l* A2*DEN/H I GHS 

~~5L A M T* E M 1 r *B TJBR * *Z7TW3 3 

SLAMZ*-A 1* SLAMl 

5LAM2--SLAMZ/ A2 

SUM3 *SUM S*F R3 S 

FRIS*(FK2 S4FR3S)+FR2S*FR3S — - • 


■soar 

ADDZ 

“RR7TR 

B4*$ 


‘FR1 S*FR2$*FR3S 

: l.-ltMl3**2/EMH4E M2 3**2/EMZ2 T7"fcM 33 


IEM13*¥2/EMITT/EM33 

B4*R4 /BB AR 

B2=ADt)2»- (2.*PR2S*fcMI3*Tl3/EHl l+2.*FR!S*EH23*T23/f 
123**2/EM22J/EM33 

B2=B2 / BB AR 

BZ=ADDZ-(FR2S*Tl3**2/EMll«-FRlSn23**2/EM2 2l/EM33 
BZ*BT/BBAR 



SUPP L 

51 


SUPP l 

52 


SUPP L 

53 


SUPP L 

54 


SUPP L 

55 

) 

SUPP L 

56 


SUPP L 

57 


SUP® L 

50 


SUPP 1 

60 


SUPP L 

61 


SUPP L 

62 


SUPP L 

63 


SUPP L 

64 


SUPP L 

65 


SUPP L 

66 


SUPP L 

67 


SUPP L 

68 


SUPP L 

69 


SUPP L 

70 


SUPP L 

71 


SUP PL 

72 


SUPP L 

73 


SUPP L 

74 


SUPP L 

75 


SUPP L 

76 


SUPP L 

77 


SUPP L 

78 


SUPP L 

79 


SUPP L 

80 


SUPP L 

81 


SUPP L 

82 


SUPP L 

83 


SUPP L 

84 


SUPP L 

85 


SUPP L 

86 


SUPP L 

87 


SUPP L 

88 


SUPP L 

89 


SUPP L 

90 


SUPPL 

91 


SUPP L 

92 


SUPPL 

93 


SUPPL 

94 


SUPPL 

95 


SUPPL 

96 


SUP® L 

97 


SUPPL 

98 

2 S* SUPP L 

99 


SUPPllOO 


SUPP LI 01 

1- 

TSUPPL102 


SUP® L 1 03 


SUPPllCK 


SUPP LI 05 


SUPP LI 06 
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C 6*<EM11*A1**2 + FM2 2*A2**2)/EM33 
F4*$UP3 

C4=(FR2S*FM11*AL**2+FR1S*FM22*A2**2)/ FM33 

GA=?.*EMU*AJL/FM33 

GB=2.*FM22*A2/EM33 

F2 “A002 

HA*EMLL/FM33 

HB«EH22/EM33 

FZ*ADDZ 

Ri=-HA-HB*tGA/GBt**2 
R2=HA«(FR2S7friS-U 1 
ZLAM=F 

TVlL AM=B4 *C6— C4 
F ZHA T=HB*{ ZL AM/GB^**2 

F2H AT*B2-F2+FRi S *ZL AM+2 • * 7 LAM *T WL AM*H B/GB**2 
F4HAT=-C6»B2»FR1 S*TWLAM+HB»I T WL AM/GB) *»2 
G2HAT=B2-F2+ (FZ-BZ ) /FR1S+FR1S ♦ZL AM 
Qi»HAT*-C6» tB2-B Z/FRlS)»FRlS*T WL AM 
SIGZ-2# ♦HB*ZLAM*G A/ GB**2 
SfG2 ^GA»(FRi$-FR2 S»2,»HB»TWLAM/GB*»2) 
GAM2*GA^(FR1S~FR2S) 

UZ=-R2*F ZHAT 

Ul-R 1 *G2 HA T— R2 ♦FZHAT 

U2-R 1*G4HA T— R2 *F4HA T 

U3=-R2*SIGZ 

UA=R1*GAM2 -R2*SIG2 


U5=SIGZ*G2HAT-GAM2*FZHAT 

U6=SIGZ*G*HAT*S[G2*G2HAT-GAM2*F2HAT 

U7 JC $IG2*G4H AT-G AM2*F4HAT 

AAX(U*UZ**2 

AAXI 2)*2 • ♦ UZ*U1 *U3*U5 

AAX<3)*Ui**2*2.*UZ*U2+U3*U6HJ4»U5 


SUPPL107 
SUPPL108 
SUPP L 109 
SUPPLUO 
SUPP L 1 1 1 
SUPPL11 2 
SUPP LI i 3 
SUPP L 11 4 
SUPP L 11 5 
SUPPL116 
SUPP LI I / 
SUPP LI 1 8 
SIJPPL11 9 
SUPPL120 
SUP^ L 12 i 
SUPP L 12 2 
SUPPL123 
SUPPL124 
SUPPL125 
SUPP H 2 6 
SUPP L 12 7 
SUPP L 12 8 
SUPP L 12 9 
SUPP L 1 30 
SUPPU31 
SUPPL132 
SUPPL133 
SUPP L 13 A 
SUPPL135 
SUPPL136 
SUPPL137 
SUPPL138 



AAXI4) =2.*Ul*U2«-U3*U7«-J4*U& 


SUPP LI 3 9 


AAXl 5) =U2**2>U4*U7 


SUPPLUO 


CALL POLLY (4tBBSfRELfANSXtAAX) 


SUPP L 141 


XBAR=l • E 25 


SUPP L 142 


DO 86 1*1,4 


SUPPL143 


1 P=2*l 


SUPPL144 


IM*IP-i 


SUPPL145 


IFIDABS(ANSXUM) I.GT. 1.0-10) GO 

TO 86 

SUPPL146 


IF (ANSX(IP).LE.O.) GO TO 86 


SUPP L 14 7 


XBAR T=D SORT (ANSXUPII 


SUPPL148 


IF(XBART*LT.XBAR) X8AR*XBART 


SUPPL149 

86 

CONTINUE 


SUPPL150 


IFIXBAR.LT. .5625) GO TO 88 
WRITE (6 f 87) 


SUPPL151 
SUPPLl 52 

" &r 

FORMATIIHl ,10X,'N(3 SOLUTtCN FOR 

XBAR' ) 

SUPP L i 53 

6 8 

STOP 


SUPP L 154 

CONTINUE 


SUPP H5 5 

1 5 

AL0W*IRI*Z2I XBAR)-R2*Zl (X BAR ) )/ 1 R1 *S2 ( X BAR l-R2*S 11 XBAR ) 1 

SUPPL156 


ALOW*ALOW/XBAR 


SUPPL157 


BL0W*(CF4(XBAR)-GA*AL0W«XBAR)/(XBAR*€B) 

x i »-~a lov-b cnw 

E TA« t BLOW*A l-AL0W*A2 1 / 1 Al- A2 I 


SUPPL 15 8 
SUPRU59 
SUPPLUO 
SURPL161 


62 



S1L*( XI -RM*HIGHS*S2L) ♦HIGHS/ IFR1S*FR2S) 

WRI TF ( 6 * 4 ) FR1 ,ER2,EF3,RM 

WR ! TF ( 6 » 72 l ) FR1 ,FR2 ,FR3 , ALOW , BLOW 

WRI TF ( 6 , 5 ) EMI i,E‘12 2,EM33, FM13, EM23 

WR l TF (6,61 Hi i , H 22, H3 3, HI 3, H 2 3 

C 13* ALOW/BDBR 

C23*BLQW/BDBR 

XRA8*XP AR/BOBR 

SI LB = SI L /BDB R 

S2L9*S2L/BDBR 

WRITE (6,41) BOBR , RROBR 

WRITE (6,71 XB AR, X9AB, Si L,S1LB,S2L,S2LB,SMALS, HIGHS 
A A ( 1 ) *B7 
AA ( 2 ) *B2 
A A ( 3 ) *B4 
A A ( 4 ) =1 • 

XALL POLL Y l NPOL , 8B$ ,R£L,ANB,4AI 


SUPP 1162 
SIJPP L 1 6 3 
SUPPLL64 
SUPP L i 65 
SUPP L 166 
SUPPl 167 
SUPP L I 6 3 
SUPP L 16 9 
SUPP L 1 70 
SUPP L 1 71 
SUPP L 1 72 
SUPPl 173 
SUPP L 1 74 
SUPP LI 75 
SUPP LI 76 
SUPP L 1 77 
SUPP LI 7 8 


S S X* SLA MZ* XB AB 

D I V = 1 . - $ L A M Z * X B A R ♦ * 2 

BETA (3 ,1 ) = (SLAM1#C13+SSX*FR1S)/0IV 

BE tA (3,2 ) = I SL A M 2 * C 2 3 ♦ S S X^FR2S ) / 01 V 

BETA (3,3 )*(FR3$*SSXMC13*C23))/DIV 

A XB"=AT*'XBA’B 

BETA ( 1 , l > *FRl S— AX3*8ETA (3,1) 

B E TSTr,TT^ A XB *8 ET A (3 ,T> 

BETA ( 1 *3 )=C13— AXB*8ETA(3 ,3) 

A A X B * A2 ♦ XB A B 

BETA (2 , 1 ) *— A AXB*BET A (3 , 1 1 
BETA (2, 2 »*F B2S-A AXB*fiEt A (3,2 i 
BETA(2,3)*C23-AAXB*BETA<3,3) 

A B ( 4 I * 1 * 

AB(3)*BETA(1 , l ) 4-BET A (2 ,2 I «*BET A( 3, 3 ) 

ABTZT=BE TA ( 1 , 1) * I BETA ( 2 , 2 ) ♦BETATT7 3 f ) ♦ BET A ( 2 , 2 ) *BETA( 3, 3)~B£TA( 
1I*BETA(2,3)-BETAU,2)*BETA(2, 1)-BETAI l, 3 1 *BETA( 3,1) 

JTB (1) ^BTTA ( 1 , 1) * ( BE T A (2 , 2 ) +B ET A ( 3 , 3 1- 8 ET'AT 37 JWE T A ( 2, 3 ) )-B ETA( 

l)*(BETA(l,2) ♦BETA (3 ,3)-BETA(3,2)*BETA(l,3))«*BETA( 3, i)*( BETA ( 1,2 
1 FTAT2",TT-BETAf I ,3) *BETAC 27211 
CALL POLLY ( NPOL ,BB$ , REL, ANT, AB) 

WRITE (6, 44") 

00 45 1 *1,4 

“ IM=(r-Tl“*2 

45 WRI TE ( 6 , 46 ) I M, A A( I ) , AB ( I ) 

WRITR6;4T1 
DO 48 1*1,3 
TTT*T*1 
I TM=ITT-1 

WR I TE ( ATI ANB'iiTT ) , ANB'MTM) , ANT ( I M ), ANK I TFT 

DO 301 I *1 ,3 
TT^2r*l 

301 OMS( l ) *— ANT (II) 

~ Minrr*3 

DO 70 1*1,2 

FF l wen . GT . OHS ( H A XI ) ) MAX f * I 

70 CONTINUE 

GO TO (TT , 7 2 ,'7 3 T» MA XI 


SUPP L 17 9 
SUPP L 190 
SUPP LI 81 
SUPP LI 02 
SUPP Li 8 3 
SUPP L 1 84 
SUPP L 1 8 5 
SUPP L 186 
SUPP L 1 87 
SUPP L 1 8 8 
SUPPL 189 
SUPPL190 
SUPP L 191 
SUPP L 192 
SUPPL193 
3, 2SUPP L 194 
SUPPL195 
2,1 SUPPL 196 
) *8 SUPP L 197 
SUPPL 198 
SUPPL199 
SUPPL200 
SUPPL201 
SUPP L2 02 
SUPPL203 
SUPPL204 
SUPPL205 
SUPPL206 
SUPPL207 
SUPPL 20 8 
SUPPL209 
SUPPL210 
SUPPL21 l 
SUP 3 L212 
SUPPL213 
SUPPL214 
SUPPL215 
SUPPL216 


I 


63 


71 


I i =2 
12=3 

GO TO 74 

I L =1 

12=3 

GO TO 74 

I 1 = 1 

12=2 

IFtOMSU 1) .GT*3MSU2)) GO TO 75 
MINI = 11 

Mini *12 

GO TO 76 

Ml NI =1 2 

M 1 01 * l 1 _ _ 

SORTm=nMS(MINI ) 

$0RT(2)=0M$(M10n 

SORT (3)=0MS!MAXI) 

00 77 1=1*3 

oMsm=soRTU> 

OMEG A ( I ) =0SQRT tC MS ! I ) ) 

00 302 1=1 *3 
ALPHA! 1 * t ) = 1 • 

0EN9 *RE TA(2fi)*8ETA(3*2 ) -BET A ( 3 * i ) *C BET Al 2, 2 )-OMS( 1 > ) 
ALPHA (l ,2) =<BETA(l ,2 >*BETA(3, II -BETA! 3,2)*( BETA I l* 1 l-OM 
ALPHA (t *3) *1 (BETA (2 *2)-0M$ (1) )*! B ETA! l* 1)-0MS( 1) )-BETA( 
L*l) ) /DENB 

CHK ( 1)=BETA(1*3) * ALPHA < l * l )+BET A! 2*3 ) ♦ALPHA! 1* 2 )♦( BETA! 
l)*ALPHAtl,3) 

DFNB=BE TA ( 3 72T*T5eT A!l*i)-CMS(2 FT— BET AC 3 * I) * BET At 1*2 ) 
ALPHA(2*l)*(BETA(3,l)*( BETA ( 2 ,2 )-OMS( 2I)-BETA(2, 1)*BETA 
ALPHA (2 * 3) * ( BETA (2 * 1 ) *BET Ml * 21 - (BET A I 1, 1 )— OMSt 2 I )* ( BET 
12)) ) /DENB 

CHK! 2)=8ETA (1*3) ♦ALPHA! 2*1 I *BET A! 2*3 ) ♦ALPHA! 2, 2 )*(BETA( 
1 ) *ALPHA (2*3) 

OENB=BETA!2*3)*( BET A (t * 1 ) - CMS (3 ) I -BET A! 1* 3) *BET A! 2* 1 ) 
ALPHA (3*1) = (BETA(2fi)*( BETA! 3*3 )-0M$( 3U-8ET A< 3,1 )*BETA 
ALPHA ! 3 , 2T^i BETA! 3 * lT*BET A ! 1 * 3) -! BETA ! 1, 1 )-nMS( 3 ) ) * ( BET 
131)) /DENB 

CHK(3) = BETA(1,2)*AL PH A! 3 , 1 ) + 1 BETA (2,2 )-OMS( 3))* ALPHA! 3, 
1)*ALPHA(3,3J 

WRltn6*4B6r ' ~ ' 

HR I TE (6,489) ( I , OMEGA! I ) * BET A ( I , l ) , BET A< I * 2) , BE T A! 1 , 3) * 
l ALPHA ( I ,2) *ALPHA(1 *3) *CHK (!) *1*1*3 ) 

SORTCII-U 

" Sort T? ) =0* “ 

SORT<3)=0. 

DO 4*2 J = i ,3 

GO TO (381 ,382 *383) , J 

SORT ( 1T=0. 

$0RT<2)*1. 


SUPPL217 
SUPPL2I 8 
SUPP 1219 
SUPPL220 
SUPPL221 
SUPOL222 
SUPP L 22 3 
SUPPL224 
SUPPL225 
SUPPL226 
SUPPL227 
SUP® L228 
SUP® L229 
SUPPL230 
SUPPL231 
SUPPL232 
SUPPL233 
SUPPL234 
SUP® L23 5 
SUPPL236 
SUPPL237 
SUPP 123 8 
SUPPL239 
S( 1 ) ) ) /DENB SUPP L 24 0 

1.2) *I3ETA(2SUPPL241 

SUPPL242 
3 , 3 ) -QMS ( 1) SUPPL243 
SUPPL244 
SUPPL245 

(3.2) )/DEN3SUPPL246 
A( 2, 2 )-0 M $( SUPPL247 

SUPPL248 

3. 3 ) -OMS ( 2 ) SUPP L 24 9 

SUPPL250 

SUPPL251 

( 2. 3 ) )/0EN3 SUPP L252 
A ( 3 , 3 )~0M$( SUPPL25 3 

SUPPL254 
2)*BETAt 3,2SUPPL255 
SUPP L2 56 
SUPP L2 57 
ALPHA!!, I ) * SUPP L25 8 


GO TO 381 


SORT ( 2 ) *0* 

50ATOF«T* 


64 


361 on 3 84 1=1,3 
DO 384 K*1 ,3 

3 84 DELTA t I , Kl* ALPHA (l ,K) 

CALL ALSOL (3 ,DELT A , SORT ,3 ) 
on 431 1 =1 1 3 

431 GAMMA (ltJ!=S ORT C 1 ) 

432 CONTINUE 


SUP 3 L ? 72 
S U 3 3 L 2 7 3 
SUPP L2 74 
SUPPt 275 
SUPP L 2 76 
SUP 3 L277 
SUPPL278 


WR I TF ( 6 , ll ) 
WRITE(6, 12) 


(I f G AMMA ( 1 » 1 ) , GAMMA ( 1,2), GAMMA! 1,3) , 1 = 1,3) 


SUP 3 L 27 9 
SUPPL280 


Cl. - 


AMPLU - XMUAVR * 

SA = SMALS * S1LB + RM 
SB = SMALS * S 1 L R**2 ♦ 


1*3 133 3 3 33 333DOSUPPL2R4 
SUPPL2B5 


DEL ( 1 


»1) 

* 


DELI 1,2) 
DEL ( l , 3 ) 


* XMUVB *"(l. 

RRDBR * F Ml l ) 
= 2. * SLAMZ 

* Al * (SLAMZ 


R0VB**3 ) / (1, - ROV B**4 ) * 

* S2LB * HIGHS 

RM * S 2L B **2 * HIGHS 

- ROVB**4) / (4. * (l* - SLAMZ * X!3\3**2) 


X BAB 
X BAB 


DELCl, l) 
SB - SA ) 


/ ( 1. - SLAMZ * 


A ♦ B * HIGHS * S2LB 
DEL ( 2 , 1 ) * A2 / Al * DEL(1,1) 

DEL (2, 2) = A2 / Al * DEL (1,2 ) 
DEL(2,3^= A2 _* (SLAMZ * X BAB 
A - R" * SMALS V S2LB 
DEL (3,1) * - SLAMZ * XBA8 * DEL(l,l) 

DELI 1, 1) 


♦ SB - SA) / (l* - SLAMZ * 


/ Ai 


DEL (3,3) = 

A (1. 

C MPA ( 2 ) * CMPAVB 


-2* * SLAMZ 
(BOBR / RRDBR ) **2 
SLAMZ * XBAB**2) 


/ Al 
SLAMZ * 


( X BAB * SA - SP ) / 


CM2) ’ = CL VB 

NDIMC* 63 

CMpTT 


SUPPL286 
SUPPL287 
SUPP L288 
SUPPL289 
X3AB **2) SUPP L290 
SUPPL291 
SUPPL292 
SUPPL293 
X 3 AB** 2 ) SUPPL294 
SUPPL295 
SUPPL296 
SUPP L? 97 
SUPPL298 
SUPPL299 
MAIN 
MAIN 


SI NP S I * 
TO = AT 0 VB 
T OT ( 1 ) * 
DO 


T7 

0. 


4- ATCVB * 
TO - ATOVB 


COS PS I + ATSV B * SIN PS l 


50 


50 1*1,3 
SMALLGC I ) * 


DEL ( I , 1) * CLVB ♦ DELC I, 2 ) * CMPAVB 


DO 51 1*1,3 

GCAP(I ,11*0. 


on 


52 


51 


52 J^TfT 

Y PR ( I , J) = 


0. 


GCAPU ,1) * GCA PiiflT 

GCA P (I , 2 ) * GCAPU, l) 

gcapu ,1 1 / cMsrrr 


♦ ALPHA! I, J) * SMALLG(J) 


YTT7TT 
Y(I ,2) 
T 


TV 

WRITE! 


YU,t) 
PL 0 TOP .LfT 
6,9000) TO, 


9 CUT 


071 

Z, TOPR, ZPR, Y, YPR, DEL, SMALLG 

F0RMATT77* TO* 1 , tPIEI5.'S7'^~Z^ r ’, WU.6,» TOPR=«, IP1E13.6 

* ZPR* f , 1P3E13.6 / • Y*« , lP9E13*6/ t YPR* • , 1P9EL3.6 

1P9 r 13 • 6/ / I 


7 7 1 — DFC* 
RETURN 


» , IP9E13;67 _1 — SMALLG* *T 


ETslTRY ”5 UPPrCOrNFT 


CMPA ( 3 ) * 


CMPAf 2 ) 


CMPA ( 2 ) 
CMPAVB 


SUPPL300 

SUPPL301 

SUPPL302 

SUPPL303 

SUPPL304 

SUPPL305 

SUPPL306 

MAIN 
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CMPA ( 3 ) 


C MPA { l ) = 2. * CMP A (2 > - CMPA<3) 

C L ( 3 ) = CL ( 2 ) 

CLI2) * CL VB 

CLUJ = 2 . ♦ CL (2) - CL I 3 ) 

PSI = I BDRR / RRDBR ) * NT I ME * 
SIN PSI a S I N I P S I ) 

cos psi a cosipsii 

TOT I 2 ) = TOT II) 

TO = ATOVB ♦ ATCVB ♦ COS PS I + 

TO TH)* TO - ATOV B 

TO PR = IBDBR/ RRDBR) • IATSVB * 

00 60 K * It 2, 

00 64 I a l, i 

SMALL G UI = UINF *»2 ♦ IDELIl, II 
A «■ DELI 1 ,3) * TOT (K) 

DO 65 I a i f 3 

G CAP ( I t K) = 0, 


ATSVB 


SIN PSI 


ATCVB 


CLIK) 


DELI 1,2) 


CMPAIKl 1 


DO 65 I a i f 3 

G CAP ( I t K) = 0, 

DO 65 J * It 3 

G CAPII, K) * GCAPUt K) ♦ ALPHA! I » J 1 ♦ SMALLGIJ) 

C ONTINUE 

DO 62 I * It 3 

YU t 2> a YII ,1) 

YPR (I f 2 ) = YPR 1 1 «1) 

WD X I a OMEGA (I) ♦ DX I 
SWOXI a SINIWDXI) 

C WD XI a COS(MDXn 

Ydtll = Yd *21 * CWOXI ♦ YPRI I » 2 ) * SWDXI /OMEGA! I) 

A fr I IGCAPI I »21 - GCAPUtlD ♦ (SWOX I - WDXI ♦ CWDX II / 

B ♦ GCAPIItll * (1 • - CWDXI11 / OMEGA! I l**2 
YPR(I t l) a YPR ( I ,2) ♦ CWOXI - OMEGA!!) * YII f 2) * SWDXI 


A ♦ ! (GC API l 1 2) - GCAPIItll I ♦ (WDX 

B / W0X1 ♦ GCAPIItll ♦ SWDXI) / OM 

DC 61 I » It 3 

zm » o. 

ZPRII) a 0, 

00 61 J * It 3 

ZII) = Z(II ♦ GAMMA I 1 1 J) * Y ( j , 1 1 
ZPRII) = ZPRII) ♦ GAMMA 1 1 1 J) * YPRIJ, 1) 


I WDXI ♦ SWDX I 
/ OMEGA! I ) 


CWDX I 


ALPH1 * TO «- ZC3I 
ALPH2 * TO PR ♦ ZPRI3) 

HEAVE « - ZPRII) ZPR12) 

IF I PLOTOP .LT. 0.) 

1 WRITE I 6 1 9000) TO, Z, TOPR, ZPRt Y, YPR, DEL , SMALL G 

2 , TOT __ 

RETURN 

FORMATI5Ei0.4) 

FORMATrSF i 0*6) 

F0RMATI1H1 ,10Xt* ITERATION FOR XBAR DIVERGED*) 

FORMAtriHl , 5 X, 4 HE 1 »E13. 5 , 5X , 6HF2 *E f3.5t 5X 4 HE 3 '« E 13. 3//5X t 4 HRM 

IE13.5////I 

EORMTnTx tTHMTI^*ETT TSTSXt 5H M 2 2“~*Fr3 75t CTt ^HM3 3 * ET 13 • 5 « 5X f 5HM 13 "« 
1 13. 5 »5X,5HM23 -E13.5/I 


SUPPL308 

SUPPL309 

MAIN 

SUPPL3L1 

SUPPL312 

SUPPL313 

SUPPL314 

SUPPL315 

SUPPL31 6 

SUPPL31 7 


SUPPL32 2 


113. 5 t5X,5HT23 -E13.5///I 


*Tl3.5tidX 


SUPPL360 

t6SUPPL361 
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At 
A A 

A 6 
A 7 


A9 

11 

12 
A 88 


A69 

721 


1H11/B -E13 *5/20X t6HL2/R *E 13 . 5 , 10X ,6HL2/B *E 13 . 5/9X t 7HK 
1 /9X , 7HK2 /M2 *E13.5) 

FCRMATC //i0X,5HR/R =E13 . 5 ,20X ,6HRR/R =E13*5///) 

FORMAT ( l HI , 2 0 X , • POL YNC1MI At COEFF 1C 1ENT S 1 / //7X f 5HPOWER,! 
1 26Xt 3H2-D7T 

FORMAT ( l 10 *2030*9) 

FORMATU HI , 2 OX , * R HOT SHF POLYNOMIALS* ///30X, * BL AOF # ,60X 
lAHREAL,2iX,AHIMAG,3lX,AHREAL, 21X,AHtMAG/) 

FORMA T ('2*02 5.9,1 OX §2 025 • 91 

FORMAT ( /////9X f 1HI , 15X, l OHGAMMA (I f II, 15X, 10HGAMMAI 1,21, 
lA II ,3) /) 

FORMATU 10 ,3E25# 5 ) 

F OR M A T(i H 1 , 8 X , l H I , 7 X , 5H 0 MEGA , AX , 9 HBET At I, 1), AX, 9HBETA ( I 
1TAU ,3) , 3 X , 1 OH ALPHA f I ,11 ,3X, l OH ALPHA! 1,21, 3X, 10H&LP FA( I 
ik//i 

FORMATU 10,8E13.5) 

“ F 0RMATT7/7I 0 X , 5'ff F k l *£ 1 3V5 ,‘t ff X , $ HF R2 *E13.5, 10X/5HFR3 = 
1AHSA *E13. 5,10X,AHSB =E13.5///I 

END' 


mi *F13.5SUPPL362 
SUPPL363 
SUPPL36A 
2X, 5H3LADE , SUP 3 t 36 5 
SUPPL366 
SUPPL367 
, , 2-T*/20X,SUPPL36fl 
SUPPL369 
SUPP L 370 
15X, 10HG AMMSUPP L 371 
SUP 3 L 372 
SUPP L373 
, 2 ) , AX , 9HBE SUPP L 37 A 
, 3), 3X,3HCHSUPPL375 
SUPPL376 
SUPPL377 
E 13 • 5//10X, SUPPL378 
SUPPL379 
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r 

c 


SUBROUTINE SETUPS 
IMPLICIT REAL *8 !A-H,0-Z) 


REAL FTVBt _ FPVfl f __ FPPRVB, DIDRVB, XNV 8, 

A FOVB, XMUAVB, ATOVB, ATCVB, ATSVB, 

REAL ELSIGt OXI , REB f RDBB, FRZ , ARR, 

ALPH2 * HEAVE, AROT, FREQF, PHIH. 

XL, YL, ER L, ER2, ER 3, RDBR, 


A 

B 

C 

H, 
I f 


ALPH1 
Y, TEST, 


DELVR, XMUVB, 
ROVR, RV8, MVB 
AMPLU, FREQU, 

RYU ORY, 


L 

2 

3 

4 

5 

6 
7 


UPRl M , XU, Y J i 


RROBR 

C MPA , CMPAS, BARG , E Ml f HVOR, 
PLOTOP, P SI LOW , PS I UP 


SSPA, SVOR, TURF, XlVDR 


INTEGER TABLE <7, 80) /560 * • •/ 


SETUPS 
SETUPS 
SFTUPS 
SETUPS 
SETUPS 
SFTUPS 
SETUPS 
SETUPS 8 
SETUPS 9 
SE TUP SI 0 
SE TU P SI l 


SE TUP SI 2 


f 








SETUPSI 5 


COMMON /BLi / 

NTIME 





SETUP SI 6 

c 








SETUP SI 7 


COMMON /INPTVB/ 

F T VB! 64 ) , FPVB ( 64 1 

, FPPRVB! 64 ) , 0 I ORVB ! 64) , 

SETUPS18 


A 

XMVBI64) , 

DELVB , 

XMUVB, 

FOVB, 

XMUAVB, 


SETUPS19 


B 

ATOVB, 

ATCVB, 

ATSVB, 

ROVR 

, RVBI64), 

SETUPS2 0 


C 

MVB ( 64 ) , 

NVB 





SETUPS2 1 

c 








SE TUPS22 


ccmmcn /inputs/ 

"HTBL, 

NZ, 

NOFF, 

NGAM f 

NSIG, 

SETUPS23 


A 

Ncm , 

NCORO, 

LOWER, 

MSTOP, 

MAXT, 

MQTR, 

SETUPS24 


B 

NOTBL , 

! NOV, 

ELS 1 G, 

OXI, 

RER, 

ROBB, 

SETUPS2 5 


C 

FRZ, 

ARR, 

AMPLU, 

FREQU, 

ALPHl, 

ALPH2, 

SETUPS26 


D 

HEAVE, 

AROT, 

FREQF, 

PHIH, 

NY, 

RY1, 

SETUPS27 


E 

DRY, 

Y ! 100 ) , 

TEST, 

UPRIM, 

XU! 30), 

Y U ( 30 I 

, SETUPS28 


F 

XL ! 30 ) , 

YL 133 ) , 

ER1 , 

ER2, 

ER 3, 

BOBR, 

SETUPS29 


G 

RRDBR 






SETUPS30 


H, CMPA, CMPAS* BARG, EMi, HVOR, NVOR, SSPA, 

SVOR, TORF, 

X IVOR 



I, PLOTOP, PSILOW, 

PSIUP 







J , NCUT 







c 








SETUPS31 

c 








SETUPS32 

c 








SETUPS33 


CALL 

WHERE (TABLE) 






SETUPS34 


CALL 

ZEROIN 






SETUP S3 5 

c 








“3FWPS36 









SETUP S3 7 


CALL 

SETUP! • ALPHl 

',4, 

ALPHl 

) 



SETUPS38 


CALL 

SETUP! *ALPHAl 1 ,4 , 

ALPHl 




SETUPS3 9 


CALL 

SETUP! • ALPH2 

' ,4 . 

ALPH2 

T 



SETUPS40 


CALL 

SETUP!' ALPHA2 ',4, 

ALPH2 

) 



SETUPS4 1 


caLL 

SETUP! * AMPLU 

• ,4, 

A MPLU 

) 



SETUPS42 


CALL 

SETUP! • ARR 

* t4. 

ARR 

) 



SETUPS43 


CALL 

SETUP! * AROT 

•,v. 

AROT 

i 



SETUPS44 


CALL 

SETUP!' ATOVB 

* f4. 

ATOVB 

) 



SETUPS45 


CALL 

SETUPl'ATCVB 

* t4t 

TfTCVB 




SETUPS46 


CALL 

SETUP!* ATSVB 

',4, 

ATSVB 

» 



SETUP S47 


CALL 

SETUP ! * BARG 


ft ARG 






CALL 

SETUP!* BOBR 

' ,4 , 

BDBR 

) 



SETUPS48 


CALL 

SETUP! 'CMPA 

f ,4, 

“CWPA " 

— j— 
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CALL 

SF TUP * 

CMP1\S 

• ,4 , CMP AS ) 


CALL 

SC T UP t 

DFLVB 

',4, OELVB ) 

SETUPS49 

CALL 

SETUP* 

DT DRVB 

*,4, OIDRV8* 64 ) 

SETUPSSO 

CALI 

SF TUP ( 

DRY 

• ?4 f ORY ) 

SETUPS51 

CALL 

SF TUP ( 

DXI 

• ,4, nxi ) 

SETUPS52 

CALL 

SE TUP * 

FLSIG 

• *4 » ELSIG ) 

SFTJPS53 

CALL 

SE TUP ( 

EMI 

•,4, EMI ) 


CALL 

SE TUP ( 

ER1 

* t4 » ERI ) 

SE TUP $54 

CALL 

SF TUP ( 

ER2 

* 1 4 1 FR2 ) 

SFTUPS55 

CALL 

SETUP ( 

ER3 

* t4 ♦ ER3 ) 

SE TU° S56 

CALL 

SE TUP ( 

FPVB 

• »4 » FPVBt 64 ) 

SFTJPS57 

CALL 

SE TUP ( 

FPPRVB 

# *4 » FPPRVB, 64 ) 

SF TU P S5 8 

CALL 

SE TUP ( 

FRZ 

* »4 ? FRZ > 

SETUPS59 

CALL 

SETUP * 

FREQU 

* ,4, FREQU ) 

SFTUPS60 

CALL 

SETUP* 

FREQF 

• 1 4 t FREQF ) 

SETUPS61 

CALL 

SE T UP ( 

FT VB 

• ,4, ETVBt 64 ) 

SETUPS62 

CALL 

SE T UP 1 

FOVft 

_r »47 FOVB ) 

SETUPS63 

call 

CALL 

SE TUP C 
SETUP ( 

HEAVE 

HVOR 

• f 4, HEAVE ) 

1 »4 t HVOR ) 

SETUP S64 

CALL 

SE TUP ( 

INDV 

',4, INDV ) 

SETUPS6B 

CALL 

SETUP ( 

LOWER 

' )4» LOWER ) 

SETUPS66 

CALL 

SE TUP * 

MAXT 

',4, MAXT ) 

SETUP S6 7 

CTOT 

sf tup ( 

MOT R 

',4, MCTR ) 

SFTUPS6 8 

CALL 

SE TUP * 

MSTOP 

• 1 4 * >STOP ) 

SETUPS6 9 

TALI 

SF TUP( 


^T4i MVB , 64 > 

SETJPS70 

CALL 

SE TUP * 

NCOI 

• t4 t NCOI ) 

SETUPS? I 

CALL 

SETUP ( 

NCORD 

• »4 » NCORD ) 

SETUP ST2 

CALL 

SE TUP* 

NGAM 

* ?4 * NGAM ) 

SETUPS73 

CALL 

SETUP! 

NOFF 

■ t4 1 NOF^ ) 

SETUP S74 

CALL 

SFTUPI 

NOT 0L 

•,4, NCTBL ) 

SF TU P $75 

CALL SE TUP * * NQUT 

• , 4» NO'JT j 


CALL 

SETUP ( 

• NS8L 

• »4 » KSBL ) 

SETUPS76 

CATlT 

SE'TUPT 

• NSIG 

• t4 t NS I G ) 

SETUP S77 

CALL 

SETUP* 

• NVB 

•, 4, NVBI 

SFTUPS78 

CALL 

CALL 

SETUP* 

SETUP* 

• NVOR 
•NY 

• *4 » NVOR ) 

•*4, NY ) 

SETUPS79 

C~AIT“ 

SETUP! 


',4, NZ ) 

SETUPS80 

CALL 

SETUP* 

•PHIH 

',4, PHIH ) 

SETUPS8I 

TSIT 

SETUP! 

* PLOT Op 

• , 4, PLOTOPJ 


CALL 

SETUP* 

•PSILOW 

• , 4, PS ILOW ) 


CALL 

SE t UP ( 

f PSIUP 

_r , 4, PSIUP) 


CALL 

SETUP* 

• RVB 

• , 4, RVB, 64 » 

SETUPS82 

CALL 

SETUP! 

• ROBB 

’,4, ROBB T 

, SETUPS83 

CALL 

SETUP* 

•RFB 

',4, REB ) 

SETUPS84 

call 

SETUP! 

•RRD&ft 

• f 4 y RRDSR T 

SETUPS85 

CALL 

SETUP* 

•ROVB 

• ,4 » ROVB ) 

SETUPS86 

CALL 

SETUP! 

• RVl 

• ,4 , RYl r 

SETUPS87 

CALL 

SETUP* 

• SSP A 

',4, SSPA I 


CALTT 

CALL 

SF TUP I 
SETUP! 

' SVOR 
•TEST 

' ,4," 5VUR ' T 

• ,4, TEST ) 

SETUPS88 

CAL I 

“SETUP! 

■ TORE 

• ,4 , TORT T 


CALL 

SETUP* 

•UPRIM 

',4, UPRIM ) 

SETUPS89 

CALL 

CALL 

st Tup ! 
SFTUPI 

• xivOr 

' XL 

' ,4, Xl VCR J 

',4, XL, 30 ) 

SETUPS90 

CTO7 

SETUP* 

■ XMVB 

• ,4 , XMVB, 64 1 

SETUPS91 
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CALL SETUP! 'XMUVR • ,4 , XMUVB ) SETUPS92 

CALL SFTUP! • XMUAVB • ,4 , XMJAVB 1 SETUPS93 

CALL SETUP ( ' XU ',4, XU, 30 ) SETUP S94 

CALL S ETUP!' Y _IOO ) SETUPS95 

CALL SETUP CYL • *4 , YL, 30 ) SETUPS96 

CALL SETUP.** YU _*,4, YJ, 30 ) SETUPS97 

C SETJPS98 

C SE TUPS99 

C SETUP100 

C . . _ _ _ _ SFTUPIOI 

C SETUP 102 

C SETUP103 

C SETUP104 

PSILOW* 1 « El 0 _ 

PSIUP- —l • El 0~ 

PLOTOP « 1. 

NOUT * 0 " 

RETURN SETUP105 

C SETUP106 

END 
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SUBPOUT l NF BLCTX ,Y, MST . MEND, NY, RY , DRY , DX I , RER, JPP IM f PL AM f X FL AM, TESftlC 
1 T » U , SCALF, HE ,UC , V,XSE*, US EP, D IS S * THET S , L OW ER , t AMQ, M S BP, XC, J SA V, SCARCE 
ILSt NITS, NT I VE, NCTBL, XTEST, NZ , NCUT ) 


PROGRAM FOR ANALYZING LAMINAR AND TURBULENT BOUNPARY LAYERS 
BY THE METHOD OF FINITE DIFFERENCES. IF THE INTEGER L AMQ 
IS GREATFR THAN ZERC, THE BOUNDARY LAYER IS LAMINAR. 

COMMON /RLI / N DUMMY 9 NDIMC , ISTD 

DIMENSION U$AV(300,100) ,SCAL$<300 ) 

DIMENSION X (300 ) ,Y<100) , UE ( 30 0, 3 I , UC ( 100, 3), VI 100, 2), XC ( 300) 
DIMENSION SD(iOO) ,SE(i001 f SF(100),VISC( 100,2), GRAD! 100) 
DIMENSION A (100 ) ,8(100) ,0(100 ),0( 100) , Ft 100) 

DIMENSION ALPHA (100) , BET A( 100 ), GAMMA ( 100 ) , DELTA! 100) 

DIMENSION SC ALE! 300,2) ,V ARK 100 ),VAR2( 100 ) 

DIMENSION FLAM(LO) , XFLAM < 10 ) , Yf* 1 ( 100 ) ,YB2( 100) 

DIMENSION U(300 A >0,21 

DIMENSION CAPGUOO) ,CAPH( 100 ) , CAP J t 100 ) , C APK ( 10 0 ) 

DOUBLE PREC I SIGN API 100) , BP! 100 ),CPl 100) , DPI 100 ), Fp( 100 ),UP ( l 
FORMAT! IHl ,4iX,36H ANALYSIS OF LAMINAR BOUN PAR Y L A Y ER/ // 5 U , 
l ME STE P NO! 3/75 1 X ,1 2 HI TER ATI ON NOT 3/ / / 4X , 1HM , 8X , iHX, L3X ,2HXC, 
IHUE ,10X,6H-DP/DX,9X,5HDELTA,9X,5HDISPL,9X, 5HTHETA, 9X, 5H SHEAR / 
FURTTSTTlHl ,4lX ,36 HA^LYsTS OF YURBUL £ NT BOUND AR Y LAYER/ //5i X , 
1ME STEP N0I 3//51X,12HIT ERATI0N NO 13// / 4X , 1HM , 8X , IHX, 1 3X , 2HXC , 
1H UE ,T0XV6H^DP/D X ,9 X ,5HCTITA, 9 X, SHOTS PL , 9X , 5HTHE T A, 9 X , 5H S HEA R , 
3 MV) 

F ORMA T ( I 5 ,8 £ 1 *. 4,1 3 ) 

FORMATdHl ,2X,3HM *I4//2X,3HX *E14.5//2X, 4Hl)E *E14.5,i 
IHDXDP/DX) =Fr4.5~,lOX,5H*EB *E14 .5 , 19X, 4HU T *E14.5///) 
F0RM AT!2X,25HPHYSICAL DELTA = E14 . 5, 8X , 12HDELTA 

1LTA STAR *E14« 5 ,8X,7HTHET A *£14.5///) 

FORMAT! 25 X,lHY,i9X,lHU,19X,lH\r f TOT, 5HDU/ DY , 1 4X , 6HNUE/NU 
F0RMAT!10X,5E20.5) 

FDR MATT / /30X,17HSEPARAT ION AT X *E13.T,“6hT XC « £ 13 * 5 ) 
F0RMAT(///40X,12HWALL SHEAR *E14.5//) 

F ORMAT I /750XiT7HTR4NSlT I CN AT X =F14 .5) 

F ORMAT ( /20X,35HSCALE CHANGE - Y-MAX INCREASED FR0ME12.4 

in 

F0RMATI10X»7HAT STEPI3,22H, THE WALL GRADIENT ISE12.4I 
BC ON = 


10 


ii 


12 

20 

24 


21 

22 

23 

25 

30 

35 

810 


rl: 

BLC 

BLC 

RLC 

BLC 

BLC 

BLC 

BLC 


90) B 
12HTIB 
12 X ,2R 
) B 
12HTI ft 
12X ,2B 
4X , B 


I.b/UXI 

FCON * l./(2.*DXl) 
TFT I STD.NE. 

D XI * 1 . E 30 


OX, 1 7H-! 1/RBLC 
BLC 

STAR =E 14. BLC 
• 5, 8 X , 12HDEBLC 
BLC 

/) BLC 

bl: 

BLC 

BLC 

BLC 

3H TQE12.48LC 
BLC 
BLC 

bl: 

BLC 


1) GO TO 900 


FCON^O. 


TCTJ 


nTNTrmJF - 

MOUT =6 

MTRATFr*^! 

Y SUB 2= Y! 2 ) 

W5T2 * 


MSTT 


MST - 
MST-1 


BLC 

BLC 

BLC 

BLC 

BLC 


NOIITI *“ 
MST1 MO* 


WOT 

MOD! 


n 

MSTl 


NOUTi) 




i 

? 

4 

5 

6 

7 

8 


10 
1 i 
12 
1 3 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

51 

52 


39 

40 

4 1 

42 

43 
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o no 


543 

544 

545 
550 


GO TO (543,550) ,11WER 
If (LAMQ) 544,544,545 
WRI TF(M0UT,lI) NTIMF ,NITS 
GO TP 550 

WRITE (M0UT,10) NT I ME , NI T S 
CONTINUE 
YTR = SORT ( REB) 

UC (I ,1) = 0. 

V(lfl) * 0. 

N V = NY 2 

NVMl = N V - l 
NVPl « NV ♦ 1 

CALL YD IFF (NY, ALPHA, 8ET A, GAMMA, DELTA, SD, S E, SF, C 2, C3, C4, Y ) 
DO 41 N-l , NVPJ. 

visc(N,n * i. 

V I SC ( N ♦ 2 1 = U 

DC 42 M = MST2 , MS Ti 
L a MSTl -M4* 2 
DO 50 N*1 , NV 

5 0 GRAD(Nfj) * SP(N+ U *UCJ_N*2 ,L! *S E( N*1 ) *UC(N«-1 ,L l-SFI N*l> *1 JC 
GRAD ( X ) = C2*UC(2,L> +C3*UC(3 ,L)+C4*UC(4,L) 

MM=M— 1 


41 


CALL PGRAD(MM,X,UE,DXI , PR ESS , S A , S B, SC , SR, S S ) 

DO 456 N*l ,NY 

4 56 UC ( N , 1 ) * UC (N,L) 

CALL SET I T (LAMQ,M,NV, RE B,X,Y,UCt PRESS, GRAD, Dl 
IRAN) 

42 CONTINUE 


MENDi a MEND - 1 
GRADS*GR AOf i ) 
GRAOSS*GRAD C l ) 


THE MAIN CALCULATION STARTS HERE. 


DO 99 M: 
I TER*0 


MSTl , MENDI 


WALLG^O. 
MP1*M+ l 
DELTP * 
DISPT * 


"wrr 

SHEAR 


DE LT/ YTR 
DI SP*YTR 

THETA*Ytft 

GRA0<1)/YTR 


MSTl MD) GO TO 225 


3“6X 


IF ( MODIM, NOUTU.NE. 

GO TO (561 ,562) , LOWER 
WRITE(MOUt ,12) R,X< Ml, XC(M),UE(M, l ), PRESS, DELTP, DISP, THETA, SHEAR BLC 




BLC 

44 


BLC 

45 


BLC 

46 


BLC 

47 


BLC 

43 


BLC 

49 


BLC 

50 


BLC 

53 


BLC 

54 


BLC 

55 


BLC 

56 


BLC 

57 


BLC 

58 


BLC 

59 


BLC 

60 


BLC 

61 


BLC 

62 


BLC 

63 


BLC 

64 

L) 

BLC 

65 


BLC 

66 


BLC 

67 


BLC 

68 


BLC 

69 


BLC 

70 

SC.MTBL3 

71 


BLC 

72 


BLC 

73 


BLC 

74 


BLC 

75 


BLC 

76 


BLC 

77 


BLC 

78 


BLC 

79 


BLC 

80 


BLC 

Bl 


BLC 

82 


BLC 

83 


BLC 

84 


BLC 

85 


BLC 

86 


BLC 

87 


BLC 

88 

(EAR 

BLC 

89 

— 

BLC 

90 


BLC 

91 


BLC 

92 


BLC 

93 

»i> 

BLC 

94 


BLC 

95 


"IFWSX-XfHIT 


xnr 

BLC 

BLC 


"96T 

97 

98 


72 


u o o 


4C5 


QLC 

OLC 


99 

100 


WF S=< XSX-XIM-l) I / IXIM)-Xl M-l ) ) 

GO TO 2?.<t 

229 IF ( GRADU) ) 227, 227, 273 

273 IF (OISP .GT. 0. • A Nil. THFTA .GT. 0.) GO TO 2?3 

2 83 CONTINUE 

XSFP* XC (M-l) 

U SE P = UE (M-l ,1) 

XRL=XI M-l) 

WRITE IMOUT, 23) XBL , XSEP 
RF TURN 

2 2 7 WE S=GRAO S7 (GRAOS-GR AO(l ) ) 

2 24 WF SI =) . -WF S 

XSEP=WFS1*XC(M-1 ) fWFSOXCCM) 

XRL = WF SI * X ( M— l )-fWFS*X(M) 

USEP=WFSl*UE (M-l ,1 ) +WFS *UE I M* 1) 

WFP=(XRL-XIM-2) I /IX I M-l )— X(M-2) ) 

WFPi=i.-WFP 

OISS*Ol SSS*WFPl+OISS*WFP 

THE TS=THE TSS*WFPl*-THETS *WFP 

WRITE(MOUT,23) XBL, XSEP 

IF ( LAMO.EQ.O. ANO.M. Lf . MTR AN+5 » L AMQ= l 

GO TO 222 

223 CONTINUE ; ' 

IF! N0T8L .EQ. 2 .AND. NITS .GT. 1 .ANO. M. GT. NZ .AND. 
1 Xf. IP) .GT. XTEST) "GO TO 283 
IF(LAMQ) 801,801,302 
8C2 IFIN0T8L . EQ. 2) GO TO 801 

CALL TRANS ( IJPRI M, PRESS, THETA, RED, UC, NY, FL AM, XFL AM, L AMO) 

1 FITAMQT 805,805,301 

8C5 WRI TE ( MOUT ,30) X(M> 

MtRAN * Mvl 

801 CONTINUE 

tFTYfNV) -15FLTV 620,641,641 

620 R Y=R Y+ORY 


RE SC AL I N G CALCULATION STARTS HERE. 

00 632 N=l , NY 
YB 1 (NT = Y(N) 

VAR l (N) * UCIN,2) 

632 WR2ITCr * UC(N,3J ~ 

CALL Y SE T ( R Y , YSU 82 , NY ,Y ) 

WRITOWTUT735 )“" Yll ( NYT ,YTNYF 

00 633 N=2,NVPl 

YIN * YIN) 

CALL TERPI YIN,YBl,VARl,NY,UPASl) 

UC (N,2Y‘ = UPASl 

CALL TERPI YIN, YB1,VAR2, NY, UPAS2) 

633 UC (NT, 3) = UP AS2 

CALL YDI FF ( NY, ALPHA , BET A, GAMMA, OELTA, SO, SE,SF,C2,C3,CA, Y ) 

1 FIX AWT) 700,700,701 — — — 

700 DO 635 N*2 , NVPl 

vsRKwr * wscinyi) 

635 VAR2IN) = VISC(N,2) 

0Tr&T6“N=2,NVPl 


BLC 

LO? 

blc 

1 03 

blc 

104 

BLC 

105 

BLC 

1 06 

BLC 

107 

BLC 

108 

BLC 

109 

BLC 

110 

BLC 

111 

BLC 

112 

BLC 

113 

BLC 

114 


BLC 

115 

BLC 

116 

BLC 

117 

BLC 

118 

bl: 

119 

BLC 

120 

BLC 

121 

BLC 

122 

BLC 

123 

BLC 

124 

BLC 

125 

BLC 

126 

BLC 

127 

BLC 

128 

BLC 

129 

BLC 

130 

BLC 

L31 

BLC 

132 

BLC 

133 

BLC 

134 

BLC 

135 

BLC 

136 

BLC 

137 

BLC 

138 

BLC 

139 

BLC 

140 

BLC 

141 

BLC 

142 

BLC 

143 


73 






YIN = YIN) 

bl: 

144 


CALL TERPI YIN.YBl , V ARl , NV PI , U PAS 1 1 

BLC 

145 


VISCIN.l) = UPASl 

bl: 

146 


CALL TER PI YIN ,Y31 , V AR2 , NV PI , U P AS 2 I 

BLC 

147 

636 

VI SC IN, 2) = UPAS2 

BLC 

14 fl 

7C1 

DC 637 N=*2 tNVPl 

BLC 

149 


VARIIN) = VIN.l) 

bl: 

150 

637 

VAR2IN) « VI N 1 2 ) 

BLC 

151 


00 638 N*2,NVP1 

BLC 

152 


YIN * Y I N) 

BLC 

153 


CAtL TERPI YIN.YBl , V ARl , NV PI , U PAS 11 

BLC 

154 


VINtl) = UPASl 

BLC 

155 


CALL TERPI YIN.YBl, VAR2.NV PI, UPAS 2) 

BLC 

156 

638 

V (Nt 2) = UPAS2 

bl: 

157 

641 

CONTINUE 

rl; 

158 

C 


BLC 

159 

C RESCALING CALCULATION EADS HERE. 

BLC 

160 

c 


BLC 

161 


CALL PGRADIM.X.UE.DXI .PRESS , S A, SB, SC, SR, SS ) 

BLC 

162 

c 


bl: 

163 

C RECURSION RELATIONS ARE SET UP HERE. 

bl: 

164 

c 


bl: 

165 


IF IISTO.EQ. 1) GO TC 820 




IF(SCALECM+l«l)-U) 522,522*521 

BLC 

166 

5 IT 

If f ISCALF(M+-l,2)-i.) 522 ,522, 523 

BLC 

16 7 

522 

LACKU=1 

rl: 

168 


F ACU1*UE I M*l ,2>/UE(M+l,il 

BLC 

169 


FAC U2*UE I M*l ,3)/UE(M*l,l) 

BLC 

170 


GO YE) A20 

BLC 

171 

523 

LACKU*2 

BLC 

172 


nTTSTO NN*1 »NY 

BLC 

173 


VARl INN> » U( M*l ,NN ,1 ) 

BLC 

174 

6l0" 

VAR2INNI « U( M+l ,NN,2) 

BLC 

175 


CALL YSETISCALEIPU ,1) , YSUB2 , NY , Y Bl» 

BLC 

176 


CALL YSE T ( SCALE I M+l ,2 1 , YSUB2 , NY , Y B2) 

BLC 

177 

820 

DO 88 N*2 » NV 

BLC 

178 


CALL CAPS! ITER, N.CAPG.CAPH/CAPJ, CAP*. s*, SS, SO,S e, SF , V ISC , V, UC I 

BLC 

179 


A IN) *-SF I Nl *C APG IN)— DELTAIN)*CAPH(N)+SF(N)*CAPJ(N) 

BLC 

180 


' B ( N) *BCON«- SA*CAPk(TD *SF I N I *C AP6(“RV-6AWWif NI*CAP HI N 1-SfI N )*CAP J(N I 

BLC 

181 


C INI -SOI N) *C APG ( NI— BETA! NI*CA PHI NI-SO(NI*CAPJ{ M 1 

BLC 

182 


OIN)~ -ALPHA(N) *CAPHTN) 

BLC 

183 


IF (I STD .EQ* i) GO TO 576 




GO TO 1574,575) ,LACKU ~ 

BLC 

184 

574 

UPAS1*FACU1*UC(N,1) 

BLC 

185 


UPAS2«FACU2*UCIN,1J 

T) LC 

186 


GO TO 576 

bl: 

187 

5TT 

YIN » VlN) — 

Err 

108 


CALL TERPI YIN.YBl, VARl, NY, UPASl) 

BLC 

189 


CALL TERPI YIN,YB2,VAR2,NY,UPAS2I _ 

BLC 

190 

516 

FIN) = PRES$*FC0N*14.»UPAS1— UPAS2) ♦CAPKINIPI SB*UCIN , 21- $C*UC (N, 3 1 IBLC 

191 

tnr~ 

CONTINUE 

BLC 

192 

c 


bl: 

193 

c 

SOLUTIONS OR VELOCITY PROFILE STARTS HERE. 

BLw 

“T9£ ' 



bl: 

195 


DO 89 N*Z , NV ~ 

BLC 

196 


.74 


CPC NV)))/RP(NV) 


4 P(N) ■= A ( N) 

RP(N) = B ( N) 

C. P ( N ) = C(N) 

DP(N) = D ( M) 

FP(N> = Fill " 

DO 77 N=2,NVNl 
CP(N) » CP(N)/RP(1) 

DP(N) = DP(N)/BP(N) 

FP(N) = FPINl/BPIN) 

RP(N-t-l) = BP(N+ll - CP(N)*AP(N«-l) 
r.P(N+l) = CP(N*l) - DPIN)*AP(N+l) 

FP(N-H) = F P (N+l I - FP(N)*AP(N«-1) 

UP (NY) = UEIMU ,1) 

UP(NVPl) = UP (NY) 

UP ( N V) * ( F P ( N V I -U P ( NY) * ( DPI NV ) ♦ CP( NV ) ) )/RP(NV) 

DO 6 6 N=3 i NV 

'NN=NW2-N •' 

UP(NN) = FP(NN) - f)P(NN) *UP( NN+2 ) - CP( NN )*UP( NN + 1 ) 

DO 65 N=2 »"NY 

UC ( N 1 1 ) * UP ( Nl 

IF (I TER) 843 , 84 i, 84 3 

DO 8A2 N=2 , NVP1 

vin, 2 ) * v(N,n 

VISC (N,2) = VISC(N,l) 

DfS$S=07T5 

diss=disp 
THE t SS=THETS 
THETS=THETA 


GRAD S=GRAD II) 

DO S5"R=2 , NVP1 

VtN.l) - V ( N— 1 (1 ) -• 5*1 Y (N)-Y( N-l ) I *1 S A*IUCI N» l) *UC( N-l* 
l N, 2T+UC T N-l ,2 )) ♦ SC* (UC I N , 3 ) +UC I N-TV3) H 
DO 56 N=1 , N V 

~5RADIN*-I) "= 5DINM) *UC(N+2,1USE(N+H*UCIN«,'TT-'SF IN ♦ 1 ) 
GRA D (I) = C2»UC(2tl 1)*C»*U C( A, II 

1MTRAN) 


*(JClN,l) 


tter=i ter*i ; 

GO TO (330,809) , LOWER 
WRITE (MIJUI ,8101 l TERTGRAC 
I F ( I TF R— 9) 811,811,812 

FPW=ABS(GRAD(l J-WALLGJ 

I F ( WALLG-1 • ) 120,120,119 

E PW=E PW/WALLG 

IF(EPW-TEST) 812 ,B1A,81A 

WAILG=GRADI l) 

GO TO 820 

DCT 44 N=1 , NY 

UC (N,3) = UC IN ,2 ) 

UC (N,2) = UC ( N , 1 ) 

CONTINUE 


BLC 

19 7 

BLC 

19;} 

BLC 

19 * 

BLC 

209 

BLC 

201 

BLC 

20? 

BLC 

203 

BLC 

204 

BLC 

209 

BLC 

206 

BLC 

207 

BLC 

208 

BLC 

209 

BLC 

210 

BLC 

211 

BLC 

212 

BLC 

213 

BLC 

214 

BLC 

215 

BLC 

216 

BLC 

217 

BLC 

21 8 

BLC 

219 

3LC 

22 0 

BLC 

22 1 

BLC 

222 

BLC 

223 

BLC 

224 

BLC 

22 5 

BLC 

226 

BLC 

227 

BLC 

228 

BLC 

22 9 

BLC 

230 

BLC 

231 

BLC 

232 

) BLC 

233 

BLC 

234 

BLC 

235 

BLC 

236 

BLC 

237 

BLC 

238 

BLC 

239 

BLC 

240 

BLC 

241 

BLC 

242 

BLC 

243 

BLC 

244 

BLC 

245 

BLC 

246 

BLC 

247 


IF ( I STD .EQ. 1) GO TO 


75 



48 liSAVtM+l ,N»*UC(N,1| 



SCALS(M*1>,ry 


BLC 

249 

99 

CONTINUE 


bl: 

250 


X$EP*Ul 


BLC 

251 


U$EP«UEIMX*l) 


BLC 

252 

2 22 

CONTINUE 


BLC 

253 


RETURN 


BLC 

254 


END 
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sup rout i NF’FlnYsBi Piornp , p , l ) 

REAL * 3 ORD(6 1 

DIMENSION" 'P(200,7), TITK56I , NF(5,4) 

I , NF P 1 6 ) 

DATA Nl" , N2 , NO , N42 

1/1 , 2 , 0 , 42 / 

DATA ORO / ' THETA-p* TORS • FLAP-H •, • RENO-H 

1 • n • CM-A • / 

IF ( PLDTOP .EQ. 0. ) RETURN 
I F ( L .LT. 2) RFTURN 

IF ( Piornp .EQ. 2.) GO TO 2 
PLOTOP = 2. 

CALL lOFRMV { • CR I MI -PETE • , • 30', '5100* ) 

2 C INTI NUE 

3 NL =1 

DO l J = 1, 6 

CALL E Z PLOT ( 9 . , Nl , Nl, P » PI 1, J+l ) , L , -N 1 , M2 

1 , N42 * 1^ , 12 , • PSI-DEGREES' , 8 , ORDI J) 

2 , Nl , Nl ", XL ", XU , Nl , YL , YU ,Nl , NO , NL) 

1 CONTINUE 

NF P (11* -l 
NFP ( 2 ) = 66 
NFP ( 3 ) = 50 


NF P ( 4 ) * 50 



NF P ( $ I = 686 
CALL FZPLHT19. 

, Nl 

, Nl, 

P 

t PC l » 2 


It L 

t “Nl , N 2 

1 

f N42 * i ,* 

' , 12 

, • PS I 

-OEGREES • 

t 

8 , 

ORD< n 

2 

, NFP , Nl , 

XL , 

XU ♦ 

Nl 

9 YL 

t 

YU 

, N 1, NO, Nil 


NFP (H* -2 
NF P ( 2 ) * 66 









~NF P ( 4 ) * 350 
NFP(5)» 380 



p 

, pi 1,6 





CALL EZPLOT (9. 

, HI" 

9 Nlf 


It L 

t “Nl 9 N2 

1 

, N42 , l ,• 

• , 12 

, ' 


« 

9 

8 9 

ORO ( 5) 

2 

f NFP t Nl , 

XL f 

XU 9 

NT 

9 YL 

9 

YU 

t N It NOt Nil 


NFP ( 2 1 = 50 









NFPlA) = 690 









NF P ( 5 ) * 40 






1 , L 



CATt EZPLOT (9. 

t Nl 

t Nl 9 


, >T1,7 


, -Nl , N2 

l 

• N*2 t l »* 

• , 12 

, * 


• 

f 

8 9 

0R0( 6) 


, NF P , Nl , 

XL 9 

XU t 

wr 

9 YL 

9 

YU* 

;NT, NO, Nl) 


NF P ( 1) =-l 









NF P( 2 ) = 50 
NF P( 3 ) *50 









NF P(4) *50 









NFP( 5) = 630 








1 

Call FZplotO. 
, N42 , l 

, Nl 
' , 12 

9 Nt 9 
, • PS I 

p , ? < 1 * 3 
-DEGREES' 

9 

)» L 
8 , 

T~-TT ~» N2 
ORDI 2 ) 

Z 

* NFP , Nl , 

XL , ‘ 

XJ 9 

Nl 

9 YL 

t 

~ YU 

• Nl, NO, Nl) 


NF P( l ) *-2 









NF P ( 2 ) * 66 









NF P ( 4 ) * 350 









NF P t 5 ) = 380 
CALL EZPLOT (9. 

, Nl 

9 Nl, 

P 

, pi 1 ,4 


It L 

, -Nl , N2 

T - 

t 9 l 

* f tz 

9 * 


9 

9 

IT 9 

ORDI 3 ) 



2 * NFP , N1 , XL , XU , Nl , YL , YIJ ,Nl, NO. Nl) 

KiF P 1 2 ) = 50 
NFP(<.)= 690 

NF P(5)- 40 

CALL EZPLOT (9. , Nl ' Nl. P , >11,5 ), L , -N 1 , N2 

1 . NA2 ,' * » 12 , • ’,8, QRO f 4) 

2 , NFP , Nl , XL , XU , Nl , YL , YU ,Nl, NO. Nl) 

RETURN 

END 
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SUPPflUTINF ST AC I M X* NY* MS TOP, MST ,0X1, RY , DRY , X , Y, U 6 , UC, V, IJSAV, SCAL S 

» S T A G 

1 


1 1 SF P ) 

STAG 

Z 

C PR r G P A ^ CALCULATING THE PCUNCA»Y LAYER PROFILE NEAR 

STAG 

3 

C 1HE STACK* TION POINT 

STAG 

4 

c 

COMMON /HI 1 / NT IMF, NQIMC , ISTO 

STAG 

5 


0 IMF* SI ON US A V (300, 100) ,SCALS!300) 

STAG 

6 


DIMENSION PHI 7(24) , PH l P ( 24 ) , FTA P( ?4 ) 

STAG 

7 


DIMENSION X ( 300 ) *Y(100) *UE(300,3),UC( 100,3J,V( 100,2) 

STAG 

9 


DIMENSION FF(IOO) tEFPClOO) 

STAG 

9 


HAT A ETAP / 0 • t.2 v«4f«6t *89 l. f 1, 2 r l.4» I«6fl.8f2.»2.2f?.4f2.6v?.8f3 

. STAG 

10 


1 >3.2 t3 • 4 * 3 • 6 1 3 • B f 4 » t 4*2 t 4 • 4i 4 * 6 / 

STAG 

1 i 


DATA PHIZ /O. ,.0233 ,.0381, .1867, . 3124, .4592, .622, .7967, . 079 3, 1. 1 63ST AG 

12 


19, L. 36 2, l. 5578, 1.75 53,1.95 38 ,2. 15 3,2. 3526, 2. 552 3,2. 75 22 , 2.9 5 21, 3. 

1STAG 

i 3 


1521,3.3521 ,3.5521,3 .7521 ,3.9521/ 

STAG 

14 


DATA PHI P /0.,. 2 266*. 41 4 5,. 5663, .6859, .7779, .8467, .8968, .932 3 , .9 56 STAS 

15 


18,. 9732 , .9839, .9995 f .9946 ,.997, .9984, .999 2 , . 999 6 , . 999 8 , .999 9. 1. , 1 

. STAG 

16 


itl., 1 ./ 

STAG 

17 


BAG =. 08 

STAG 

1 8 


IF(ISFP) 10,10,5 

STAG 

19 

5 

B AG = • 5 

STAG 

20 

10 

EFtl) = 0 . 

STAG 

21 


FFPT 1 T = o. ' “ ’ ' 

STAG 

ZZ 


DO 20 M= 1 , MX 

STAG 

23 


IFIUFIM.U) 20,20,19 

STAG 

2 4 

IS 

MSP = M 

STAG 

25 


GC TO 21 

STAG 

26 

20 

CONTINUE 

STAG 

27 

21 

A STAG = (UF (MSP+2,1 )-UE fMSP+l ,1) )/(X( MS P + 2 )-X(MSP+l) ) 

STAG 

28 


IF(ASTAG) 22,22,23 

STAG 

29 

22 

ASTAG=IUEI MSP, i ) -UE OSP-1,1) ) / 1 X ( MSP ) -XI MS P- 1 ) ) 

STAG 

30 

23 

SQAS = SORT(ASTAG) 

STAG 

31 


OflT 2.6 7 SO AS 

STAG 

32 

3 C 9 

tF(0ELT-YIN Y-3 ) ) 311,310,310 

STAG 

33 

310 

R Y=R Y+DRY 

STAG 

34 


CALL YSFTIRY,Y<2) ,NY,Y) 

STAG 

35 


go rno9 

STAG 

36 

311 

CONTINUE 

STAG 

37 


DC 80 N=2fNY~ " 

STAG 

38 


YET = Y ( N) *SQAS 

STAG 

39 


Da"T3"W^T7?^ 

STAG 

40 


IF(YET-ETAPINN) ) 408*408,33 

STAG 

41 

4C 2 

mark = NN ~ 

STAG 

42 


GO TO 410 

STAG 

43 

33 

CONTINUE' ‘ 

STAG 

44 


FF(N> = YET-. 6479 

STAG 

45 


FTP(K)' = 1. ' ‘ " 

STAG 

46 


GO Tr 80 

STAG 

47 

410 

FRACT = I YET-FTAP (MARK-1) )/( ET API M ARK )-ET API MARK- l ) ) 

STAG 

48 


FRAC. 1 = 1. -FRACT 

STAG 

49 


EFIN) = PHTZTM7SRK-I1*FRAC1+PHTZ (MARK) *FR ACT 

STAG 

50 


EFP(N) = PHIP(MARK-l)*FRACl+PHIP(MARK)*FRACT 

STAG 

51 

8 C 

CONTINUE ~ “ 

STAG 

52 


Ml = MSP-MST9P 

STAG 

53 


M2 * MSP+MSTDP' 

STAG 

54 


79 




STAG 

55 

50 

f M=M+l 

STAS 

56 


MST-Hfl 

STAS 

57 


SCALS(M) =RY 

STAS 

58 


00 71 N*l,NY 

STAG 

59 


UCIN,3) * UC ( N * 2 ) 

STAG 

60 


UCIN.2) » UE(M,1}*EFP(N) 

STAG 

61 


V (Nt 2) = V { N »l ) 

STAG 

62 


V(N,1) * -SQAS*EF(N» 

If II STD . EQ. U GO TO 71 

STAG 

63 


USAV(M,N)=UC(N,2) 

$TAS 

64 

71 

CONTINUE 

STAS 

65 


IE (M-H2) 50,55,55 

STAG 

66 

55 

IF(UFIM,1) —BAG ) 50,50,81 

STAG 

67 

81 

CONTINUE 

STAG 

68 


RETURN 

STAG 

69 


END 

STAG 

70 


80 


25 

50 


75 


SLRROUTINF A ft PR (PRFC ,XS l G, NS IG, A Si?, AS * AR» CM At , R M AT , MGA V , nj f 

» ACAP t TATTPR 

l 

lMlCK,RDBB,GANAW,UINF t l 001 , 0 X 1 ,Br.AP) 

ATTPk 

2 

nil-ENSl 3N XSIG I 100) ♦ ASZ (30 1, AS t 30,30) » ARI 30) , BC API 10 0, 3 ) 

ATTPR 

3 

0 1 MF NSI DN AC AP ( 30 ,3 ) «THICK(24)« GAM AM ( 1000) 

A TTPR 

4 

nCURIF PRECISION CM AT (60*60 ),RMATI130) 

ATTPP 

5 

PI*3.14150 

ATTPR 

6 

NGPl = NGA M+ 1 

ATTPR 

*T 

00 50 M =1 f NGPl 

ATTPk 

ft 

CMAT(M, 1 )«ASZ(M) 

ATTPR 

9 

RPAT t M) = AR 

ATTPR 

10 

DC 25 N=i 1 NGAM 

ATTPR 

1 1 

(!RAT(M,N+1)=AS(M,N) 

ATTPR 

12 

CONTINUE 

ATTPR 

13 

CALL ALSOL ( NGPi ,CMAT , RMAT) 

ATTPt 

14 

00 75 M=l,NGPl 

ATTPR 

15 

ACAP(M,1) = RMAT (M) 

ATTPR 

16 

G AMA K( i 1 sGAMi t AC AP , OX I , PIT 

ATTPR 

17 

SAVE^XSIGC NSIG+l ) 

ATTPR 

1 ft 

XSIG(N$IG«*1>=2. 

ATTPR 

19 

CALL CPC (Ot NGAM,NPf XSIGtNS IG« XSIGt NS IG, XS IG, NS I 3, AC AP , B CAP , 

MICK tRATTPR 

20 

l OBB , GA PA W, UI NF , UHHT , 1* , S AVE, DXI , PREC ) 

A T T 3 P 

21 

XS!G(NSIG*1)*$AVE 

ATTPR 

2 ? 

RE TO** 

ATTPR 

«_ _s 

END 

ATTPR 

2 4 


81 


SUBRDUTI NE UNPOP ( NG AM f AR , ALAM , A FACT, R M AT , 
IP, XSIG,8CAP,THICK,RORB,UINF,XC,UF) 

D I MENS! ON AR(30) , ALA M( 3 0 ) , XGAMI 30 ) , AS ( 30, 
l RCA P (100 , 3) t THICK ( 2 A U X C ( 3 00 ) , J E ( 3 00 , 3 ) 
DOUBLE PRECISION R M ATI 130 V , C M AT { 6 0, 60 ) 

NGPl «NG AMf 1 
DO 5 M*i,NGPl 

SU8=AR ( M) - ALAM ( M ) *AF ACT / 3 . 

R PAT I M) * SUB 
CMAT<M,1)*1.__ 

C MA T ( M, 2 ) = XG AM ( M j 
DO 5 N*2 , NG AM 
C MA T ( M, N+l ) ® AS ( M ,N) 

CALL ALSOU NG_P1 , CM AT , RMAT ) 

DO 10 N= 1 , NG P 1 
ACA P(N,1)*RMAT(N) _ 

DO is M*1 , MX 

SI GN*1 » 

IF(M-NZ) 12,14,14 
SIGN=- SIG N 

CALL QECAL t 0 » NGAM, NGAM, NF,X$ IG, ACAP,BCAP, 

1 UE t M , 1 ) tSIGNl 

"CONTINUE 

RETURN 

END 


CMAT,XGAP, AS, A CAP, MX,NZ,NUNPOP l 

UNPOP 2 

30), AC API 30, 3) , XSIG( I JO), UNPOP 3 

UNPOP 4 
UNPOP 5 
UNPOP 6 
UNPOP 7 
UN POP 8 
UNPOP 9 
UNPOP 10 
UNPOP 11 
UNPOP 12 
UNPOP 13 
UNPOP 14 
UNPOP IS 
UNPOP 16 
UNPOP 17 
UNPOP 18 
UNPOP 19 
UNPOP 20 

THICK, RDBB,0., U INF, XC< Ml, UNPOP 21 

UNPOP 22 
UNPOP 23 
UNPOP 24 
UNPOP 25 
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SUBROUTINE ALSOL (NT , C, R ) 

OCUB LE PRECISION C NDIMC.NOI MC I , R ( 130 > ALSO!. 

DOUBLE PRECISION CM AX f SAVE, SUM 
COMMON /RL1 / NTI ME_, NDl_MC 

NT1 = N'T-1 * ’ 

DO 99 J=l,NTl 
C MAX = C(NT,JJ 
L*NT 

DC 1C I *J,NT1 

IF (DARS(CMAX)-DABS (CII ,JI I) 5,10,10 
5 C MAX « (T(r,j) 

L = I 

10 CONTINUE 

DC 15 JJ=J,NT 
SAVE = CIL.JJI 
C (L, JJ) * C ( J,JJ) 

15 C l J , JJ) = ‘SAVE/CMAX 
SAVE = R(L> 

R(L) = R(J) 

R ( J) = SAVE/CMAX 
JP1 * jVi 
DO 25 I=JP1 , NT 

DO 70 J J=JP1 ,NT 

20 C(I,JJ! = C(I,JJJ - C ( I , J ) *C ( J , J J ) 

2 5 R (I J = RTH - ”RTJ ) * CT I7JT 

99 CONTINUE 

R I NT ) = R ( NT >7 C I NT , NT )' 

DO 150 K = 1 » NTI 

l *NT -k 

IPl * I+l 

sun - o; ' . 

DO 125 J*IP1,NT 

125 $ UN ~ SUTT ftrJT*C(T7Jr " — ' 

150 R(I| * RII) - SUM 

RETURN 

END 
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SUBPnuri NF CPCIISF.P ,NGAM f NF f XSIG, NSIG,XS IGA, NS I G A , X S IG8 ,N S I G3 , ACAPCPC 


l,BCAP f THtCK,R0B8,GAMAW,UINF,UQQT, S IGN, XC, OX I, CP ) CPC 

DIMENSION XSIG(IQO) , XS l GA UOO J , X S IG8< 100), AC API 30, 3 ) , BC AP { 1 00 , 3 ) CPC 
0 1 ME NS I O N GAM A W(IO OO), TH I CK ( 2 41 CPC 

THFT'a=4RCTIXC1 ~ CPC 

R EC I P- 1 • / I U! NF*U I NF I CPC 

SUM=0. CPC 

ANGLE -0. CPC 

DC 5 N*i , NF CPC 

ANGLE S A NG L E «• TH E T A CPC 

5 SUM- SUM* TH I CKf N) *COS f ANGLE ) CPC 

CP-UDPT*RFClP*(THICKIi) +2.*( l .-XC) *SUM) CPC 

CALL OECALU SEP,NGAM,NS IG,NF, XS IG, ACAP, BCAP, THICK, RCBB, GAMA d { U , UICPC 
INF, XC,U, SIGN) CPC 

CP<P*2. *< SIGN*U/UINF-1. ) CPC 

CAL Li GA JU CitNGAM, ACAP, BCAPI 1,1 ),XS IG 11) , XS I G( N S IG» 1), 3 AMAWI l), XC,CPC 
1VAL1 I CPC 

CALL FGAMI ( 2 ,NG AM , AC AP, BC API 1 , 2 ) , X S IGAC1) , XS IGA ( NS IGA + 1 ) , G AMAWI 2) ,CPC 
l XC , VAL2 ) CPC 

CALL EGAMI (3,NGA M,ACAP,BCAP( 1 ,3 ) , XS IG B( l ) , XS IGBINSIGB+l ),GAMAWf 3) ,CPC 
i XC , VAL3) CPC 

JC P*CP»SIGN»RECIP*IU5*VALl-2.»\f AL 2 ♦ *5 AL 3 ) / OX I CPC 

it (I SEP) 2 0,20,10 ~ CPC 

10 CALL FSIGI CitNSIGf X SIGt BCAP, XCtV ALII CPC 

CALL ESIGl (2 tNSIGA, XSfGA, BCAP,XC , VAL2 ) CPC 

CALL E SIGI (3,NSIGB y XS IGB , BCAP,X C, V AL 3 ) CPC 

CP*CP*RECIP*I1.5*‘VAL1-2.<‘VAL2 + .5*VAL3I/DXI CPC 

2 0 CP*-CP CPC 

RETURN ” " CPC 

ENO 


1 

2 

3 

4 

5 

6 
7 
B 
9 

I 0 

II 
12 

13 

14 

15 

16 
1 7 
18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


r. 
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SUSP OUT INF CLC M1NC01 ,ISFP,NGAM f XSIG,NSIG, XSIGA, MS T GA f X S IGB , MS IGR f ACL C* 
ICAP,BCAP,THICK,RDBB,C< MAW , U I NF, UOOT , OX I , AROT f °A ) 

COMMON /CLCMBL/ CLVB, CPVB, CMPAVB 
01 MF NSl ON ARGL(2l) , ARGMI21) 

01 ME NS I 3 N G AM A W ( 100 0 ) * t H I C K ( 24 ) 

0 IMF NS I ON XSIGC100) , XSI GA ( 100 ) , XS IGR( 100)»ACAP( 30, 3 NBC API 103,3) 

FORMAT I / /40X ,4HCL * E 13. 5/40X , 4HC M =F13 .5 , 17H (ABOUT M IGCHOR )) /40X, CLCM 
14HCM =F13.5,24H (ABOUT PITCH AXIS - A =1=7.4, 1H) ) 

M0LT=6 

SAVe = THtCK(U 
THICK(I) =0. 

0 T=3. 141 59/FLOAT (NCOI) 

C L =0. 

CM = 0. 

Xt=-i. 

ANGLF =0. 

FLI=0. 

F M I =0. 

IF(I SFP) 5,5,7 
XATT=XSIG[NSIG+1) 

IF ( X ATT- • 95 ) 8,5,5 
XAQ=XATT+5.F-4 

XA1^VxXQr.“0T5 ” 

C1=-.5*(K + XATT) 

C2=C 1+ XAtT 
C 1P=.5*(1.-XAP> 

C2P=CI PVXAP 
DC 10 I =1 , NCOI 
A NG L F = A NGL f ♦ 0 T 
XIP1=C1*C0S (ANGLF) + C2 

C All C PCrrSFPTNG A M, 1 ,XSTG , NS t G, X S I G A, NS I G A , X S I GB , N S I GB , AC AP , BC 
lHICK 5 PD3B ? GAMAW.niNF^U0nT ,1.0 ,X I PI, DX I, CPU I 
C ALL "C PC rrSFP,“NGA3T,I ,X$ ! G, NS 1 G, XS IGA* NS IGA,XS IG B,NS IGB, A CAP , BC 
1HICK,RDBB,GAMAW,UINF,UD0T,-1. , X I P l , OX I , CPL ) 

FLIP 1=C PL -CPU 
FMIP1 = XI P1*F LI PI 
CLiCL+TXm-XTRTFLT PUFLT) 

CM=CMM XIP1-XI ) *<FMIPH-FMI> 

X l =X T PI 
FLI-FLI PI 

FMI=FMI Pi 

XI*1. 

flt=o. 

FMI=0. 


nc 15 I =1 , NCOI 

ANGLF=ANGtF^DT 

XIP 1=C1P*C0S( ANGLE) +C2P 


lHlCK f ROaBtGAMAW,UlNFfUOOT, 1.0, XI PI, OX I, CPU) 

lHICK,ROBB,GAMAW,UINF!uOOT|-l.,XIPi,OXI,CPL) 
F LI P1=C P'L-CPU 


FMIPl = Xl P1*FLIP1 

rL=XT-T XTFi-XTT* (FLI Pl + FL I r 


a r. l c y 

L 

Cl CM 
MAIN 

2 

CLCM 

3 

C LC M 

4 

CLCM 

5 

* C L C M 

6 

CLCM 

7 

CLCM 

B 

CLCM 

9 

CLCM 

10 

CLCM 

l 1 

CLCM 

12 

CLCM 

13 

CLCM 

14 

CLCM 

15 

CLCM 

15 

CLCM 

17 

CLCM 

IS 

CLCM 

19 

CLCM 

20 

Cl CM 

21 

CLCM 

c — 

CLCM 

23 

CLCM 

25 

CLCM 

25 

CLCM 

26 

CLCM 

27 

CLCM 

28 

CLCM 

29 

TCLCM 

30 

CLCM 

31 

TCLCM 

32 

CLCM 

33 

CLCM 

34 

CLCM 

35 

CLCM 

36 

CLCM 

37 

CLCM 

38 

CLCM 

39 

CLCM 

40 

CLCM 

41 

CLCM 

42 

CLCM 

43 

CLCM 

44 

CLCM 

45 

CLCM 

46 

CLCM 

47 

.TCLCM 

48 

CLCM 

49 

.TCLCM 

50 

CLCM 

51 

CLCM 

52 

CLCM 

53 

CLCM 

54 
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CM=CM-( XIPl-XI)*(FMIPl+FPI) 

CLCM 

55 


XI=XIPl 

CLCM 

56 


FLI=FHPi 

CLCM 

57 

15 

fmi=fmi Pi 

CLCM 

53 


XI PI =XAQ 

CLCM 

59 


on 16 1 * 1,21 

CLCM 

60 


CALL CPS (I SFP.NGAM, 1 ,XS I G, NS I G, XS I GA, NS I G A, X S I GB, NS 1GB, ACAP 

, 3CAP , TCLCM 

61 


1HICK ,Rf)BB,GAMAW,UINF ,UOOT ,l.O,X I Pl.DX I .CPU) 

CLCM 

62 


CALL CPC (I SEP.NGAM.l ,XS tG.NSIG.XSIGA, NSIGA,XSIGB,NSl GB. ACAP 

,BCAP, TCLCM 

63 


1HICK.RDBB ,G A MAW .UINF.UDOT.-l. .XIPl.OX I.CPL) 

CLCM 

64 


ARGLU )=CPl-CPU 

CLCM 

65 


ARGM ( I ) = XI Pl*ARGL( I ) 

CLCM 

66 

16 

XIPl =XI PI ♦ *00125 

CLCM 

67 


SUML=0. 

CLCM 

68 


SUMM=0. 

CLCM 

69 


00 17 l =1,19,2 

CLCM 

70 


SUML*SUML + 2.*ARGLU )+4.*ARGL( 1+1) 

CLCM 

71 

17 

SUMM= SUMM + 2 • * ARG MCI 1+4. *ARGMt !♦!) 

CLCM 

72 


CL=CU-0.833333E-3*(SUMH-ARGL(2l)-ARGt (1)1 

CLCM 

73 


C M=C M+ 0. 833333 F-3*(SUMM+ A RGP< 21 )- ARGM ( 1 ) ) 

CLCM 

74 


BC0N=16.*BCAP(1 ,l)*SQRT <5. £-<♦*< XATT-XS IG( 1) ) l/UINF 

CLCM 

75 


C L*C L+BC ON 

CLCM 

76 


CM<M+XATT*8C0N 

CLCM 

77 


GO TO 100 

CLCM 

78 

5 

00 99 I =1 , NCOI 

CLCM 

79 


ANGLE=ANGL E+OT 

CLCM 

80 


XIPl^-COSl ANGLE) 

CLCM 

81 


CALL CPC (I SEP, NG AM, 1 ,XS I G , NS I G, XS I GA, NS IGA, XS IGB, N S f GB, ACAP 

tBCAP, TCLCM 

82 


1HICK ,RDBB , GAMAW , UINF ,UDOT ,1.5 ,X l P l ,DX I , CPU ) 

CLCM 

33 


CALL CPC ( I SEP, NG AM, l , XS l G , N$ I G, XS I GA, NS IGA, XS IGB,N S IGB, ACAP 

,BCAP, TCLCM 

84 


1HICK , ROBB , GAMAW ,UI NF ,U0 OT , -l • ,X I P 1 , DX I , CPL 1 

CLCM 

85 


FLIP1=CPL— CPU 

CLCM 

86 


F Ml P 1= XI P1*FL1P1 

CLCM 

87 


CL=C L+ ( XI PI— XI 1 * <FL I Pl+FLI) 

CLCM 

88 


CM=CM+<XIPl-XI)*(FMIPl+FMn 

CLCM 

89 


X! =X!P1 

CLCM 

90 


F L I * F L I P l 

CLCM 

91 

99 

FM1=FMI PI 

CLCM 

92 

100 

CL*. 25*CL 

CLCM 

93 


CM*-.125*CM 

CLCM 

94 


cmpa»cm<-arot*Cl*.5 

CLCM 

95 


WRI TE ( MGUT ,A ) CL .CM.CMPA , AROT 

CLCM 

96 


THIC K( 1 ) *S A VE 

CLCM 

97 


CL VB * CL 

MAIN 



CMVB = CM 

MAIN 



CMPAVB * C MPA 

MAIN 



RE TURN ■*" 

CLCM 

98 


END 
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SU3RCUTI NE QFCAL (ISEP,NGAM,NS 1G, NF ,X S IG, AC AP , RC AP, T HIC< ,RDH3 ,GA 1MAQFCA L 


1 ,UINF,XC ,U,SIGNI OFCAL 

DIMFASION ACAP(30,3) , RC AP 1 100 , 3 ) , XS l G (100 ) OFCAL 

DIMENSION THJ^CK (24) OFCAL 

F PS = 1 . F - 6 OFCAL 

CORR=. 7071 07/(1.-. 63662 *S CRT < RDBB ) +. 2 5*R CBB ) OFCAL 

SI NT* SORT! 1 . -XC*XC) OFCAL 

THFTA=ARCT 1XC) OECAL 

C CHIN T=0. OFCAL 

SUM = 0. _ OFCAL 

SlNT2*St N( ,5*THETA) OECAL 

CCST2=C DS( .5*THFTA) OECAL 

IF (SINT -EPS) A, 6, 6 OECAL 

A F AC T = THF TA* . 5 OECAL 

GO TP 8 OECAL 


6 FACT*(1.-XC) /SINT OECAL 

8 DO 10 N=1,NF OECAL 

COUNT=COUNm. OECAL 

A NG L F = TH E T A *C OUN T OECAL 

SUM* SUM+TH ICK(N)*(COUNT*FACT*STN( ANGL E )-CDS ( ANGL E ) ) OFCAL 

10 CONTINUE OECAL 

U=2.*SIGN*UINF*COST2*SUM+ACAP<l,ll*SINT2*.2 5*C3ST2*( l. + XCIM 3.*XC-QECAL 

11. )*GAHH"A QECAL 

SUM*0. QECAL 

ANGLE=0. ‘ QECAL 

DO 12 N= l * NG AM OECAL 

ANGLF=ANGLE+THETA ' ' QECAL 

12 SUM=SUM*ACAP(N«-l , l ) *S IN ( ANGLE ) QECAL 

(J=U«-C0ST2*SUM " ' OFCAL 


25 


40 “ 

46 


IFUSFP) 25,99,25 
SUM=0. 

XSFP* XSI G ( 1 ) 


X ATT=X SIGTNSTGFl) 

DD 40 N=2,NSIG 

~ SUM=SUM*BC AP(N,1 l*FB<X$tG(N-l),XSIG(NT, XS IG( N+l ) , XC. I 
IF(XC— XA TT— EPS) 45^45^46 

IXATT-4.*XC)-SIGN*( l!-SQRT( (XSEP-XC)/( XATT-XCM) 


OECAL 
QECAL 
OECAL 
OECAL 
OECAL 
QECAL 
QECAL 

XC) /(XC-X4TT >)*( 3.*QFC4 L 

QECAL 


GO TC 55 T " “ QECAL 

45 IF(XSEP-XC) 49,49,48 QECAL 

“ 4T F AC T=- 51 GN* ( l .—SORT ( IXS EP-XC I / (X ATT-XOTl QFCAL 

GO TO 55 QECAL 

49 FACT*- SIGN QECAL 

55 U=U*C0ST2*(BCAP(1 ,1 l*FACT+SIGN*SJM) QECAL 

99 Tf= I S IGN*UI NF +SQRT ( 1 . *XC ) ♦ ^^R1^JT7SQRTn:.'+XC*.5*RCBR) QECAL 

RETURN QECAL 

“ FND 


1 

? 

3 

4 

5 

6 

7 

8 
n 

1 0 
1 1 
1 2 

13 

14 

15 
L 6 
1 7 
18 

19 

20 
21 
22 
I j 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
3ft 

39 

40 

41 

42 

43 

44 

45 

46 
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r> 3* > 


SLRRPUTINF YVRtY* !) Y VR 

REAL Y<10) YVR 

REAL MVR YVfl 

CCMMQN /INPTVR/ FTVft(64), FPVB(64)t FPPRVRl64)t PI OR V8 ( 64) ♦ YVB 

XMVB<64)* OELVB t XMUVB* FOVR, XMUAVR» YVB 

ATOVB t ATCVB » ATSVB, ROVB, RVB(64)t YV3 

MVB ( 64 ) f NVB YVB 

Ytn = ( R VB ( I ) - OELVB ) **2 * MVBUI YVB 

Y(2) = FPVBm**2 * HVB(I) YVB 

Y(3) - FTV8n)**2 * D I QRVBC I ) YVB 

YU) = (nFLVB - PVBU)I ♦ FTVBU) ♦XMVB(l) * MVRC I I YVB 

Y (5) * FPVB(I) * FTVSm * XMVB(l) * MVBCi) YVB 

Y { 6 ) = RVB(I) * (DELVB - RVB ( I ) ) * MVB(I) YVB 

Y ( 8) = (RVBU) - OELVB) * FPPRVR(I) * FTV R( I ) * XMVBU > * 1 VB< l ) YVB 

IPl ■ IM YVB 

1 F < IP l .GF . NVB) GO TO 1 2 YVB 

SUM * 0. YVB 

DO 10 J » IPLi^NVB YVB 

10 SUM * SUM — t RVB t4*i ) - RVB(4)> * IRVBI4 + 1J * MVBU + ll YVB 

A _♦ RV BCJ) * MVR(J) ) YVB 

12 Y<7> '« FPPRVBU) ** 2 * SUM / 2. YVB 

RETURN YVB 

END ' 


1 

2 

3 

4 

5 

6 
? 
8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 
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SURRCUTI NP POLLY IN, BBS , R EL , AN , AA ) 

IMPLICIT RF AL *8 <A-H,0-Z) 

r complex roots n c a polynomial bairstows method 

DIMENSION A [3 3) , ANf 6 Ol , _C ( 26 ) , AR A R< 26 ) , B( 30 ) , AA( 30 
I 1 1 = 1 
7 NP1 = N* 1 
NPPl=Nf 2 
DO 6 Cl 1=1, NPI 
LLL=NPPI-I 
601 A ( I ) =A A ( LLL) 

13 nn 14 K = l , NP1 

14 A BAR ( K> = A ( K ) 

AB SSQ = BB S*BBS 
REL SQ=RE L*REL 
NB A R=N 

8 1 1 ) =A { 1 ) 

cm=A(U 

15 IFIN8AR-2) 200,210,17 

17 PI =• 2 
Ql = . 1 

18 I TF R = 0 

19 PI =P1*5 • 

QT =m*T(T. 

33 P=P1 

Q=Ql 

NRPi =NB AR+1 

34 1=1" 

LAST = NB AR 
DTF ST=9. 99036 

C BA1 RSTOW ITERATION 

37 Bm**BAR(2»-P*Bm 

DO 40 K=3t NBPl 

40 "B (K) =ABAR< K)-P*B (K-"n-Q*BTK-2J 
45 C (2) =B(2)-P*Cin 

Do 501 IC»3*rRT ‘ 

50 C<K)=Bm-P*C(K-l)-Q*C(K-2) 

CTLA ST)=C(LASTI-B(L ASTI 

D=C ( LAST-l ) «C ( LAST- 1 )-C< L AST ) *C< LAST-21 

D'SOR -D*D 

I F (DSQR-1# D— 36 ) 19,1 9,60 

67TDELP = lH { LAS T) *0 (LAST-1) -Bl IAS 1+11*01 L AST-2)) HT 
DE LQ = ( B ( LAST<*1 ) *C ( L AST- 1 ) -B( L AS T )*C( L AST ) I /D 
C TEST FOR CONVERGENCE 


RELP=DELP/P 


RELQ=DELQ/Q 


RELPS=R6LP*RELP 


KbLUS=KbLH*K¥lU 


DE LSQ=RE LPS+RELQS 


P =F*DE L P 


0=0+0E LtJ 



1 F (RFLPS— RELSU) 70,70,65 
65 IF(DELP*DELP-ABSSQ> 70,70,80 

7 0 I F f RFLOS— HFLS ' O) 1 20,120, 7 5 

75 IF(DELQ*DELQ-ABSSQ) 120,120,80 
80 GET TC T9 0,100) ,L 


POLLY 1 
POLLY 2 
POLLY 3 
POLLY 4 
POLLY 5 
POLLY 6 
POLLY 7 
POLLY 3 
POLLY 9 
POLLY 10 
POLLY 11 
POLLY 12 
POLLY 13 
POLL Y 14 
POLLY 15 
POLLY 16 
POLLY 1 7 
POLLY 18 
POLLY 19 
POLLY 20 
POLLY 21 
POLLY 22 
POLLY 
POLLY 24 
POLLY ? 5 
POLLY 26 
POLLY 2 7 
POLLY 2 8 
POLLY 29 
POLLY 30 
POLLY 31 
POLLY 32 
POLLY 33 
POLLY 34 
POLLY 35 
POLLY 36 
POLLY 37 
POLLY 3 8 
POLLY 3 9 
POLLY 40 
POLLY 41 
POLLY 42 
POLLY 43 
POLLY 44 
POLLY 45 
POLLY 46 
POLLY 47 
POLLY 48 
POLLY 49 
POLLY 50 
POLLY 51 
POLLY 52 
POLLY 53 
POLLY 54 
POLLY 55 
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SO ITFR=ITERU 

IF ( 250- ITER)310,37,37 
ICO IF ( DTE ST-OELSQ) 34,34,110 

110 DTEST=DELSQ 

0 (2) =A(2 )— P*B(l ) 

DO 115 K*3,NPl 

115 BIK) =A(KI-P*B(K-l)-0*9<K-2) 
GO TO 45 

C ITERATION HAS CONVERGED 

120 GO TO 1133, 143), l_ 

130 L =2 


LAST =N 
GO TO 110 

C FACTOR OUT QUADRATIC 

140 N8AR*NBAR-2 
NBP1=NBAR+1 

A BAR 12) 3 AB ARC 2 1 — P*ABAR I 1 ) 

DO 150 K*3,NBP1 

150 ABAR (K) =>ABAR(K|-P*ABAR( K-l I— QAABARCK-2 I 
GO TO 250 

C SCLVF LINEAR EQUATION 

200 NBAR=NBAR-l 

Rl*-ABAR (2 ) /ABAR (1 ) 

R2=0. 

GO 70262 

C NORMALIZE QUADRATIC 

210 P»AB AR ( 2 ) /ABAR I 1 ) 

Q=ATAR(3)/ABAR(l ) 

NBAR»NBAR-2 — 


C SOLVE NORMALIZED QUADRATIC 

250 Rl=-P/2. 

C1*RI*R1-Q 
IF (Cll 270,280,260 
260 C 1*0 SORT (Cl) 

R2*R 1-C1 ' "" 

R1*R1+Cl 
262 Cl*0.' 

GO TO 290 

270 Cl=-Ci ~ 

C1*D SORT CC1) 

2 90 _ r2»R’i - 

290 C2=-C1 

AN! 1 1 1 ) *C1 ~ 

ANUini)*Rl 

SNfm*2)«C2 " ' - 

ANI 1 11*3 ) *R2 

111*1 IT* 4 

I F ( NBAR-1) 4 ,200 ,15 

C SPEC ISl CONDITIONS' ~ 

310 WRITE (6,600) 

600 FORWATTI X,'50Hlira CCNVFRCENCn N 25ff ITERATIONS .POLLY HAS SPOKEN) 
A CONTINUE 

RETURN” 


END 


POLLY 

56 

POLLY 

57 

POLLY 

58 

POLLY 

59 

POLLY 

60 

POLLY 

61 

POLLY 

62 

POLLY 

63 

POLLY 

64 

POLLY 

65 

POLLY 

66 

POLLY 

67 

POLLY 

68 

POLLY 

69 

POLLY 

70 

POLLY 

71 

POLLY 

72 

POLLY 

73 

POLLY 

74 

POLLY 

75 

POLLY 

76 

POLLY 

77 

POLLY 

78 

POLLY 

79 

POLLY 

80 

POLLY 

81 

POLLY 

82 

POLLY 

83 

POLLY 

84 

POLLY 

85 

POLLY 

86 

POLLY 

87 

POLLY 

88 

POLLY 

89 

POLLY 

90 

POLLY 

91 

POLLY 

92 

POLLY 

93 

POLLY 

94 

POLLY 

95 

POLLY 

96 

POLLY 

97 

POLLY 

98 

POLLY 

99 

POLL Y100 
POLL Y101 
POLL Yltf2 
POLL Y103 
POLL Y104 
POLL Y105 
POLL Y106 
POLL Y10T 

TOLIYIO0 
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o o o a 


41 

l 0 
20 


50 


60 


53 


32 

33 
31 


SUP R CUT l NT $FfTT(LGQ,M,NV,RF8 ,X,Y,UC, PRESS , GRAD , DEL T, 01 SP# THE 
1 SC ,MTRAN) 

SUBROUTINE FOR CALCULATICN OF 80UNCARY LAYER TH!C<NE$S, 
DISPLACEMENT THiCKNIFSSt MOMENTUM THICKNESS AND EDDY V I SC n S l TY • 


DIMFNSI IN X 13001 ,Y( 100 ) t UC ( 100 , 3 I ,V!S C< 100 * 2 I « GRAOC 1001 
RTR = SORT (REB) 

NY = NV ♦ 2 

UEOGE = . 995 *UC(NJY t 1 J 

DO 10 N=l, NV 

IF (UEDGE-UC ( N+l , 1 ) ) 41 , 41,10 
NOELT = N 
GO TP 20 
CONTINUE 

DELT = YC NOELT) 4* (UEOGE-UC (NOELT , 1 > ) ♦ t Y ( NOELT +11 -Y( NOFLT I ) / I JC 
1T> 1 , 1 ) -UCT NOEL T ,U) 

SUM = 0 . 

DO 50 N- 2 ,NY 

SUM = SUM* ( Y INJ^-Y ( N-l) )_* ( UCC N, 1 1 *UC{ N— It II) 

Dl SP = ( YIN Y ) - . 5 *SU M/UC ( NY , 1 H / RT R 
SUM = 0. 

uFo of =~ uc" r ny a r i 

DO 60 N «2 , NY 

SUM = SUMf TY(N) -YIN-in *rCUFOOE-UC IN , 1 ) ) *UCC N, 1 )MUEDGE-UC( N- 
1 *UC ( N-l , 1 ) ) 

theta = . 5 * sum/ c rtr*u~edg~e **2 ) 


TA , VI SETUP 
SETUP 


1 

2 

3 

4 

5 

6 
7 
9 
9 


IFCLGOI 53,53,56 


NVFl=NV*I 
EASE * l. 

IF ( M-MTRAN7 


1T73273T 


IFCMTRAN+5-M) 31,31 ,33 

EASE = ( Xt M)-X< MTRANn/TXTRmmTT-XT'MTR'AN ) ) 
CONTINUE 

INNER =0 * " 

F AC 1 * • 16 *RTR*E ASF 

TGI 6 8 * U FOG E *T5T SF^REWFASTE 


FACZ~= 

FFAC1 

FFAC2 

TAUW 


-RTR/26. 


= GRAD ( 1 ) /RTR 
N*2"VN'VPT 


402 
4 01^ 

7CL 

702 


"DO" T57T 

ALTER = U*FAC2/U.+5.5MYINI/DELT)**6I 

TFT I NNER) 402 ,"4"01I4 _ 02 

VISC (N,i)=ALTFR 

GTT T O T6'0 

CONTINUE 

T A UM Y = T AUW-YT Nil * FF2TC2 
IF(TAUMY) 701,701,702 


nr rr 

521 


VI SC IN, 11*1. 

GO TP 703 

FK* YIN T*FF AC I ♦ f S"QRT (TSUM Yl' 

VI SC ( N, l ) = 1*+FAC1*Y(N)*YCN)*A6SIGRAD(N) l*( U- EXPC EX )1 «2 
IF (VISC rST,T I— ALTFRJ" 

VI SC (N, 1 1-ALTER 
IRNFF^T 


160 ,160, 3“ZI " 


SETUP 
St TUP 

sr tup 

SETUP 
SETUP 
SETUP 
SETUP 
SETUP 10 
SETUP ll 
SETUP 12 
SETUP 13 
SETUP 14 
SETUP 15 
(NDELSE TUP 16 
SETUP 17 
SETUP 18 
SETUP L 9 
SETUP 20 
SETUP 21 
SETUP 22 
SETUP !j 
SETUP 24 
l, 1)1 SETUP ZZ 
SETUP 26 
SETUP 2 7 
SF TUP 2 R 
SETUP 29 
SETUP 30 
SETUP 31 
SETUP 32 
SETUP 33 
SETUP 34 
SETUP 35 
SETUP 36 
SETUP 37 
SETUP 38 
SETUP 39 
SETUP 4 0 
SETUP 41 
SETUP 42 
SETUP 43 
SETUP 44 
SETUP 45 
SETUP 46 
SETUP 4 7 
SETUP 48 
SETUP 49 
SETUP 50 
SETUP 51 
SETUP 52 
SETUP 53 
SETUP 54 
SETUP 55 


91 


16C 

CONTINUE 

SETUP 

56 


SAVE*l. 

SETUP 

57 


DO 16? N*2,NV 

SETUP 

58 


R A VF = VI SC ( N 1 1 ) 

SETUP 

59 


ViSC(N,l)«t Vise (N*l ,1)*R AVE+SAVEJ/3. 

SETUP 

60 

162 

SAVF «RAVE 

SETUP 

6 i 

56 

CONTINUE 

SETUP 

62 


RETURN 

SETUP 

63 


END 


o o o r> 


Sl*KfHJTt NF MIXER (FPRFS, PRECtli INF , unOT t T H 1 CK f NF, XB S T G,NS IG, I Mi) 

1 1 , THE Ti , REB ,USE P , X4 , :P4) 

DIMENSION FPRES(IOO) , THICK (24 1 , X5S IG( 100) 

FCAP ( X) =-19.556*X+107.535*X*X-336.33*X**3*508.1*X**4-295.96* <** 5 
UIK X) =- .46532* X**6 8425*X*X-.45293*X**3> *6592*X **4 
UI2fX)=-. 045929* X-i.916l5*X*X *2 .9 1843 *X **3-5 .42 I 25* X* *4 
0 1 S T -• 5 * ( X B S f G (2 )-X 9S IG < 11 ) 

XSEP=XBSIG(l)-DI$T 
XATT = XB SIG ( NSIG) *01 ST 


IF 

AT 


25 


24 


6 

7 

27 


TO 


I NOT IS NONZERO, 
SEPARATI ON. 


THE BOUNDARY LAYER IS TURBUL EM T 


CALL H4X4 ( I NOT, XSFP , DEL 1 ,T HET 1, X ATT , R EB, USEP , X 3 , H 3, X4,H4) 
IF ( X SF" P- 1 . J 24, 25 ,2 5 
C P4-0. 

GO TO 27 ' ‘ 

URAT=EXP(-.08712-UI l(H4)-.24723*{ .325 5HJI2( H4H) 

C P4 = 1 . - i iV-PREC ) /UR AT** 2 
DEAOL^XA TT-XSEP 
IF iOEADL-2 • ) 5,6,6 " 

G« ( • 5*DE ADL I **2 


G O TCP? 

G*l. 

C P4-PREC ♦( C P4— PREC) *(U-G*XSEP) 
CONTINUE 

C OE F =( PR EC-CP41 /( X ATT-X4 1 
C 7-2 • *UOOT /Ul NF 

T2 = -2 • *01 NF 

DO 20 1 , NSIG 

SUM*0. 

C0UNT*0. 


x^xTSrGTRT 


I F ( X-l . ) 
ThFTTT 


2,2,3 
ARC T I XT 


C"(TUNT*T. 

C OUNT*THET A 


TANT * S I N ( • 5*THFT A ) /COS ( .5*THET A) 

Cl 

DO 10 N= 1 , NF 
COJNT 
ANGLES 

"‘TUM* SUM¥TH"rriCrNl“*1 Z LOCUST ANGLE! 4X7*rCOOW'*TANT3 
IF J I 1 

“SLM= SUM- . 5 *C Z*T HI CKTTT" 

GO TO 35 

cT=cz*n;-xi 

XRAD=1 • / ( X* SQRT C X*X— 1 • ) ) 

CUCTFIT-1.1 

RF = SGRT( CX-i.)/(X+l.n 

SUM^TFTrc KTIT'*l(R'ffO*'iX7*T , ITF-T» ) HCZ *T"r»-T5^XR’A0Tl” 
F RAD*XR AD 

" C OUNT^T. 

DO 30 N*2,NF 

~ C CUNT =C'DUNT* 1 . 


FRAD=FRAD* XRAD 

TO 1C RT NT *F R AD*('C2 * (COUNTERFEIT* T*CTT 


'EL M I Xf. R 

1 

MI XE R 

2 

M I X-R 

3 

. MIXER 

4 

MIXc R 

5 

Ml XER 

6 

MI XF R 

7 

M I XF R 

8 

MI XFR 

9 

MI XER 

10 

MIXER 

11 

MI XFR 

12 

MI XFR 

13 

MI XFR 

14 

M I XE R 

15 

MIXER 

16 

MIXER 

17 

MI XFR 

1 8 

MI XER 

19 

MIXER 

20 

MIXER 

21 

MI XFR 

2? 

MIXER 

£ 3 

MIXER 

24 

MI XER 

Z 5 

MIXER 

26 

MI XER 

27 

MIXER 

28 

MI XER 

29 

MI XER 

30 

MI XE R 

31 

MIXER 

32 

MIXER 

33 

MIXER 

34 

MIXER 

35 

MIXER 

36 

MIXER 

37 

MIXER 

38 

MIXER 

39 

MIXER 

40 

^GLMIXER 

41 

MIXER 

42 

MIXER 

43 

MIXER 

44 

MI XE R 

45 

MIXER 

46 

MIXER 

47 

MIXER 

48 

MIXER 

49 

MIXER 

50 

MIXER 

51 

MIXER 

52 

MIXER 

53 

MI XF R 

54 

MIXER 

55 
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35 

C P -C P4 

MIXFR 

56 


IF<X-X4) 55 ,50,50 

MI XE R 

57 

50 

CP*CPMX-X4)*C0EF 

MIXER 

58 

55 

CONTINUE 

MI XER 

59 


FPRES(M) *-UI NF *C P+SUM 

MI XER 

60 

20 

CONTINUE 

MI XER 

61 


RETURN 

MIXER 

62 


END 

MIXER 

63 
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IT O 


25 
3 0 


4 


40 

41 


16 

7 

8 


SUP ROUT I NO BUBBt OFL 1 ,THE T 1 , KO 3# XC 1 ,Ul , XC5 , OCP, 1 EL 5, X, XC , MX , N 7 . , X 5 , tjQ UBB 
1 5 tUO *ALTCtRONFL ,USTOP) HIJBR 

DIMENSION X(300) , XC ( 300 I » UE ( 300 , 3 ) Su*? 

FC AP ( X ) = -l 9. 55 6* Xf 1 07.535*X*X-336.33*X**3 + 538.i*X**4-2 15 .96* X** 5 BUB 3 
urn X) = -• 46 532 *X+.6 8425 *X*X-.45293*X**3+ .6592*X **4 mm3 

UI2< X)=-.045929*X-i .91615*X*X+2.9 1843*X**3-5.42 125*X**4 BURP 

FDELT( X)=F XP<2. 5 7^3 3425 2+X-.4 379 *X*X- . 0 765 1 1* X** 3- . )0 39 7 J 7*X**4) BURR 
FA ICH { X)=FXP(-3.748i + .038 77?*X+.41967*X*X + .0 710 46*X** 3+ . 00 3 2 1 62* K* 3 UBR 
1*4) BURR 


DEL! < X) = -. 045 32 9*UQG{X)-3.9 242*X + .54535*X*X-1.39l47*X**3-i:K 
1 X**4 

FCRMAT11H1 , 44 X*31H ANALYSIS OF LE ADINC-FDGE BUBBL E////34X, IH <♦ 
1HU , 19X , 1HH » 18X , 4HDI S ?/) 

FORMAT ( 20X f 4F20»5) 

M 0UT=6 

H 1 2 5 

H5=.429 

nr 5 M = N Z , M X 

IF ( XC 1 -XC ( Ml I 4,4,5 

M 1 = M 

GO TO 6 

CONTINUE 

X1=X {Ml -“1 ] +T7CT HI T~XTNI-ITT*TXC1— X C( Ml — 1 ) ) / { XCC M 1 l-XCCMl- 1) ) 

X4-X 1-f RE NF L / IUI *REB I 

ARG=ALOG(( X4-X1 T/TRE B*nFLl *0EL1 *U 1 ) ) 

H4= . 25*FA ICH ( ARG) 

DFL4=.53*FDELT(ARG) *OEL 1 
X5=X4+10.5*DEL4* (1 . - (H4 / .429 ) **2 ) 

rFrui-usTrrPT“4TT4r»vo 

ALTL=ALTC*DELl 

IF ( X5-XI . LT* ALTLI X5=Xi+AlTL 

URA T =F XP <-. 087 12-UI1<H4)-. 24723*1 • 325 5+U I 2( H4 ) ) ) 

DCP = Ul*ijl*U.-U T RAT**2T 

ORAT=FXP (-2.243 74-FC API H4 ) + .247 23*<2 . 02 14+DEL It H4 ) ) ) 

DEL5 =DR~AT*TXFr4 

DC 7 M=NZ,MX 
IFIX5-XCMI ) 16,16,7 

GO TO 8 ' 

CONTINUE _____ _____ 

FAC T1 = 1 . -F AC T 

XC 5 = XC C M5-TT*FACTr+ XC< M5 J *T9TCT” 

U5=UF < M5— 1 ,l)*FACTi*UECM5,ll*FACT 

WRITE (MTTUT, 251 

WRT TE ( M3UT ,30) XI ,(J 1 ,Hl , OE Li 

WRITETMOUI ,30) X4 ,U 1 ,H4 , DEL4 

WRI TE t MOUT ,3 0) X5 ,U5 ,H5 , 0EL5 

RETURN " - 

FKO 


34?5*3UB1 
BURR 
19X, 1PUBK 

BURR 
RUB* 
BURR 
BUB 3 
BURR 
3 UD 3 
6 UP 3 
BU93 
PURR 
B US 3 
BUB3 
BURR 
BURR 
BURR 
BUftR 
3UR3 
B UR 3 
BURR 
BURR 
BURR 
BURR 
8U83 
BUB 3 
BURR 
BUB 3 
BUP^ 
BUBB 
BUBB 
BUBB 
BUB3 
BUBB 
BUBB 
BUBB 
BUBB 
BUBB 
BUBB 
BUBB 


1 

2 

3 

4 

c 

6 

7 

n 

9 

10 
1 1 
12 

13 

14 

15 

16 
l 7 

1 B 

19 

20 
21 
22 

2 Q 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 
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SUBROUTINE YSET ( R , A , NY, Y) 


YSET 

l 

01 ME NSI1N Y( 100) 


i YSET 

2 

RP1-1.+R 


YSET 

3 

Y 1 1 ) «0* 


YSET 

4 

Y (2 ) *A 


YSET 

5 

DO 10 N=3 , NY 


YSET 

6 

Y (N) =RPl *Y ( N-l ) - R*Y 1 N-2 ) 


YSFT 

7 

RETURN 


YSET 

8 

END 


YSET 

9 
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o u> o o 


SURROUTI NE H4X4U NOT, XI , DELI » THET 1, X5 , RFB, Ul, X3 , H3, X4 f H4 ) H4X4 

CURLF (H) =2 6. 703/H + 3 05. 03 * ALOG < H ) -2 1 i 1 ♦ 3*H+33 27. 8*H*H-2403. 9* 3** 3 H4X4 
FOELTI X) =EXP(2.5 773-.34252*X-.4379*X*X-.0765 11*X**3-.00 3 f 570 7* X**4)H4X4 
FAICHI X)=EXP(-3 .7431*.038772*X+.41967*X*X+.071046*X**3+ .0032162* X*H 4X4 

1*4) H 4X4 

10 FCRMAT<//20X,54HA SOLUTION FOR X4 COULD NOT RE OBTAINED IN 1000 TRH4X4 
1IAL5) ~ H4X4 

M0UT=6 H 4X4 

H4X4 

IF I NOT IS NONZERO, THF BHUNOARY LAYER IS TURRULFNT H4X4 

AT SEPARATION. H4X4 

H 4X4 

IF (I NOT) 2 ,5,2 H4X4 

2 H3*THE T1 /DELI H4X4 

X3*X 1 ~ H 4X4 

DEL3*DE L 1 H4X4 

GO TO 20 " “ " H4X4 

5 X3=X1+5.F4/(U1*REB) H4X4 

ARG* A LOG (TX3-X1 ) /IRE B*D £ LI *0E LI ) I H4X4 

H3 = THET1 *F AT CO ( ARG) /DELI H4X4 

DE L 3*. 59*FDELT I ARG) *DEL 1 H4X4 

I F ( X3-X5 ) 20,15,15 H4X4 

15 H4=.4l^ H4X4 

X4*X5 H4X4 

GO TO 5 0 ' ' " """" H 4X4 

20 CONTINUE H4X4 

l GO* 6 * ' H4X4 

D I ST-X5— XI H4X4 

UNDER*!). " H4X4 

H4*H3+H3 H4X4 

C OF FT *D t L 3 *H 3 " H4X4 

COEF 2*1 0.5*0EL3*H3 H4X4 

= F 1 *C UR LF ( H 31 H4X4 

95 0 VE R*H4 H4X4 

“ H4=.5*M4*UN0ER) ’ H4X4 

X 4*CUR LF (H4)*C0EF1»SU B H4X4 

I GO* IGO+1 * H4X4 

IF I XV^LTCR) *I _ ,507S2 H4X4 

41 IFUGO-IOOO) 95,61,61 H4X4 

T*Z I F < AB5 1 X4-ALTFR ) /HI 5T-. 001 I 50,50,43 H4X4 

43 UNDER=H4 H4X4 

H4=^ 5* ( 0 VER+H4 ) ' ~ H4X4 

X4*CURLF<H4)*C0FFUSUB H4X4 

~AETE R*X5-C DEFT - * TI •- IH4/ .429 1 **2T7fl¥ H4X4 

IGG*IGO+ 1 H4X4 

IF I X4-ALTE R ) "57, 50,51 " H4X4 

51 IF ( IGO-IOOO) 43,61,61 H4X4 

52 I FI ABS( X4-ALTFR) /DI5T-.0D1 ) 50,50755 H4X4 

61 H4*. 429 H4X4 

' X4*X5 H 4X4 

WRI TE ( MOUT ,10) H4X4 

50 CONTINUE : H'4X4 

RETURN H4X4 

FND " ’ H4X4 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
Z, 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 
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SURRCUT I NE SETSX (MS PI , XS EP t X ATT , XS IGf ANGLE) 
DIMENSION XSIGCIOO) 

A«.5*(XSEP+XATT) 

B s » 5» ( XATT— XSEP) 

ARG*0. 

DO 5 N«l»NSPl 
X SIG (N) *A-B*Cf)S ( ARG) 

5 ARG*ARG+ ANGLE 
RETURN 

ENO _ 


SETSX l 
SFTSX 2 
SETS X 3 
SETSX 4 
SETSX 5 
SETSX 6 
SFTSX 7 
SETSX 8 
SETSX 9 
SETSX 10 
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FUNCTION ARCTTXT 

ARCT 

1 


PI -3 • 141 59 

ARCT 

2 


IF (ABSI XJ-I.F-61 1,2,2 

ARCT 

3 

l 

ARCT = .5*PI 

ARCT 

A 


GO TO 6 

ARCT 

5 

2 

IMX+. 99999) 3,4,4 

ARCT 

6 

3 

ARC T*P I 

ARCT 

7 


GO TO 6 

ARCT 

3 

4 

A RC T * A T 4 N ( S O R 1(1 • - X * X l? X i 

ARCT 

9 


IP(ARCT) 5,6,6 

ARCT 

10 

5 

ARC T *ARC T4- PI 

ARCT 

11 

6 

CONTINUE 

ARCT 

12 


RETURN 

ARCT 

13 


END 

ARCT 

14 
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FUNCTION 3AMi(ACAP,DXI f PI) 

OIMFNSION AC AP ( 30 »3 ) 

GAMl=PI*(~i # 5*ACAPfl»l)~.75*ACAPI 2 , 1 ) + 2 . * ACAP | i 
1 A P ( l f 3 ) - , 2 5 * AC A P 12, 3) )/D XI 
RETURN 
END 


GAMI 

GAML 

f 2 ) +ACAP ( 2, 2)-.5*4CGAMl 
GAMI 
GAMI 
GAMI 


i 

? 

3 

4 

5 

6 


100 




FUNCTION FBIXl ,X2,X3,Y) 

F R 

1 


01 *1 • / l X2-XI ) 

FH 

2 


D 2= 1 • / ( X3— X2 ) 

FR 

3 


Tl =ABS( Y-Xl ) 

F ft 



T2 * A B S ( Y- X2T " ' " ' "" 

FR 

5 


T3*ABS( Y-X3 ) 

FO 

6 


EP$=l.E-6 

FR 

7 


IF ( Tl-F PS) 2,3,3 

FB 

8 

2 

Fl=0. 

FO 

9 


F2=A10G l T2 ) 

FB 

10 


F 3=AL0G ( T3 ) 

FB 

11 


GO TO 10 

FB 

12 

3 

Fi=ALOG(Tl) 

FO 

13 


IFIT2-EPS) 4,5,5 

FB 

14 

4 

F 2*0* 

FB 

15 


F3*A LOG ( T3 ) 

FO 

16 


GO TO 10 

FB 

17 

5 

F2=ALOG(T2) 

FR 

18 


IF (T3-EPSr 6,7,7 

FB 

19 

6 

F 3-0 • 

FB 

20 


GO TO 10 ‘ 

FB 

21 

T 

F 3=A LOG ( T3 ) 

FB 

22 

10 

FB=( (Y-xl)*FI*DUf 01*02 )*(X2-Y)*F2»(Y-X3f*F3*D2 173.14159 

FB 

Z j 


RETURN 

FB 

24 


END ‘ ' 

FB 

ZZ 
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O' CO 


6 


7 


SUBROUT I NE FGAMMMU« NG, A , 8, XS EP, X ATT, GAMMA, Y , GI I EG AM I 

DIMENSION A ( 30 ,3 ) EGAMI 

SI NT -SQRTd .-Y*Y> EGAMI 

THE T A*ARCT t_Y) EGAMI 

SUM*0.‘" ' ' EGAMI 

CCUNT=l • _ __ EGAMI 

no 6 N-2 ,NG EGAMI 

CCUNT«C OUNT> I • EGAMI 

SUM= SUM* A ( N^i , NU) *f $ I Nf ( COUNT +1 . ) ♦THET A) /{ COUNT +1 . )-SlN ( (CO JNT- l. )EGAM I 

I* THE TAJ /(COUNT- UJJ EGAMI 

GI*( 3* 14159— THE TA^SI NT) * C A(1 YNU ) ♦«5* A( 2 , NU ) > ♦ . 5*SUM- . 25*GAMMA*( W ♦EGAMI 
1 Y) * SINT* SI NT EGAMI 

IF(Y-XATT) 8,8,7 EGAMI 

DIFF«1.-XATT _ EGAMI 

IF (0IFF-1.E-6I 8,879 EGAMI 

GI*GU2.»B*r)IFF ♦*(-!. 5) *SQRT( ( X ATT^XS FPI*(l«-Yl*< Y-XATT ) ) EGAMI 

Continue " egami 

RETURN _ EGAMI 

END ~ - EGAMI 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 
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SUBROUTINE ESIGI (NU , M » X S f B v Y f S II 
DIMENSION X$tlOO)f*mO,3) 

SUM = 0. 

DC 1C 1=2, NX 

10 suM=suMfR(i ,NU)*GB(xsn-n , xs i n , xsn ♦ n , y i 

SI=*U,NU1 *Rl NT ( XS ( l J » X S ( NX* l I » Y 1 ♦SUM 
return 

END 


r SIS I 1 

ESISI 2 
ESISI 3 

esisi 4 

ESISI 5 
ESISI 6 
ESISI 7 
FSI3I fi 
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FUNCTION GBIXl»X2,X3,X» GB l 

GB=ABINT(X1 »X2 » X J - A8I NT (X3»X2»X ) GB 2 

GB-GB/3.H159 GB 3 

RF TURN . . GB * 

F NO ” GB 5 
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PUNCTIf N ABINTI A,n,X) 

ARINT 

1 


ARG A = AR S ( X— A ) 

A H I Nl T 

p 


ARGR=ABS<X-B> 

ABIN T 

3 


CPE C =2 . * ( B-A) 

ABINT 

4 


AP1=A+1. 

ABINT 

5 


BP1=R*1. 

ARINT 

6 


IF (ARGA-l.E-6) 2,1,3 

ARINT 

7 

2 

CA = 0. 

ARINT 

ft 


GO TO 5 

ABINT 

9 

3 

CA=ALOG( ARGA) 

ABINT 

1 0 


I F ( ARGB-1# F-6)~ 4,5,5 

ARINT 

I 1 

4 

CR = 0. 

ARINT 

12 


GO TO 6 

ABINT 

13 

5 

CB=AlOGl ARGB) 

ABIN T 

14 

6 

ABINT=(CA-.5)*ARGA**2-( CB-.5) *ARGB**2-( ALOG< API . 

5 )*AP l**2MALOC,< ABINT 

15 


1BP11-.5) *BP1**2-C0EF*(( X-BJ*(CB-1 .)+BPl*l ALOG1 BPl > 

- 1 • > ) ABINT 

16 


A8INT=AB!NT/COEF 

ABIN T 

1 7 


RETURN 

ABINT 

ia 


FND 

ABINT 

19 


105 



FUNCTION BINTI XS ,XZ,X) 

BINT 

1 


RTS=SQRT II , + XS) 

BINT 

2 


RTZ = SQRT{U + XZ) 

BINT 

3 


RINT=-i.-X+RTS*RTZ 

BINT 

4 


IFIXZ-X) 2,3,3 

BINT 

5 

2 

RTSX = SQRTt X-XS) 

BINT 

6 


RTZX- SQR T ( X-XZ) 

BINT 

7 


BINT =BI NT+ I XZ-XS) *ALCGC t RTSX+RTZX ) / ( RTS +RTZ ))+RTSX*RTZX 

BINT 

8 


GO TO 50 

BINT 

9 

3 

IFIX-XS) 5,5,4 

BINT 

10 

4 

BINT*BINT+ ( XZ-XS ) *A LGGl S QRT (XZ-XS ) / ( RTS+RTZ » ) 

BINT 

11 


GO TP 50 

BINT 

12 

5 

RTSXxSQRTI xs-xi 

8 I NT 

13 


RTZX*SQRT< XZ-X) 

BINT 

14 


BINT*BI NT+ ( XZ-XS) *41CG( (RTSX+RTZX)/! RTS+RTZ ) )-RTSX*RTZX 

BINT 

15 

50 

continue 

B INT 

16 


RETURN 

BINT 

17 


END 

BINT 

18 
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SUBROUTINE SCltTSBLtNSBL.FRZ, ARR, ROBB I 
DIMENSION SSL (300) 

0ELZ«FRZ*RD8B 

£N«ARR/FRZ 

00' 5 'N- 1*300 ' “ 

I F IFN-N) 4,4,5 

4 NE«N 
GO TO 6 

5 CONTINUE 

6 NG-NSBL-NE 

EN«FLOAf (NG'I * . 

NGMl*NG-l 

SBL(1)*0." 

00 7 N*2,NE 

7 SBU NI ■SBL I N-U ♦ CELZ 
FRACT-2.2/DELZ 

FR AC 1«F RAC T-I7 

R*FR AC T** 1 1 , /FLOAT ( NGMl ) ) 

B SAVF-R 

R*R— I R**NG—FRACT*R'f FRAC 11/ (EN*R**NGMI -FRACT1 
I F I A BITS AW-R I -T.ir-6 r 979.8 

9 RPl*Rf 1 . 

00 10 N-NFVN5BL 

10 SBL I N*1 ) »RRl *SBL(N) “R*SBL IN-1 ) 

RETURN" * 

■ END 


SCAL 1 
SCAt 2 
SCAL 3 
SCAL A 
SCAL S 
SCAL 6 
SCAL 7 
SCAL A 
SCAL 9 
SCAL 10 
SCAL 11 
SCAL 12 
SCAL 13 
SCAL 14 
SCAL 15 
SCAL 16 
SCAL 17 
SCAL IB 
SCAL 19 
SCAL 2 0 
SCAL 21 
SCAL 22 
SCAL 1 * 
SCAL 29 
SCAL 25 
SCAL 26 
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SUBROUTINE TERPF t XI , J ,T ABl,T A82 , T AB3, T AB4, X l TAB , FP ) 

TERPF 

1 


DIMENSION TABl (24) , T AB2 ( 24 ) , T AB3 < 241 , T AB4< 24 ) , X IT AB ( 24) 

TERPF 

2 


Iff XI-.00011 2*2,10 

terpf 

3 

2 

GC TV (3,4,5,61 ,J 

TERPF 

4 

3 

FP=2.53-2.439*4L0G( XI) 

TERPF 

5 


GO TO 99 

TERPF 

6 

4 

FP=3*54-1*725*4L0G( #707 l*Xl ) 

TERPF 

7 


GO TO 99 

TERPF 

B 

5 

FP*4.5R-l. 2 195* A LOG t »5*X I ) 

TERPF 

9 


GO TO 99 

TERPF 

10 

6 

F P*1 0. 1 2 

TERPF 

11 


GO TC 99 

TERPF 

12 

10 

DO 12 N*1 ,24 

TERPF 

13 


IF t X I-XI TAB ( N) ) 11,11,12 

TERPF 

14 

1 l 

NX*N 

TERPF 

15 


GO TO 13 

TERPF 

16 

12 

CONTINUE 

TERPF 

17 

13 

TX*(XI-XITAB( NX- 11 ) / (XIT AB (NX l-X IT AB( NX -11 1 

TERPF 

18 


TX1 = 1.-TX 

TERPF 

19 


GO TO (14,15,16,17) , J 

TERPF 

20 

14 

FP=TXI*TAB1(NX-1H-TX*T ABM NX » 

TERPF 

21 


GO TO 99 

TERPF 

22 

15 

FP»TXl*TAB2(NX-l!*TX*TAB2(NX) 

TERPF 

23 


GO TO 99 

TERPF 

24 

16 

FP=TXl*TAB3f NX-l H-TX9T AB3 (NX » 

TERPF 

25 


GO TO 99 

TERPF 

26 

17“ 

FP=TXl*TAB4(NX-ll+TX*T AB4 ( NX ) 

TERPF 

27 

99 

CONTINUE 

TERPF 

28 


RETURN 

TERPF 

29 


END 

TERPF 

30 
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SURROUTI NE E VSUl NNF »XX»SSC,SST, C.CR»TT 8, CCM, TTM } 

EVAL 

l 


DIMENSION SSC (50) *SS'i (50 ) 

EVAL 

2 


COST * 2 •* X X - U 

EVAL 

3 


COSTS * COST*'*? 

EVAL 

4 


IFICfSTS-l.E-m 303,304,304 

FVAL 

5 

304 

TANT = SORTU. /COSTS - l.J 

EVAL 

6 


THF * ATAN(TANf) 

EVAL 

7 


GO TO 305 

EVAL 

8 

3 03 

THE = 1,5708 

EVAL 

9 

3C5 

I F (COST) 403,404,404 

FVAL 

10 

403 

THF = 3.14159 - THF 

EVAL 

11 

404 

ARG * 0. 

FVAL 

12 


SUM! = 0# 

EVAL 

13 


SUM2 * 0. 

FVAL 

14 


00 551 N*i * NNF 

EVAL 

15 


ARC = ARG ♦ THF 

FVAL 

l 6 


SU«l * SUMl «■ SSC(N) *51 N ( ARGI 

EVAL 

17 

551 

SUM2 = SUM2 ♦ SST(N) *SIN( ARC) 

FVAL 

18 


r.CR = SUMl * STN< THE 1 *f.CM 

EVAL 

19 


TTB * (l. - COS (THE! >*SIJM2*TTM 

FVAL 

20 


return 

EVAL 

21 


END 

EVAL 

?? 


i 


{ 


I 
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SUBROUTINE SI MP < NS , OX ,0R0 , FI NO) 

SIMP 

I 

DIMENSION ORD ( 50 1 


SIMP 

2 

INTEGRATION OF NS ♦ 

1 EQUALLY SPACFO ORDINATE VALUES 

SIMP 

3 

BY SIMPSON'S RULE. 

^NS MUST BE EVEN 

SIMP 

4 

SUM » 0. 


SIMP 

5 

DC 88 I *2 1 NS t? 


SIMP 

6 

SUM * SUM ♦ 2. *010(1-1 1 

♦ 4.*oro<i) 

SIMP 

7 

FIND » DX* ( SUM - URDU) 

♦ GROINS +1 ) )/ 3 • 

SIMP 

$ 

RETURN 


SIMP 

9 

ENP 


SIMP 

10 
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SUn ROUTINE SECT (XU, YU , XL , YL, NOFF , NF, RCDRC, TM AX, CMAX,ST, SC) SECT 

PROGRAM TO COMPUTE COEFFICIENTS TN AND CN OF THE FOURIER SERIES SECT 

REPRESENTATION OF SECTION THICKNESS AND CAM8ER DISTRIBUTIONS SECT 

OIMFNSION XU I 00) t YU ( 30 ) , X L 130 ) , Y L ( 30 ) , YUC( 30 ) , YL C< 3 0 ) , S T< 24 ) • SC ( ^SFCT 
1) ,OUM(50) ,TBAR(5(j> ,CPAR(50) SECT 

12 FORMAT ( / ///47X,26HI NPUT AND COMPUTED OFFSFTS/ ) SECT 

13 FORMAT (19X,4HX1/C,12X ,4HYU/C , 11 X, 5HYUC/C , 20X , 4HX 1/C, 12X ,4HYL/C, 11XSECT 


1 , 5H Y LC /C /) 

14 FORMAT! 3X,3F16.5,0X,3F16.5) 

N A-6 
R NA=6 • 

ft NF = r L OAT( NF) 

MCUT =6 

PI = 3* l 41 59 
DELT = PI/ ( 2 • ♦R NF ) 

NTC = 2*NF - l 

N I N T = NTC ♦ Y 

NSIMP = NTC ♦ 1 
RDRC =• 5 *RCDBC 
VARY = 0. 

CB = 0. 

TB * 0. 

THFTA=~0. 

DO 89 K = 1 , NTC 

thfta = thettt cflt “ 

Xi = . 5* ( 1 • * CaSJTHFTAn 

on 90 L A M= 2 , NOF F 

IFIX1-XUILAM) ) 110,90,90 

11C YU TNT - = YU I LAM-1 ) 4. f xl X\3 T LAM- 1 ) l* I YU ( L AM ) 

1) - XUILAM-IH 

gc re in 

SC CONTINUE 

111 no 80 LA M *=2 ,~NOF F" ■ — 

I F ( Xl-XL ( LAM) ) 210,80,80 

210 YITNT = YL(LAM-l) + (XI - XL( LA~R-TT f * ( YCTLAM j 
1) - XL t L AM— 1 ) ) 

GC TTTTTT 

80 CONTINUE 

112 TR A*TK+ IT = .^lYlJVWT YTr NT F 
89 C BAR ( K+ 1 ) = • 5* ( YUI NT + YL I NT ) 

" TMAX * 0, 

C MAX = 0. 

no~79T<T=* 2, NS IMP " 

IF (TBAR( K) — TMAX) 001,802,802 

eC2 TMAX = TBARIKJ 

8C1 IF (C BAR ( K) -C MAX) 79,702,702 

TO? CMAX = CBAR(K) 

75 CONTINUE 

1201 C MA X = l . 

1202 CO NTT RUE" 

I F ( T MA X- 1 *F -5 ) 1140,1141,1141 

1140 TMA X=l. 

1 141 DO 69 K=2, NSIMP 


SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

SECT 

YU(LAM- 1) ) / { XU( LAMSECT 
SECT 
SECT 
SECT 
SECT 
SECT 

YL ( LAM- 1 ) ) / ( XL( LAM SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 
SECT 

’ SECT 

SECT 
SECT 


1 

2 

4 

5 

6 
7 
P 
9 

10 

11 

12 

13 

14 

15 

16 
1 7 
i R 
i 9 
20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 
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69 

CBAR(K) = C.BAR(K)/CMAX 

SECT 

56 


TBARd) a 0. 

SECT 

57 


CRAR(l) = 0. 

SECT 

58 


TRAP (MINT) a 0. 

SECT 

59 


CBARtNINT) * 0. 

SECT 

60 


TTA * TRARINA1 

SECT 

6 1 


TTB * TB AR ( NA+1 ) 

SECT 

62 


TTC - TB AR ( NA4-2 ) 

SECT 

63 


TAA * DE LT* (RNA-l . ) 

SECT 

64 


TBB * TAA 4- DELT 
TCC = TBB + DELT 

SECT 

SECT 

65 

66 


XA = ,5*C0S!TAA) 
XB = ,5*C0S( TBB) 

SECT 

SECT 

67 

68 


XC * • 5*C0S ( TCC ) 

SECT 

69 


SLOPf * 1 ( TTC-TTBI * ! XB-X A) / IXC-XB ) ♦ <TTB-TTA)*I XC-XB)/ ( XB- 

XA II / ( XSECT 

70 


1C-X4I 

SECT 

71 


THETA - 0. 

CO SB * COS (TBB) ) 

SECT 

SECT 

72 

73 


00 456 I *2 »NA 

SECT 

74 


THFTA * THETA ♦ DELT 

SECT 

75 


COST * COS? THETA) 

SECT 

76 

456 

TBAR(I) * ( SORT tl. -COST)/! l.-COSB ) **1 .5 1 *1 TT B*l l.+COST- 2..C0S8) /II SECT 

77 


U-COSB) ♦ • 5*SLDPE* ( COS T-COS B 11 

SECT 

78 


NLE * 2*NF ♦ 1 — NA 

SECT 

79 


COSRl « 1. 4- COS ( PI -RN4*DELT 1 

SECT 

80 

— 

THETA * PI 

SI NA S*SI N(RNA*DELT) **2 

SECT 

SECT 

81 

82 


COSAS=COS( RNA*OELT) 

SECT 

83 


ANG*0. 

SECT 

84 


DO 457 l*2 f NA 

SECT 

85 


I NO * 2*NF ♦ 2 - I 

SECT 

86 


THETA » THETA - OELT 

SECT 

87 


CCST1 *1*4- COS ( THETAI 

SECT 

88 


AKG*ANG4-DELT 

SECT 

89 

457 

COEF*l SI NAS-SI Nf ANG) **2 ) / ( C0SR1*! COS! ANGI+COSAS ) ) SECT 

TBAR ! f NO) * ! S ORT ! ROBCCOSTl ) *COEF/TM AX ♦TBAR! ML E )*! COST 1/CO SR l) **l SECT 

90 

91 


l. 51/(2. -COSTl) 

SECT 

92 


THETA « TAA 

SECT 

93 


NAPl * MA ♦ 1 

SECT 

94 


00 458 1 * NAPl t NLE 

SECT 

95 


THETA a THFTA ♦ OELT 

SECT 

96 

45e 

TBAR ( I ) * TBARUI/U.-COSCTMETAI) 

SECT 

97 


THETA * 0* 

SECT 

98 


00 459 1 *2 f NSI HP 

SECT 

99 


THETA » THETA ♦ OELT 

SECT 

100 

459 

C BAR II ) a CBAR! ( I/SINITHETA) 

SECT 

101 


RKK * 0. 

SECT 

102 


00 59 K*i,NF 

SECT 

103 


RKK * RKK ♦ i # 

SECT 

104 


THETA * 0. 

SECT 

105 


So m 1*1 »N!NT 

DUH(I) » TBAR! 1 1 *SINITHET A*RKKJ 

SECT 

SECT 

106 

107 


TFT' THE T A - TUFT A ♦ CELT SECT 108 

CALL SIMP! NSI MP»0ELT tDUMtV ARY ) SECT 109 

yriR r~rr7*VA*Ym sect uo 


112 


THFTA =0. 

0 0 8 88 I *1 ,NINT 

DIIMIIJ * C.BARU) *SIN!THETA*RKK) 

6P8 THETA = THETA + CFLT 

CALL SI H P C N S f “ P » BE Lt , DU M , V A R Y » 

59 SC(K) * 2 • * VAR Y/ PI 
DO 969 1*1, NOE F 
X * XUIII 

CALL FV4L!NF,X,SC,ST»CB,TB,CNAK,TMAX I 
565 YUC II) * CB ♦_ TB 

f)0 869 l *1,NI0FF 

x * xnn 

CALL EVAL!NF,X»SC»ST ,CB,T B ,CMAX , TMAX I 
£65 Y LC (II * CB - TB 
SUM1 * 0. 

C CUNT * 0. 
on 659 T *£ ,NF 
COUNT * COUNT ♦ 1. 

659 SUMl * SUM1 - ST CD *COUNT*l— l , ) **I 
RCDBC = 8.*<TMAX*SUMl)**2 
RC3BC=2. *RCOBC 
TMA X*2 • * TMAX 

C x «2 . *C'H A X 

WRITE! MOUT ,121 
WRITE! MOUT »13ji 

WPI TF I NOUT ,141 (XU!!) ,YU( I), YUC! I ) ,XL III t YL! I) , YLC! I), 1 = 1, N IFF > 
RF TURN ~ 

END 


SECT 111 
SECT 11 
SFCT U 
SFCT 11 
SFCT 11 
StCT 11 
SFCT 11 
SFCT 11 
SECT 1 1 
SFCT 12 
SECT 12 
SECT 12 
SFCT 12 
SECT 124 
SECT 12 5 
SECT 126 
SECT 127 
SECT 128 
SECT 129 
SECT no 
SECT 131 
SECT 13' 
SFCT J3 
SECT 134 
SECT 1C 5 
SECT 136 
SECT 137 
SECT 138 


i 
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SUBROUTINE COROX (NS BL,NZ , ROBB ,S RL , X, X C ) 

BOUNDARY LAYER COOR 01 NATES AND CORRESPONDING CHORD AL 
COORDINATES ARE COM PUTED HERE * 

DIMENSION SBL 1 310) , X (300 ) , XC ( 330 ) 

56 FORMAT(//lOX f 3lHITERATI ON TO COMPUTE XC FOR M = 15,32H 
IERGE IN 1000 STEPS,) 

3 7 FORMA T( IH1 ,25 X, IHM f 20X , 1HS ,25X, IHX , 2*X , 2HXC/ / ) 

3 8 FOR M A T ( 2 2 X,I5,3E25, 5) 

MOUT *6 

MX = NS8L ♦ NZ - l 
R ZERO * RDB8/2. 

XC(NZ) = -l, 

DO 255 M»I , N Z 

MM = NZ ♦ I - M 
55 X(M) * SBL(NZ) - SBL(MM) 

DO 256 M*NZ , MX_ 

MM * M f I - NZ 

56 X ( M ) - SRL(NZ) ♦ SBL(MM) 

00 2 57 Mil , M X 
IFINZ-M) 333,257,335 

33 K '* M ♦ l - NZ 
GO TC 33* 

35 K '* NZ - M ♦ l 
3* XC ( M) * -l* * SBL(K) 

TF(SBL(K)~RZERO) 341,3*1,3*2 
*1 XC ( M) » -l • ♦ SBL ( K) **2/ (*,*RZERO) 

*2 Continue 

00 258 L*l ,1000 
SAVE * XCIM) 

CALCl = $QRT((1.+XC (M)) /RZERO) 

C ALC 2 *^SQlTnTT( l • ♦XC ( Mf i/RZ ERO) 

XC(M)«XCCM) + CALC1*(SBL( K) - RZFRO*( CALCl *CALC2 ♦ALOGC 
l) /CALC 2 

I F ( ABS( SAVE-XCI Ml )-l .E-6 ) 257,257, 258 
58 CONTINUE 

WRITEIMDUT ,336) M 

57 CONTINUE 
WRITEIMOUT ,337) 

DO 26* M=1 , MX 

I F ( NZ— M) 261,261,262 

62 k*NZ-Md 
GO TO 263 

6 1 K»M«-1-NZ 

63 WRITEIMOUT, 3381 M, SBL IK) ,X(M) y XC(M) ________ 


CORD X 1 
CORD X 2 
CORD X 3 
CORD X * 
CORD X 5 
CORD X 6 
DID NOT CONVCORD X 7 
CORD X 8 
COROX 9 
CORD X 10 
CORD X 11 
Ct)RDX 12 
CORD X 13 
COROX l* 
CORD X 15 
CORD X 16 
CORD X 17 
CORD X 18 
CORD X 19 
CORD X 20 
CORD X 21 
CORD X 22 
CORD X 23 
CORD X 2* 
CORD X 2 5 
CORD X 26 
CORD X 27 
CORD X 28 
CORD X 2 9 
CORD X 30 
CORD X 31 
CORD X 32 
CORD X 33 
CALCl ♦C4LC2)IC0RDX 3* 
CORDX 35 


RETURN 
EW 
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SUBRPUTI NF PG1TZ?)(M, X,lJE f DX I f PRESS,SA, SB,SCtSR,SS) 


PGRAO 


SUBROUTINE FOR CALCULATION OF PRESSURE GRADIENT AND 
DERIVATIVE COEFFICIENTS. 

DIMENSION XHOO) ? UE(300,3) 

0 1 2 = X ( M4* l ) — X f M ) 

D2Z=X(M*2)-X(M) 

D21 = X< M + 2 ) — X (M-f l ) 

DIMI^XI M+l ) — X [ M— 1) 

DZM1=X(MT-xH-1) 

X I M^D I Z / <02 2*02 I ) 

FTAM=i. /OlZ-i./02i 
ZETAM-D2 1/(DIZ*02Z) 

PRESS = (3 • *UE ( M + l , 1 ) -4. *U EC M«-l y2)^UE{M^l f 3))/( 2 
1 XI M*IJE ( M+2 1 1 ) +E TAM*UE ( M+l 1 1 ) -ZET AM*JE<M, I) ) 
SA*1./01 Z+l./OLMl 
SB=DIMi/(0IZ*0ZMl) 

SC *0 1 Z / ( DlMl *6 ZMi j 
SR=DIMI/0ZM1 
S S =0 1 Z /D ZM I 
RETURN 

W " 


PGRAD 

PGRAO 

PGRAO 

PGRAO 

PGRAO 

PGRAO 

PGRAD 

PGRAD 

PGRAD 

PGRAO 

PGRAD 

PGRAO 

PGRAD 

*DX n + UE( M + lt 1 )* ( PGRAD 
PGRAD 
PGRAD 
PGRAD 
PGRA 0 
PGRAD 
PGRAD 
PGRAD 
PGRAD 


l 

? 

3 

4 

5 

6 
7 
R 
9 

1 0 

11 

12 

13 

14 

15 

16 

17 

18 
L 9 
20 
21 

2 7 
2 3 
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SURROUTI NE TRANS (UPR l M, PRESS , THET A f RE B , UC , NY , FL AM f X FLAM , L AM J ) 

SLBRnUTINE TO TEST FOR TRANSITION IN A L AM INAR BOUNOAR Y LAYER, 

0 I MF A Si W UC ( 1 6 0 » 3 ) fTTAMlO I , XF L AM ( LO ) 

F ( X > = • 1 1 746 - l.0582E~3*X - 1 • 1023E-4*X*X 
TKAY = P RE SS*RE P *TH FT A^*2/UC { NY , 2 ) 

IF CTKAY-.077) 2,2,99 

2 IF I A B S ( TKAY)— ,0101) 3,3,4 

3 ARG = TKAY»7 2.4B 

GO TO 5 

4 ARG = 0. _ ' • 

DO 6 N=l ,1000 

SAVE = ARG _ 

ARG = A RG - ( A RG * c ( ARG ) **2-TKAY) / ( F< ARG1*( • 1 1746- ARG*3. l 746 E- 3 
lRG*AR5*5.5115E-4 > ) 

~ IF<ABS(l.-SAVE/ARG)-i.E-6) 7,7,6 

6 CONTINUE 

7 IFfARGf II. ) 8,8,5 

8 EF = 1 .75 

GO tn 10 

5 _ DO 15 N S 1 , 1 0 

IF ( ARG-XFLAM (Nl ) 24,24,15 
24 N8AR = N 

GO TO 16 
15 CONTINUE 

16~ EF = FLAM( NBAR-l ) «• ( ARG-XFLAM( NBAR- 1) )*( FL AM ( NBAR ) -FLAM( NBAR- l ) ) 
IF LAM ( NBAR) -XFL AM ( NB AR-l ) ) 

10 ' B * .5*EF 

A = 3 .36 * C UPR I M/UC ( NY ,2 ) ) **2 
RTH * F( ARG)*(SQRT<B*B*9860.*A1-BI/A 
I F ( REB*THETA-RTH) 99,50,50 
50 LAMQ = 0 

95 CONTINUE 

RETURN 

END 


TRANS l 
TRANS 2 
TRANS 3 
TRANS 4 
TRANS 5 
TRANS 6 
TRANS 7 
TRANS 8 
TR AN S 9 
TRANS 10 
TRANS 11 

tTrans 12 

TRANS 13 
TRANS 14 
- A TRANS 15 
TRANS 16 
TRANS 17 
TRANS 18 
TRANS 19 
TRANS 20 
TRANS 21 
TRANS 22 
TRANS 23 
TRANS 24 
TRANS 25 
TRANS 26 
/(XTRANS 27 
TRANS 28 
TRANS 29 
TRANS 30 
TRANS 31 
TRANS 32 
TRANS 33 
TRANS 34 
TRANS 3 5 
TRANS 36 
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SUBROUTJ NE C AF§‘ntFRfN f CAPG,CAPH f CAPj f CAPK f SRf SSt SD,SE, S F f V 
1C) 

D l ME A SI HM CAPGUOO) , CAPH ( 100 ) , CAPJ I 100 ) t C APK ( 10 0) 

0 1 MENSI ON VI SC (100 ,2) ,V ( 1 00, 2 ) , UC( 1QQ , 3) , SOI 100 ) y SE( 100 ) f SF 
IF ( ME R) 4,2,4 

C APG (N) 3 SR*V ( N , 1 ) - SS *V ( N, 2 ) 

C A PH ( N f * SR* VlSC < N , l I -SS * V l SC ( N, 2 1 

CAPJ(N) =SR*(SO(N)*VISC(N*l ,l)*SE(N)*V ISC(N,ll-SFm*VISC(N- 
IS*(S0(N)*V I S t ( N*1 , 2 ) *SE ( N ) *V ! SC( N, 21 -SE I N ) *V I SC ( N- 1 , 2 ) ) 
CAPK(N)* SR*UC (N,2)-SS*UC(N,3) 

GO TO 6 

CAPG(N)*.5*JCAPT,(S) *V(N,il) 

C APH ( M ) « . 5 * ( C APHTnI W IS C C N, 1 ) I 

CAPJ(N)».5*(CAPJ(N)+S0(N)*Vl$C(N*l, i ) *SE( N > *V I S C( N , i)-SF(N) 

i-i, m 

CAPMN)*.5*(CAPMNH.UCCN,1I) 

‘TONTTWJE 

RETURN 

END 


ISC, V,UCAPS 
CAPS 
CAPS 
(100) CAPS 
CAPS 
CAPS 
CAPS 
1,1) )- SCA PS 
CAPS 
CAPS 
CAPS 
CAPS 
CAPS 
*VI SC ( NC APS 
CAPS 
CAPS 
CAPS 
CAPS 
CAPS 


1 

? 

1 

4 

5 

6 
T 
8 
9 

10 

11 

12 

13 

14 
1 5 
16 

17 

18 
19 
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SUBROUTINE TERP l Yl N , YBAS E ,V ARY, NY ,V ALUE I TERP 

r. TERP 

r. SUBROUTINF FOR DETERMINING INTERPOLATED VALUE OF THE TERP 

C FUNCTION VAR V AT Y = YIN. TERP 

C TERP 

DIMENSION YBASEIIOO) .VARYtlOO) TFRP 

IFIYtN-YBASEtNY-in 2,3.3 TERP 

3 VALUE * VARYINYI TERP 

GO TO 10 TERP 

2 DO 15 N*l,NJT _ TERP 

IFIYIN-YBASECNII 24,24,15 TERP 

24 NBAR*N TERP 

GO TO 16 TERP 

15 CONTINUE TERP 

16 021*YBASEINBARI-YBASE<NBAR-ll TERP 

D31 *VBASE < NBAR»1 )— YBASE I NBAR-l > TERP 

032*031-021 " TERP 

D3A*YBASEt NBAR^l I- YIN TER® 

D2A* YBASEI NBARI— Yl N TERP 

0A1*YIN- YBASE(N3AR-H _ _ TERP 

VALUE=D3A*D2A*VARY( NBAR-l I / ( D21*D3 1 MD3 A*OAl*V AR Y l N BAR I /1 02 1*032I-TERP 

ID2A*DA1*VARY(N8AR+1 )/<D31*D32) TERP 

10 CONTINUE ” TERP 

RETURN TERP 

END ' " '■ TERP 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 

ie 

19 

20 
21 
22 

23 

24 

25 


118 


SljB ROUT I NE YD I FF ( NY , ALPH A , BET A, GAMMA, DELT A » SD f S E , SF , C 2, C 3, C 4 , Y) 
DIMENSION ALPHA (100) f '»FT A( 100 ) , GAMMA( 100), DELTA! 100) 

OTMENSION SD(IOO) ,SF(1Q0),$F< 100),Y(100) 

N V-NY-2 


YDI FF 
YDIFF 
YDIFF 
YDIFF 
YD 1 c F 
YDIFF 


NVP1=NV+1 
DC 40 N=2,NV 

ALPHA (N) = 2.*(2.*Y(N)-Y<N-1 )- 1 { N + l) )/( (Y(N + 2)-Y(N~l))MY(N*-2)-Y(NY0IFF 
l*lll*IY(N»2)-YI_NI)l YDIFF 

DELTA (N) = 2>mN+2)*Y(N*l)-2.*Y( N) ) / (< Y t N+2 )- Y( N- 1 ) )* ( Y< N + 1 )- Y < N YD I F F 
i-ll)*IY(N)-Y(N-in) 


40 


35 


1 

2 

3 

4 

5 

6 

7 

8 
9 


YDIFF 10 

BE TAIN) “« (DELTA fN) *( Yf N) -YCN— 1 ALPHA! N )*( Yt N + 2 )— Y (N ) ) ** 3 > / ( Y YD H- F 1 1 

i(N+l)-Y(Nn**3 YDI=F 12 

GAPMA(N) * -At PH A ( N t -BE T A ( N) - CELT A ( N ) YDIFF 13 

CONTINUE YDIFF 14 

DO 39 N*2tNVPl YDIFF 15 

SD(N) = (Y(N)-Y(N-i))/( (Y(N^l)-Y(N-l) )*(Y(N4i)“Y(N))) YDT=F 16 

YDIFF 17 
YDIFF 18 
YDIFF 19 
YDIFF 20 
YDIFF 21 
YDIFF 22 
YDIFF 13 
YDIFF 24 


SEtNT - i./( Y(N)-Y(N-l) »—!•'/€ YfN^ll-VfNI ) 

SF(N) = <Y(N+l)-Y<N) )/< <Y(N)-Y(N-i))*(Y(N + i)-Y( N-l)) ) 
CONTINUE 

C 2 = Y(3)*Y(4) /(Y(2)*(Y(3)-Y( 2) )*<Y(4)-Y( 2) ) ) 

C 3 = - Y ( 2 ) ♦ Y (4 IT (Y ( 3 ) *( Y ( 4 )-Y ( 3 ) ) ♦( Y( 3 )— Y ( 2 ) ) ) 

C4 = Y(2)*Y(31/(Y(4)*(Y(4|-Yl3n*(Y(41-Y(2n) 

RETUftN 

END 
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SUBROUTINE FLOE RIBC AP,XS I G,NS IG,U INF, ELO,Y,YMAX ) 


ELDER 

1 


0 1 ME NS I ON RCAPC100.3) ,XSIG(100» 


ELDER 

2 


BCAP (NSIG«-1 ,11=0. 


ELDER 

3 


XS=XSIG( 1 ) 


ELDER 

4 


XZ=XSI 3( NSIG+l t 


ELDER 

5 


IFIXZ-l.l 16,16,1 


ELDER 

6 

l 

OEADl -X2-X S 


ELDER 

7 


YMAX =1 * E— 1 0 


ELDER 

8 


SUH= •5*tXSIGt2)-X$) *RCAP(2,il 


ELDER 

9 


no 10 N*2 , NSI G 


ELDER 

10 


X = XSIG(N~4*i) 


ELDER 

11 


SUM=SUMf # 5*(X~X$IG(N))*(RCAP(N*i, l)«-BCAP< N» 1) ) 


ELDER 

12 


IF(N-NSIG) 4,2,4 


ELDER 

13 

2 

A NGLF *1 * 5708 


ELDER 

14 


GO TO 6 


ELDER 

15 

4 

ANGIE=ATAN( SQRT ( ( X-XS ) / ( XZ-X ) )) 


ELDER 

16 

6 

Y=SUM+BCAPU ,1 ) * (OEAOL* ANGLE— SORT I ( X-XSI *1 XZ-X 1 I ) 


ELDER 

17 


IFIY-YMAX1 10,10,8 


ELDER 

18 

8 

YMAX=Y 


ELDER 

19 

10 

CCNT1NUE 


ELDER 

20 


ELO»Y/YMAX 


ELDER 

21 


IF(ABS(ELD)-UINF) 20,20,12 


ELDER 

22 

12 

IFIFLO) 16,16,16 


ELDER 

23 

14 

ELO*-UlNF 


ELDER 

24 


GO TG 20 


ELDER 

25 

16 

FLD-U! NF 


ELDER 

26 

20 

CONTINUE 


ELDER 

27 


RETURN 


ELDER 

28 


END 


ELDER 

29 


12a 


8 

9 

10 


11 

12 


SUBROUTI NE RE vrt CJC , V, X , Y , MX , NY , RV ,[)R Y , U E , X5 . DEL5 » M ST, R EB) 

DIME NSI ON UCU00 ,3) ,V (100,2) , Y< 100 ) 

DIMENSION X(300> ,’IE (300 ,3) 

DIMENSION TARK74) ,TAB2(24),TAB3(24),TAB4124),X I T AB (24) 

DATA TA3 1 /24. 98 ,23.29,21.0*, 19.3 3,17.61, 15.29, 13.46, 11 .5'*, 10.3 

1. 38. 8. 35. 7. 3 2. 6. 29. 5. 31. 4. 4. 3. 5 7. 2. 22. 1. 26.. 66. .31, .14, .01, 3. ,C 
DATA TAD2 /20. 05 ,18 . 85 , 1 7 . 25 , IS . 04, 14 . 8, 1 3. 1 2, 1 1 . 77 , 10. 3, 9 . 36, C 

1 ,7.95,7.2,6.43,5.46,4.9 ,4.18,2.39, 1.8 6, 1.11, .62, .32, .04,0., .)./ 
DATA TAB 3 /16. 65, 15. 8, 14. 67, 13. 8, 12. 91, 11. 66, 10. 65, 8. 43, 8. 7 1,8. 
17.59,7.01,6.41 ,J. 7 7, 5. I 3 , 4. 5, 3. 3 1,2. 28, 1.48,. 9, .51, .C9, .01, )./ 
DATA TAB4 / 1 0. 1 2 , 10 . 05, 9 . 93 , 9 . 78 , 9 . 58 , 9. 1 7, 8 . 72 , 8 .0 8, 7. 6,7. >, 6. 

16.53.6. 18.5. 79.5.36.4.91.3.98.3.05.2.21. 1 .5. . 95. .22. .03.0./ 

DATA XI TAB / . 000 1 , . 0002 , .0005 , . 00 1, . 002, .005 , . 01 , . 02 , .0 3, .0 .0 

106.. 07,. 08,. 09, . I,. 12 ,. 14,. 16 ,. 18, .2, .25, . 3, .35/ 

FORMAT! ///40X,23HAT REATTACHMENT, BETA =613.5) 

M0UT=6 

RTR= SCR T (REPf) 

UC (l ,2) =0. 

UC(1,3)=0. 

V ( l , 1 I =0. 

V (1 ,2) =0. 
or 5 M=1 f MX 

IF7X5-XTMTT _ 4,T75' 

MST=M*2 
GO TO 6 
CONTINUE 
XA=X(MST-2) 

XB =X ( MST-1 ) 

OA=UE ( M5T-2 , 1 ) 

UB=UE (MST-1 ,1) 

ZA=ALOG 1 UA*DET5*REB I 

PGRAD=2.*< UA-UB)/ I I UA*U B ) * (X B-X A ) ) 

BE TM2=( . 0 9T4-'S75 ftT 7 5 E L5 * P G R AD ) ) / ( .0 249 *.004565*7 A ) 

IF (BFTM2-1. ) 8,7,7 

BFTM2^T7 

GO TO 10 

IFfBETM2-V31 979,17) 

BETM2 = . 3 

BFTA'=I .7 ( BE TM2 *BETM2 ) 

WRI TE ( MOOT ,3 ) BETA 

AGAM=. 09 74*Bt I M2-70249/"BETA 

BGAM=. 004565/BETA 

AH*1 ;-T5.3+3.9*BETM2T*r.0974-.0249*BETM2 ) 

BH=BETM2*( 5 .3*3 . 9*BET M2 ) *.004565 
G A M A = A (JAM — BG AM* C A 

DERI V=UA*RFB*EXP(-ZA)*GAMA*GAMA*( l.*BET A*( 1 . +AH *BH* Z A ) ) /( AH + BH< 

rzon 

ZB=ZA*DER1 V* ( XB-X A) 

DELB=E XPirB ) /TOT*RTTB) 

GAMB=AGAM-BG AM*ZB 

DE LL = . 3 5"* D"FLTT*RTR*TTET M2 / G AMB 

IF(DELL-Y( NY-3) ) 14,12,12 

~TTY=RY+DRY 

CALL YSET(RY,Y(2) ,NY,Y) 

corn 


RE ATT 

1 

REATT 

? 

R F A T T 

3 

RE ATT 

4 

9REATT 

5 

REAT T 

6 

i5RF ATT 

7 

RE ATT 

8 

, REATT 

9 

RFATT 

10 

f REATT 

n 

REATT 

12 

.RFATT 

13 

REATT 

14 

REATT 

15 

REATT 

16 

REATT 

17 

REATT 

18 

REATT 

19 

RFATT 

20 

REATT 

21 

REATT 

22 

REATT 

i- 

REATT 

24 

REATT 

25 

REATT 

26 

REATT 

27 

REATT 

28 

REATT 

29 

REATT 

30 

REATT 

31 

REATT 

32 

REATT 

33 

RFATT 

34 

REATT 

35 

REATT 

36 

REATT 

37 

REATT 

38 

REATT 

39 

REATT 

40 

REATT 

41 

REATT 

42 

REATT 

43 

REATT 

44 

REATT 

45 

l*RE AT T 

46 

REATT 

47 

REATT 

48 

REATT 

49 

REATT 

50 

REATT 

51 

REATT 

52 

REATT 

53 

REATT 

54 

REATT 

55 
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I F (BETA-4. ) 10 2,101, 101 
TERP9 = I.-4./BETA 
l NOE X=3 
GO TC no 

IF (RfTA-2/r 104,103, 103 
TERPB*.5*9ETA-l. 

I NOE X=2 
GO TO 110 
TFRPB=BE TA-1 • 

I NDE X=1 
K s 0 

TFRP 1 = 1 . -TERPB 
K*K+ 1 

GO TO (16,17 ,99) ,K 
G =G A MA 
OELTA*DE L5 
UEDGE^UA 

L=3 

GO TO 18 
G*GA WB 
OELTAVOELB 
UEOGE=UB 
1*2 

X ICH*G/ ( DE L TA*RTR*BETM2 ) 

UCOW=RTR* (UEDGE*G) **2 

EFC0-G/BETM2 

NLAMsNY 

DO 75 N=2 , N Y 

XI *Y ( N) * XI CO 

IF ( X I-* 3 5) 20,19,19 

UC(N,U=UEDGE 

GO TO 75 

CALL TERPF ( XI , I NDEX , T AB1 ,f AB2 , f AB3,T AB4, X ITAB, FP1 ) 

I NDP1 =1 NDE X+l 

CALL TERPF (XI , INOPi ,T AB1 ,T AB2 ,T AB3fTA84, X ITAB, FP2 ) 
F P=TE RPl *F P1+TERP3*FP2 

UC ( N ,U * UE OGE ♦ ( 1 . -E FCO*F P f 

IF(N-NLAM) 21,75,75 

AL TER=UCOW* Y (N I ' “ “ 

I F ( A L TE R-UC ( N , L ) ) 33 ,33,32 
UC (Ni t) *AL TfR 
qq TO 75 

NLAMxKi " ~ 

CONTINUE 

GO TO 50 “ 

DO 60 K*2 , 3 

SAVE 2=0. 

DO 6 0 N=3 , NY 
SAVEl=UC (N-l,Kl 

UC ( N— 1 ,K) = ($AVE2 + SAVEH>UC (N,K))/3« 

SA Vf2 = SA VlTI 

DUDX*0. 


00 65 N*2 , NY 


REATT 56 
REATT 57 
REATT 58 
REATT 59 
REATT 60 
REATT 61 
REATT 62 
REATT 63 
REATT 64 
REATT 65 
REATT 66 
REATT 67 
REATT 68 
RFATT 69 
REATT 70 
REATT 71 
REATT 72 
REATT 73 
REATT 74 
REATT 75 
REATT 76 
REATT 77 
REATT 78 
REATT 79 
REATT 80 
REATT 81 
REATT 82 
REATT 83 
REATT 84 
REATT 85 
REATT 86 
REATT 87 
REATT 88 
REATT 89 
REATT 90 
REATT 91 
REATT 92 
REATT 93 
REATT 94 
REATT 95 
REATT 96 
REATT 97 

REATT 98 
REATT 99 
REATTIOO 
REATT101 
RE AT T 102 
REATT103 
REATT104 
REATT105 
REATT106 
REATT107 
RTATTY08 
REATT109 
REAf TltO 
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V(N,U=VCN-1 ,U-<Y(N)-Y(N-1) ) *t OUDXP+OIJDX ) REATTUl 

V(N,2)=V(N,U REATT112 

65 OUOXOUDXP REATT113 

RETURN REATTlll 

END “ REATT115 
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SUBRDUTI NE EL PIT CAL PHI f ALPH?, EM I « T ORF f T HET l « U INF* OX I , CM PA* C MP AS ) ELPI T 


SAVE T*AL PHI ELPI T 

STFP*TORF*DXI EL P I T 

S I NS=SI N C S T E P) FLPI r 

cnss=cns(ST'EP) elpi r 

CONS T=2# *EMI*(UINF/TORF)**2 ELPI T 

ALPHl*THETZ+( ALPH1-THET Z ) *COS S+AL PH2*S INS/TO RF+CONST*( 2 .*CM PA-C MPAELP I T 
1S>*( I. -C OSS) ♦CONST*! CMP AS -CMP A) *($ INS -STEP*COS$ )/< TOR F* OXI) ELPI T 

ALPH2=ALPH2*COSS-TORF*SlNS*{SAVET-THETZ)*CONST*CCMPA-CMPASI*( U-COELPl T 
I s '9 »CON ST»C HPA»TCRF»$ I NS _ ELPI T 

RETURN ELPI T 

FN0 ELPI T 


1 

2 
3 
A 

5 

6 

7 

8 
9 

10 

11 

12 
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SUPRrUTINE VWA!THCfJARG.H,S,NVOR. X1,UINF,VZ IP.XGA'I.NGPi.nxn 

VWASH 

l 

DIMENSION VZIP(30) * XG.'M(V)) 

VWASH 

2 

DO 10 N=1 1 NG PI 

VWASH 

3 

01FF=XGAMIN»-X1 

VWASH 

4 

Sl!M=0. ' 

VWASH 

5 

DC 5 K = l ,NVOR 

VWASH 

6 

Sli Ms SU M + lit FF/ I DI FF*DIFF+H) 

VWASH 

7 

OIFF=DIFF-S 

VWASH 

8 

VZIPINMVZIP(NM-SIM*e.ARG 

VWASH 

9 

RETURN 

VWASH 

10 

END 

VWASH 

l 1 
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SUBROUTINE W ASH ( XGAM.NGAM, T I ME, AL PHI, ALPH2, HEAV E, AROT, F REQF ,PHI H 

*UWA S3 

1 


iINF,CMRR,NF,VZIPfHOTR,INDV) 

WASH 

2 


DIMENSION XGAMI30) , VZI P < 30 ) , C AMP R ( 24 ) 

WASH 

3 


NGPl = NGAM+1 

WASH 

4 


ANGLF = FREOF*TIMF 

WASH 

5 


GO TO (108,120) f I NOV 

WASH 

6 

1C8 

GO TO ( l 10,123) f MHT R 

WASH 

7 

L 10 

CCNST =-ALPH2*C0S(ANGLE)*UINF+HEAVE*C0S( ANGLE*PHIHI*ALPHl*UINF 

WASH 

8 


FACT =-ALPH2*FREQF*S!N( ANGLE ) *J INF 

WASH 

9 


GO TO 130 

WA SH 

10 

120 

CCNST=UI NF*ALPH1+HFAVE 

WASH 

11 


FAC T--UI NF *A LPH2 

WASH 

12 

13 3 

DO 10 M= 1 , NGPl 

WASH 

13 


X =XG AM ( M ) 

WASH 

14 


THETA = ARCT(X) 

WASH 

15 


5UM=0. 

WASH 

16 


CCUNT*0. 

WA SH 

17 


DO 20 N* 1 » NF 

WA SH 

18 


COUNT=COUNT+1. 

WASH 

19 

20 

SUM=SUM*C0UNT*CAM9R I N) *CCS I COUNT *THET A) 

WASH 

20 


IFtM-lI 2,4,2 

WASH 

21 

2 

IF(NGPl-M) 3,4,3 

WA SH 

22 

4 

SUM = SUM ♦ SUM 

WA SH 

23 


GO TO 50 

WA SH 

24 

3 

COUNT = 0. 

WA SH 

25 


COTT * X/SI N CTHET A I 

WASH 

26 


DO 30 N*1 , NF 

WASH 

27 


COUNT * COUNT+THETA 

WASH 

28 

30“ 

SUM=SUMi-COTT*CAMBR(N>*S INI COUNT) 

WA SH 

29 

50 

VZIPIM) * UINF*SUM«-CONST*EACT*I AROT-X 1 

WASH 

30 

10 

CONTINUE 

WA SH 

31 


RFTURN 

WA SH 

32 


END 

WASH 

33 
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APPENDIX B 


DETERMINATION OF COUPLING PARAMETERS 



APPENDIX B 


DETERMINATION OP COUPLING PARAMETERS 

The characteristic equation for the rotor blade is 

3 


B, 


2k X 


2k 


k=0 


= 0 


where 


B 0 “ f 0 ‘ M 


— 2 2 

w * T (e 

W M ee 


_ 2 2 


M 00 M ee 


Br 


f 2 4* 2 


M g e e 
M ee 


+ 2 


w g 2 %e T gfe 

M 00 M ee 


■g © 


V© 


M #j M ee M ©e 


B. 


f / M g e . “8 2 M efe> 2 

M e© 


M 0/J M ee 


M 


+ 2 


g © T g © + 2 M #e T #e 


M e© M jz$ M e© 


B< = 1 - 


M g © _ M 0© 2 

M ©©. ' % M e© 
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in which 


0 


_ 2 _ 2 _ 2 
“0 w e 


_ 2 _ 2 

0)0 OJ co 9 


2 2 _ 2 _ 2 

+ w ^ oj e 


_ 2 


0)0 + CO 0 + O) 


_ 2 

e 


The characteristic equation for the two-dimensional system 
is found to be 


£ - 

k=0 


2k 


x 2k - 0 


where 


D, 


2 2 
f - (j j h 3-n 

0 0 a 1 


- 2 2 
<°0 h b b 1 


_ 2 


D 2 - f e - £>0 g a x a x - 


_ 2 — 

0)0 S b x b i 


- h a a^ 2 - h b b 2 


Dj, = f;, " c 4 x 2 " ® a x a l g b x b l 


Dg = 1 - C 6 x 


-2 


in which 
h„ 


M 


00 


R M, 


©e 


h = 

b M ee 
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“ 2 h a fl l 


®b 2 h b A 2 


"/ h a A l 2 + S P 2 h b V 


C 6 - h a + h fc a 2 ‘ 


a l = A 1 < “0 2 1 b 1 + r m “/ 1 a 2 ) - B "/ l s . 


b ! - A 2 ( “0 2 1 s 1 + r m “/ X s > + B "/ l a . 


_ 2 


_ 2 


Equating Do/Dg to Bq/B^, D 2 /Dg to B 2 /Bg and Dji/Dg to Bii/Bg 
provides three relations in the three unknowns ° 

x t Isia. and 1 s 2 . If ai and bn are eliminated, the 
following equation for Y is obtained: 

( r l ^2 “ r 2 ^l) + ( r l s 2 ” r 2 s l )(^2 s ]l ' s 2 ) = 0 


r h g 2 ■ 

1 = - h a + — 7~ 

L 6 b J 


r o = Z&v - 1 h 


s 2 = ( - z>0 2 ) g a x , S;L = s 2 + 


2 h b Sa P 
g b^ x 


H - (1 - C 6 X ) B 2 /B 6 - fg + s 2 p + Jl 


hw P 


6 k x 
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2 


^2 _ _ ^6 X ~ ^c/ )/®6 " ^2 + w {3 F + ^ q / 

in which 

F = f 4 - B 4 /B 6 + (B 4 Cg/Bg - C 4 ) x 2 

With some algebraic manipulation, a polynomial of fourth 
degree in x^can be extracted from that equation. The value 
of x is taken to be the square root of the smallest positive 
root of that polynomial. The original equations are then 
used to solve for a^ and b]_ , from which 1 s -^ and lsg are 
readily obtained. 
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